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ABSTRACT 



We develop a new systematic approach to quantum field theory that is designed 
to lead to physical states that rapidly converge in an expansion in free-particle Fock- 
space sectors. To make this possible, we use light-front field theory to isolate vacuum 
effects, and we place a smooth cutoff on the Hamiltonian to force its free-state matrix 
elements to quickly decrease as the difference of the free masses of the states in- 
creases. The cutoff violates a number of physical principles of light-front field theory, 
including Lorentz covariance and gauge covariance. This means that the operators in 
the Hamiltonian are not required to respect these physical principles. However, by 
requiring the Hamiltonian to produce cutoff-independent physical quantities and by 
requiring it to respect the unviolated physical principles of the theory, we are able to 
derive recursion relations that define the Hamiltonian to all orders in perturbation 
theory in terms of the fundamental parameters of the field theory. 

We present two applications of this method. First we work in massless (f)^ theory 
in six dimensions. We derive the recursion relations that determine the Hamilto- 
nian and demonstrate how they are used by computing and analyzing some of its 
second- and third-order matrix elements. Then we apply our method to pure-glue 
quantum chromodynamics. After deriving the recursion relations for this theory, we 
approximate all physical states as two-gluon states (thus they are relatively simple 
single-glueball states), and use our recursion relations to calculate to second order 
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the part of the Hamiltonian that is required to compute the spectrum. We diag- 
onahze the Hamiltonian using basis-function expansions for the gluons' color, spin, 
and momentum degrees of freedom. We examine the sensitivity of our results to the 
cutoff and use them to analyze the nonperturbative scale dependence of the coupling. 
We investigate the effect of the dynamical rotational symmetry of light-front field 
theory on the rotational degeneracies of the spectrum and compare the spectrum to 
recent lattice results. Finally, we examine our wave functions and analyze the various 
sources of error in our calculation. 
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CHAPTER 1 



INTRODUCTION 



The Standard Model of particle physics represents much of our understanding of 
the fundamental forces of nature. This model is based on quantum field theory; so to 
test and improve the model, we must have methods for solving these theories. At the 
present time, the most successful and rigorous nonperturbative approach to solving 
realistic quantum field theories is the Euclidean lattice formulation. However, this 
method has its share of difficulties. Physical states tend to have very comphcated 
structures in the lattice approach. Multi-particle states and excited states are hard 
to simulate on the lattice, as are dynamic processes such as scattering. In addition, 
there are significant sources of error in lattice calculations, including those from finite 
lattice spacing, finite spacetime volume, and statistics. It is possible that some of 
the effects that are observed when field theories are simulated on a Euclidean lattice 
are not features of the field theories, but rather artifacts of the discretization of 
spacetime, the finite spacetime volume, or the unphysical Euclidean metric. Clearly 
other methods for deriving the states of quantum field theory from first principles are 
required, if not to compete with the lattice, then at least to test it. 

In this dissertation we develop a new systematic approach to quantum field theory 
that can lead to physical states that rapidly converge in an expansion in free-particle 
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Fock-space sectors (free sectors), and we present two applications of this approach. 
For the first apphcation, we demonstrate how the operators of a theory are calculated 
in our approach by calculating and analyzing some matrix elements of the invariant- 
mass operator^] (IMO) in massless 0^ theory in six dimensions. (We work in six 
dimensions so that the theory is asymptotically free.) Then we use our method to 
calculate part of the IMO in pure-glue quantum chromo dynamics (QCD), and we 
diagonalize it, obtaining glueball states and masses. These calculations appeared 
previously in Refs. [Q] and 0. The massless scalar calculation in Ref. [||] was 
extended to the massive case by Kylin, Allen, and Perry 0. 

The primary motivation for our approach is to force the physical states of a theory 
to have a relatively simple structure. We avoid some of the difficulties of the Euclidean 
lattice method by working in the continuum, using a Minkowskian metric, without 
truncating spacetime. Our formalism also has the advantage of being manifestly 
boost-invariant. Of course our method has its difficulties, but they are quite different 
from those of the lattice. Thus our method can complement the lattice approach by 
providing us with a very different view of the same field theories. 

In the most interesting field theories, it is unlikely that the physical states can 
quickly converge in an expansion in free sectors in an equal-time approach. This 
is because in equal-time field theory, the states often must be built on top of a 
complicated vacuum unless the volume of space is severely limited. For this reason, 
we work in light-front field theory (LFFT). In LFFT, it is possible to force the vacuum 

of a theory to be empty by removing from the theory all particles that have zero 

^The invariant-mass operator is given by the square of the momentum operator: V^V^ = A4^. 
See Chapter 2, Section 2 for more details. 
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longitudinal momentun]0. Any physical effects of these particles must be incorporated 
into the operators of the theory in order to obtain correct physical quantities. 

In LFFT, the Hamiltonian is trivially related to the IMO, and it is more natural 
to work with the IMO because it is manifestly boost-invariant. If the IMO satisfies 
three conditions in the basis of free-particle Fock-space states (free states), then its 
eigenstates will rapidly converge in an expansion in free sectors. First, the diagonal 
matrix elements of the IMO must be dominated by the free part of the IMO. Second, 
the off-diagonal matrix elements of the IMO must quickly decrease as the difference of 
the free masses of the states increases. If the IMO satisfies these first two conditions, 
then each of its eigenstates will be dominated by free-state components with free 
masses that are close to the mass of the eigenstate. The third condition on the IMO 
is that the free mass of a free state must quickly increase as the number of particles 
in the state increases. If the IMO satisfies all three conditions, then the number of 
particles in a free-state component that dominates an eigenstate will be limited from 
above. This means that the IMO's eigenstates will rapidly converge in an expansion 
in free sectors^. 

^This is because there are no negative longitudinal momenta and momentum conservation requires 
the three-momenta of the constituents of the vacuum to sum to zero. 

•^There are three subtleties here. The first subtlety is that the first and third conditions on the 
IMO will not be satisfied for those free states in which many of the particles have negligible center- 
of-mass transverse momentum and little or no mass. However, the contributions of these free states 
to the physical states in which we are interested are typically suppressed. For example, in QCD 
these free states have very large widths in transverse position space and are thus highly suppressed 
by confinement. In QED, the particles with negligible center-of-mass transverse momentum and no 
mass are long-transverse-wavelength photons. These photons decouple from the physical states in 
which we are typically interested, e.g. charge-singlet states like hydrogen and positronium. Thus the 
contributions to these physical states from the free states containing these photons are suppressed. 

The second subtlety is that exactly how quickly the IMO's off-diagonal matrix elements must 
decrease and the free mass of a free state must increase are not known. We assume that the rates 
that we are able to achieve are sufficient. This can be verified by diagonalizing the IMO and 
examining the rate of convergence of the free-sector expansion of its eigenstates. 

The third subtlety is that the coefficients of the expansion for highly excited eigenstates may grow 
for a number of free sectors and then peak before diminishing and becoming rapidly convergent. 
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To satisfy the first condition on the IMO, we assume that we can derive the IMO 
in perturbation theory. If this is true, then the couphngs are small and the diagonal 
matrix elements of the IMO are dominated by the free part of the IMO. To satisfy the 
second condition, we place a smooth cutoff on the IMO to force its matrix elements to 
quickly decrease as the difference of the free masses of the states increases. Once we 
have removed the particles with zero longitudinal momentum from the theory, it is 
reasonable to expect that the third and final condition will be satisfied automatically. 
This is because the free-particle dispersion relation of LFFT should force the free mass 
of a free state to quickly increase as the number of particles in the state increases. 
(See Appendix ^ for details.) 

By suppressing the matrix elements of the IMO that have large changes in free 
mass, the cutoff regulates the ultraviolet divergences of the theory. Unfortunately, it 
also violates a number of physical principles of LFFT, including Lorentz covariance 
and gauge covariance. This means that the operators in the IMO are not required 
to respect these physical principles, and renormalization is not simply a matter of 
adjusting a few canonical parameters. The simplest way to systematically determine 
the IMO in this case is in perturbation theory. In order for a perturbative computation 
of the IMO to be strictly valid, the theory must be asymptotically freeQ. If this is 
the case, then by requiring the IMO to produce cutoff-independent physical quantities 
and by requiring it to respect the unviolated physical principles of the theory^, we can 

derive recursion relations that define the IMO to all orders in perturbation theory in 

^Our method may work even if the theory is not asymptoticaUy free. For example, it works in 
QED because the scale at which the electron charge is large is astronomical. 

^Some of the physical principles, such as cluster decomposition, are violated in a very specific 
manner and can still be used to restrict the form of the IMO. However, the restriction in this situation 
is always weaker than it would have been had there been no violation of the principle. 
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terms of the fundamental parameters of the field theory. If our cutoff is large enough, 
then the couplings will be small and the perturbative approximation to the IMO may 
work well. 

The physical principles that we use to determine the IMO form a subset of the full 
set of physical principles of light-front field theory. This raises the question of how 
the remaining principles, which are violated by our cutoff, are restored in physical 
quantities. Since the IMO is uniquely determined by the principles that we use, the 
remaining principles must be automatically respected by physical quantities derived 
from our IMO, at least perturbatively. If they are not, then they contradict the 
principles that we use and no consistent theory can be built upon the complete set 
of principles. The reason that this process is possible is that there are redundancies 
among the various physical principles. 

It is possible to compute operators other than the IMO in our approach. Although 
we compute operators perturbatively, we can use these operators to compute nonper- 
turbative quantities. For example, the spectrum can be computed by diagonalizing 
the IMO (see Chapter 5, Section 4). However, there are drawbacks to computing 
operators perturbatively. It is possible that there are intrinsically nonperturbative 
effects in the theory that require nonperturbative renormalization. Any such effects 
are neglected in this approach. Another problem is that perturbative renormalization 
makes nonperturbative physical quantities somewhat cutoff-dependent. 

In general, field theories have an infinite number of degrees of freedom. How- 
ever, since our IMO will cause the physical states of the theory to rapidly converge 
in an expansion in free sectors, we can truncate this expansion. This means that 
approximate computations of physical quantities will require only a finite number of 
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finite-body matrix elements of operators. In addition, since we assume tliat we can 
compute tliese matrix elements perturbatively, we do not implement particle creation 
and annihilation nonperturbatively in these calculations. This allows us to always 
work with a finite number of degrees of freedom. 

When we perturbatively calculate matrix elements of operators, we do not trun- 
cate the space of intermediate states that can appear; so our approach does not 
use a Tamm-Dancoff truncation^ 0,11,^1. We also do not completely eliminate any 
interactions, such as those that change particle number. These strengths of our ap- 
proach allow us to better describe physical theories. However, the truncation of the 
free-sector expansion of physical states has drawbacks that are similar to those of 
perturbative renormalization. It neglects any physical effects that require an infi- 
nite number of particles and contributes to the cutoff dependence of nonperturbative 
physical quantities. 

The accuracy of our results and the strength of the cutoff dependence of our non- 
perturbative physical quantities are determined by the order in perturbation theory 
to which we calculate the operators of the theory and the number of free sectors that 
we keep in the expansion of physical states. If we use a cutoff that is too small, then 
the couplings of the theory will be large, and it will be necessary for us to keep many 
terms in the expansion of the operators. If we use a cutoff that is too large, then 
the free-sector expansion of the physical states will converge slowly, and it will be 

necessary for us to keep many sectors in the states. We assume that if the order of 

^The only degrees of freedom that we remove from the full theory are the particles with zero 
longitudinal momentum. We should be able to replace the physical effects of these particles with 
interactions without compromising the validity of the theory. This is because these particles are vac- 
uum effects or have infinite kinetic energies (or both) , and thus are not observable in the laboratory 
as particles. 
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perturbation theory and the number of free sectors are manageable, then there is a 
range of cutoff values for which the approximations work well and physical quantities 
are relatively accurate and cutoff- independent. 

As we mentioned, we remove from the theory all particles with zero longitudinal 
momentum. We should replace their physical effects with interactions. However, 
due to the limitations of our method, we can reproduce only those effects of these 
particles that can be derived with perturbative renormalization and a small number 
of particles. 

The field theory to which we are most interested in applying our method is QCD, 
the fundamental theory of the strong nuclear interactions. There are a number of 
approaches that are similar to ours [0, |[ |, 0, [12|, |T3|, and some of these 



methods have been used to calculate the physical states of QCD [0, |T5[. These 
calculations are based on nonrelativistic approximations and use sharp step-function 
cutoffs. Nonrelativistic approximations drastically simplify the diagonalization of the 
IMO, but are insufficient for states containing light quarks or gluons. Sharp cutoffs 
prevent the complete cancellation of the infrared divergences that appear in light-front 
gauge theoriesQ 

Our approach is completely relativistic and uses smooth cutoffs to ensure the 
complete cancellation of the light-front infrared divergences. It is largely based on 
the renormalization methods of Perry Perry and Wilson |17], Wilson |18], and 
Glazek and Wilson [0 , as well as the Hamiltonian-diagonalization methods of Wegner 

m 



^These appear due to the exchange of massless gauge particles with arbitrarily small longitudinal 
momentum. 
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To the best of the author's knowledge, this dissertation contains a number of 
original contributions: 

• a method for deriving recursion relations that uniquely determine the IMO (or 
equivalently, the Hamiltonian) to all orders in perturbation theory in a light- 
front field theory; 

• derivations of these recursion relations in massless 4>^ theory and pure-glue QCD; 

• a computation of part of the emission matrix element of the renormalized IMO 
to third order in light-front massless 0^ theory; 

• a computation of the two-gluon to two-gluon matrix element of the renormalized 
IMO to second order in light-front pure-glue QCD; 

• a computation of the pure-glue QCD glueball spectrum in a continuum light- 
front approach. (This is the first continuum relativistic light-front QCD bound- 
state calculation.) 

This dissertation is organized as follows. Each of the chapters and appendices 
begins with a brief outline of the material that it contains. In Chapter 2 we introduce 
LFFT and derive the canonical Hamiltonians for massless theory in six dimensions 
and pure-glue QCD. In Chapter 3 we lay the groundwork for our method and derive 
the recursion relations that determine the IMO in each of the two theories. In Chapter 
4 we calculate and analyze some of the IMO's second- and third-order matrix elements 
in massless 0^ theory. In Chapter 5 we calculate part of the pure- glue QCD IMO, and 
diagonalize it, obtaining glueball states and masses. We conclude with a summary 
and a discussion of the direction of future work in Chapter 6. In Appendix A we 
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give a qualitative argument that shows that in LFFT, the free mass of a free state 
should rapidly increase with the number of particles in the state. In Appendix B we 
prove that the transformation that we use to alter the cutoff on our IMO is unitary. 
In Appendix C we discuss some of the technical issues involved in the numerical 
calculation of the matrix elements that enter our calculation of physical states in 
pure-glue QCD. 
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CHAPTER 2 



LIGHT-FRONT FIELD THEORY AND CANONICAL 

HAMILTONIANS 



In this chapter we begin with an introduction to hght-front coordinates and a 
discussion of some of the special properties of hght-front momenta. We then discuss 
the nature of the Poincare symmetries in LFFT. We conclude by deriving the canon- 
ical Hamiltonians for massless 0^ theory in six dimensions and pure- glue QCD. The 
reader who is interested in a historical survey of LFFT or a more detailed analysis of 
its properties is encouraged to consult the recent reviews and introductions [ED . 



2.1 Light-Front Field Theory 
2.1.1 Light-Front Coordinates 

In 1949, Dirac introduced the Hamiltonian light-front approach to dynamics as 
part of his effort to construct a viable Hamiltonian formalism for relativistic systems 



21| . In light-front dynamics in d spacetime dimensions, any (i- vector a is written 

a = (a"*", a~, d±), (2.1) 
where in terms of equal-time vector components, 

d-2 d-2 

d± = ^i^i = (2.2) 

i=l i=l 
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and Cj is the unit vector pointing along the i-axi^. Since a spacetime coordinate is a 



d-vector, it is also written this way: 



X = {x^, X , x±). 



(2.3) 



Rewriting the usual rf-vectors in this way does not do anything by itself. What 
defines light-front dynamics is the choice of x+ as the "time" coordinate, x~ as the 
"longitudinal space" coordinate, and x± as the "transverse space" coordinate. This 
choice means that when we quantize fields in light-front dynamics, we quantize them 
at equal light-front time (x^) rather than at equal time This choice also implies 
that the light-front Hamiltonian is the operator that evolves states forward in light- 
front time. We will see that these differences have profound effects on the dynamics 
of LFFT. 

If a LFFT is to yield the same results as the corresponding equal-time field theory, 
it is necessary for Lorentz scalars to agree in the two approaches. Since the inner 
product in the standard equal-time approach is 



d-l 



a-b = - a'b\ 



(2.4) 



i=l 



we find that in light-front coordinates. 



a ■ b = -a^b + -a b^ — d±- b±. 



(2.5) 



This means that the nonzero components of the metric are 
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= 9 



+ 



9 




(2.6) 



Another common convention is a 



± 



(aO±a'^-i)/V2. 
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and 



9+- 
9ij 



9-^ = 2> 



-5. 



(2.7) 



where j = 1, 2, 3, . . . , d — 2. The d-dimensional volume element is 





dx° 


dx+ 


dx~ 






dx+ 


dx~ 



dx~^dx d'^ "^Xx 



-dx dx d x±, 
2 



(2.8) 



where the double bars indicate the absolute value of the Jacobian. The gradient 
operator is written just like any other rf- vector: 



(2.9) 



where 



d 



dx^ 



(2.10) 



and 



d 



dx\ 



(2.11) 



According to Eq. ( p.5| ), the inner product of a momentum and a position is 

1 ^ _ 1 



p ■ X = -p X + X ~ p± ■ x±. 



(2.12) 



This implies that the momentum conjugate to x~ is p'^, the momentum conjugate 
to x^ is p~, and the momentum conjugate to x± is p±. Since x~ is the longitudinal 
position, we refer to p'^ as the longitudinal momentum, and since x"*" is the time 
coordinate, p~ is the energy. Similarly, since x± is the transverse position, we refer 
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to pi as the transverse momentum. For a free particle, — m^. This means that 
the kinetic energy of a particle is 



2.1.2 Special Properties of Light-Front Momenta 

The behavior of light- front momenta is the source of many of the special properties 



we see that every particle has p+ > 0. Since the vacuum has no momentum, and there 
are no particles with p+ < 0, the vacuum can only contain particles with p+ = 0. 
Thus we can remove the effects of the vacuum from the physical states of the theory 
by removing from the theory all particles with p"*" = 0. This is the first step in our 
effort to construct physical states that have a relatively simple structure. 

Particles with = are either vacuum effects or have infinite kinetic energies 
(or both). These particles are not observable in the laboratory; so we should be able 
to replace their physical effects by modifying the operators of the theory, without 
compromising the validity of the theory. In principle, we should be able to modify 
the operators of the theory to include the physical effects of the infinite-kinetic-energy 

— particles using renormalization. However, there remain two open questions. 
First, can this renormalization be done perturbatively, with a finite number of par- 
ticles that have p"*" > 0? Second, what about the physical effects of the p"*" = 
particles that do not have large kinetic energies? It seems unlikely that perturbative 
renormalization will yield operators that include these effects. In this dissertation, 
we use perturbative renormalization with a finite number of particles and leave any 



of LFFT. From 




(2.14) 
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nonperturbative/infinite-particle-number renormalization effects for future consider- 
ation. 

Tlie dispersion relation shows that a particle's kinetic energy behaves differently 
under the scaling of transverse and longitudinal momenta. This means that these 
momenta must be assigned different scaling dimensions [Q. The dispersion relation 
also shows that longitudinal locality is violated in LFFT (since the longitudinal mo- 
mentum appears in the denominator), and thus cannot be used to restrict the form 
of the operators of the theory. We will see that our cutoff violates transverse locality, 
but since this violation is weak, we will be able to use an approximation to transverse 
locality to restrict the form of the operators. 

2.1.3 Poincare Symmetries in Light-Pront Field Theory 

The choice of light-front quantization has interesting consequences for the Poincare 
symmetries in LFFT. To analyze these consequences, we compare the canonical def- 
initions of the generators of these symmetries in the equal-time and light-front ap- 
proaches. We consider the case of scalar field theory for simplicity and the case of 
d = 4 dimensions for physical realism. 

The energy-momentum tensor is defined in terms of the Lagrangian C by 

dC 

In equal-time field theory, the generator of translations along the yU-axis is (/i = 
0,1,2,3) 

= j d^xT^^". (2.16) 

is the Hamiltonian. The angular-momentum tensor is 
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Ji'^ = J d^x [x'^T°'^ - x^'T'^"] , (2.17) 
in terms of which the generator of boosts along the i-axis is {i = 1,2,3) 

/C, = (2.18) 

and the generators of rotations are 
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Jl = J 



J2 = J'\ 



J3 = J', (2.19) 

where J7i generates rotations of the system around the i-axis. 

Notice that the only component of the energy-momentum tensor that appears 
in J''^'^ and contains the interaction (which is in C) is T^^. This means that the 
Hamiltonian and the boost generators move states off the quantization surface = 0. 
For this reason, they are called dynamic operators. The rotation generators do not 
move states off the quantization surface and are said to be kinematic operators. 

In LFFT, the generator of translations along the /i-axis is^ (/i = +, — , 1, 2) 

= ^J d^x^dx-T+^'. (2.20) 
The angular-momentum tensor is 

Ji'^ = d^x±dx- [x^T^" - x'^T+^J . (2.21) 

^Although we speak of translations along the longitudinal axis , longitudinal position is not an 
observable quantity in light-front dynamics and wave packets cannot be localized in x~ [p2[ . This 
is related to the fact that > 0. 
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The Lorentz generators that are contained within J^" are divided into the hght-front 
boost generators: 

Bi = J+\ 

/C3 = (2.22) 



and the hght-front rotation generators: 

^1 = J~\ 
T2 = 

Js = J'^- (2.23) 

Bi and B2 are the generators of boosts along the 1-axis and the 2-axis respectively. 
Under these transverse boosts, each particle's transverse momentum (whether the 
particle is composite or point-like) transforms according to 

P±^P±+P'^v±, (2.24) 



where v± is a boost parameter [22, 23] . /C3 is the generator of boosts along the 3- 
axis, just as in equal-time field theory. In LFFT this boost is called a longitudinal 
boost. Under a longitudinal boost, each particle's longitudinal momentum transforms 
according to 

p+ ^ eV, (2.25) 

where u is a boost parameter. 

In LFFT, the only component of the energy-momentum tensor that appears in 
J'^'^ and contains the interaction is T"*'~. This means that the transverse boost 
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generators are kinematic. The longitudinal boost generator is dynamic in general. 
However, under a boost along the 3-axis, the light-front time variable is transformed 
according to 

x+ x^er (2.26) 

This means that the longitudinal boost generator shifts states off the quantization 
surface x~^ = constant only if the constant is nonzero. To force /C3 to be kinematic, 
we use x"*" = as our quantization surface. 

jFi and JF2 are dynamic operators and generate rotations of the quantization sur- 
face x+ = around the surface of the light cone = jTs is the generator of 
rotations about the 3-axis and is kinematic, just as in equal-time field theory. 

The differences in the Lorentz generators in equal-time and light-front field the- 
ories have important implications for the construction and classification of physical 
states. In general, we would like to solve for states as functions of their momenta 
and classify them according to their mass and spin, as well as their transformation 
properties under parity and charge conjugation. In equal-time field theory, it is dif- 
ficult to solve for states as functions of their momenta because the boost generators 
are dynamic. On the other hand, classifying states is straightforward because the 
rotation generators, the parity operator, and the charge conjugation operator are all 
kinematic. 

In LFFT with the quantization surface x+ = 0, it is simple to solve for states as 
functions of their momenta since boosts are kinematic. However, classifying states 
is very difficult because rotations and parity are dynamic. (Parity interchanges the 
light-front time and longitudinal-space coordinates.) It is straightforward to partially 
classify states in LFFT since charge conjugation and rotations about the 3-axis are 
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kinematic, but for a complete classification, we must either compare theoretical re- 
sults to experimental data or diagonalize the rotation generators and parity operator 
in addition to the IMO to compute their simultaneous eigenstates. Comparing our 
theoretical results to experimental data is tricky because when we make approxima- 
tions, such as renormalizing perturbatively or truncating the free-sector expansion 
of states, we disrupt the properties of the states under dynamic symmetry transfor- 
mations. On the other hand, computing simultaneous eigenstates of the IMO, the 
rotation generators, and the parity operator is expected to be much more difficult 
than diagonalizing the IMO alone. 

2.2 The Canonical Hamiltonians 

In this section, we simultaneously derive the canonical Hamiltonians for massless 
(f)^ theory and pure-glue QCD. For massless 0^ theory we work in d — 6 dimensions 
so that the theory is asymptotically free, and for pure-glue QCD we work in d = 4 
dimensions. 

The canonical Lagrangian density for massless scalar field theory with a three- 
point interaction is 




-90 • 90 - 



(2.27) 



The canonical Lagrangian density for pure-glue QCD is 




^ cui/^ c 1 



(2.28) 



where 



(2.29) 
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Greek indices are Lorentz indices, c's are color indices, repeated indices are summed 
over, and the /'s are the SU(iVc) structure constants. For each theory, we derive the 
canonical Hamiltonian from £ by the following procedure: 

1. In pure-glue QCD, we choose the hght-cone gauge, A+ = 0, and derive the 
Euler-Lagrange equation that determines A~ in terms of A±c- 

2. Using the canonical procedure and treating the field classically (i.e. letting 
it and its derivatives commute), we derive the Hamiltonian in terms of in 
massless cf)^ theory and in terms of A±c in pure-glue QCD, dropping terms that 
are zero if the field is zero at spacetime infinity. 

3. We quantize the field by expanding it in terms of free-particle creation and 
annihilation operators. (We define the field expansions below.) 

4. In each term in the Hamiltonian, we treat the creation and annihilation opera- 
tors as if they commute and move all the creation operators to the left of all the 
annihilation operators. (This "normal ordering" drops some constants in both 
theories and the so-called "self-inertias" in pure-glue QCD.) 

5. We drop the terms in the Hamiltonian that have no effect if there are no particles 
with p"*" = 0. 

We work in the Schrodinger representation, where operators are time-independent 
and states are time-dependent. Thus we quantize the field by defining it to be a su- 
perposition of solutions to the Klein-Gordon equation (since the particles are bosons), 
with the quantization surface — 0. In massless 0^ theory, we use 

<j>{x-,x±) = [ oAaie-'P'-'' + aWP'-^] , (2.30) 

J L ^ x+=0 
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and in pure-glue QCD we use 



ai£_Lsie 



—ipi-x 



(2.31) 



where 



(2.32) 



Ci=l Si=-l,l ''^dPi 

and Ud = 2{2nY~^. Here V is the d-momentum operator, e is a positive infinitesimal, 
Si is the spin polarization of particle i, and the gluon polarization vector is defined 

by 



e±s = 



(2.33) 



In massless cj)^ theory, e = and spin and color degrees of freedom are ignored in 
all quantities (such as Di). The purpose of e is to regulate divergent effects from 
exchanged gluons (either instantaneous or real) with infinitesimal longitudinal mo- 
mentum. We take e — > before we calculate physical quantities (see Chapter 5). 

In the course of deriving the canonical Hamiltonian in pure-glue QCD, it is nec- 
essary to take the inverse longitudinal derivative of the gluon field. We do this by 
moving the derivative inside the field expansion: 

1 ^ „ 1 



Cl 



—ipi-x 



+ 



(2.34) 



x+=0 



ipi ipi 

The expansion coefficients a\ and a, are identified as particle creation and anni- 
hilation operators. They follow the convention 



Ctj (p j , S j , Cj ) , 



(2.35) 



and have the commutation relations 
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[ai,a]] = 5ij (2.36) 

and 

[ai,a,] = [aj,aj] = 0, (2.37) 

where 

= napt^^"-'^ (pi - Pj)Ss.,sA.,c, (2.38) 

and 

S^'-'\pi-pj) = S{pi-pl)S^'-'\p,^-pj^). (2.39) 

Let be the IMO. Since the momentum conjugate to x+ is p~ , the Hamiltonian 
is identified as P~, and it follows from 

= (2.40) 

that 

The Hamiltonian can be written as the sum of a free part and an interacting part: 

V-^Vf,^ + v, (2.42) 

where 

^f7ee = jDi^aia,, (2.43) 

and in massless 0^ theory 
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^-9^ [4aia25('^-^)(p3 - Pi -P2) +44«i<^^''~'nP2 +P3 -Pi)]- (2.44) 
In pure-glue QCD, we give the interaction in terms of its "modified" matrix elements: 
V = jDiD2Dsalalaid^'^-^\pi-p2-p3) ^g2g3\v\gi:\> 

f 

+ g^ jD^D2D^alaia2 5^'^-^\pi+p2-p^) ic/3k|c/i52^ 

4 

+ 5'^^ jDiD2D^Diala\aia2 5^'^~'^\pi +P2 -Pa -P4) ^ ifi'35'4kl5'i5'2^i 

4 

+ D1D2D3D4 444ai - P2 - P3 - P4) XI i5'25'35'4kl5'i:|>i 

4 

+ [D,D2DsD,alaia2asS^''-'\p 

1 +P2 +P3 ~ P4) X/ *l^fl'4klfl'l5'25'3^i) 
i=l 



(2.45) 



where a modified matrix element is defined by 

<i\v\j> 



ndS('^-^){pi-pj) 



(2.46) 



9=1 



The modified matrix elements are 



^52^3^151^ = if 



C1C2C3 



-Ss2,S3^±si - |(P2± -P3±) - ^(P2 -P^) 

P2± , 



+ ^SUS3^L2- |(Pl± +P3±) - -+^Pt+Pt)j 

+ <^si,.2Cls3- |-(^1^ +P2±) + ^(P| +P2 ) 



C1C2C3 



(^.2,s3£±sr |(P3± +P2±) - ^(Ps" +P2") 

P3± , 



+ (^5l,S3£±S2- |-(Pl± +P3±) + ^(Pl^ +P3 ) 



(2.47) 



and 
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and 



and 



, f C1C4C fC3C2C / c c _ K _ K _•) 

~ J J l"S2,S4"Sl,S3 "-'si,S4"-'S2,S3;j ) 

<[:5354kbl^2:f2 = /ciC3C^C4C2C^^ ^ ^ + Pl + Pt 

Pi -PiPi - P2 
Pi ~ Pa Pi ~ P2 

^9Mv\9i92^A = r'''''r'"'%,,,As,^ /l 7 ^{ 7 ^ (2-48) 

Pi +P2 Pa +P3 



I fClC4CfC3C2C(_! f _L X _ X ^ 

~ J J \ "S2,S4"S1,S3 ~ "S3,S4"S1,S2^ 

+ 5si ,54*^^2,^3 )] 

i525354|t'|5i^2 = -r'''J''''-^%,,s25,3,s,, 
^5253541^^15^3 = -r^^7'^'^'<^.i,.3<^.- 



+ pt pt 


-pt 


Pi - pt pt 


+ pt 


Pi + pt pt 


-pt 


Pi - pt pt 


+ pt 


Pi + pt pt 


-pt 



4:5253541^15^4 = -r"7""^^...4'^.3,.2 ^l yj^/r^ ^ (2.49) 

Pi — Pa P2 "T P3 



^54k|5l5253:fl = [r '^7'^'^' (-<^.2,.-3'^.l,.4 + <^.3,.4'^.l,.-2) 

I fClC2CfC3C4C X _ <X S -)] 

— J J \ "S2,S4"S1,S3 ~ "S1,S4"S2,S3/J ' 

4:54^1515253:^2 = -/""7"-"^...4^S2,.3 l ^ 5 f+ 7 , 

Pi — Pa Ps + P2 

^54^1515253^3 = -r-7--^5.2..45s...34^4^, 

P2 ~ P4 P3 + Pi 
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^9M9i929:^>^ = -/^3C4c^c.c.c^ ^ P3 +P4 Pi ^2.50) 

Vt -Ptvi + P2 

where s„ = — s„. 

Let l^i) = in massless 0^ theory and = \gi) in pure-glue QCD. We work 
in the free basis, the basis of eigenstates of V^^^. They are given by 

= «I4---4lo), (2.51) 

for any integer n > 0. The associated eigenvalue equation is 

Koo 166 ■ ■ ■ 6) = E p7 166 ■ ■ ■ 6) , (2.52) 

i=l 

where 

p- = 4 (2.53) 

Pi 

(since = 0) and the sum is zero if n = 0. 

The noninteracting limit of the operator dispersion relation, Eq. (|2.41| ), is 

= (2.54) 

where A^lee the free IMO. It has the eigenvalue equation 

•ML 166 ■ ■ -6) = |ei6 ■ ■ - 6) , (2.55) 

where 

n 

M'' = P+Y.pT -PI (2.56) 

and P is the total momentum of the state. 

Finally, in terms of the free states, the completeness relation is 

1 = |0) (0| + j D, 16) (61 + ^ / D,D^ |ei6) (661 + ■ ■ ■ . (2.57) 
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CHAPTER 3 



THE METHOD FOR COMPUTING FREE-STATE 
MATRIX ELEMENTS OF THE INVARIANT-MASS 

OPERATOR 

In this chapter we lay the groundwork for our method for computing the IMO 
of a field theory, and we derive the recursion relations that uniquely determine the 
IMO in each of the theories that we consider. After defining the action of our cutoff, 
we place a number of restrictions on the IMO. First we force the IMO at a given 
cutoff to be unitarily equivalent to itself at a higher cutoff. This implies that the 
IMO is unitarily equivalent to itself at an infinite cutoff, and will therefore yield 
cutoff-independent physical quantities. From the statement of unitary equivalence, 
we develop a perturbative series that relates the interactions at two different cutoffs. 
We then proceed to use physical principles to restrict the form of the IMO. We require 
it to conserve momentum and to be invariant under boosts and rotations about the 
3-axis. Although our cutoff violates exact transverse locality, we are able to require 
the IMO to respect an approximate transverse locality. In practice this means that 
the IMO's matrix elements are analytic functions of transverse momenta. The cutoff 
also violates cluster decomposition, but we show that the implications of this violation 
are simple enough that we can still use this principle to restrict the form of the IMO. 
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We show that in each of the theories that we consider that the IMO must become 
the free IMO in the noninteracting hmit, and it must be a function only of the cutoff 
and the couphng. The final physical restriction on the IMO is that it must reproduce 
the perturbative scattering amplitudes of the appropriate theory. This restriction 
specifies the form of the low-order interactions in each theory that we consider. 

We use the restrictions from physical principles and the perturbative series that 
relates the interactions at two different cutoffs to derive the recursion relations that 
determine the IMO. A crucial step in this process is the removal of the coupling from 
the perturbative series. This allows us to separate the cutoff dependences of the 
operators in the interaction from the cutoff dependences of their couplings. 

3.1 The Cutoff 

The IMO is a function of the cutoff, and can be spht into the canonical free IMO 
and an interaction: 

M'iA) = ML, + MUA). (3.1) 

At this point the interaction is undefined. To force the off-diagonal matrix elements 
of the IMO to decrease quickly for increasingly different free masses, we regulate it 
by suppressing its matrix elements between states that differ in free mass by more 
than the cutoff: 

{F\M\A)\I) = {F\MUI) + {F\M^^,{A)\I) 

= M'^{F\I)+e-^{F\V{A)\I), (3.2) 

where \F) and |/) are eigenstates of the free IMO with eigenvalues Mj, and M|, and 
Apj is the difference of these eigenvalues: 
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A^, = - Ml (3.3) 

V{A) is the interaction with the Gaussian factor removed, and we refer to it as 
the "reduced interaction." To determine the IMO, we must determine the reduced 
interaction. 

We will see that {F\ V{A) \I) does not grow exponentially as A'jpj gets large; so 
the exponential in Eq. ( |3.2| ) forces the off-diagonal matrix elements of the IMO to 
rapidly diminish as A'j^j grows. This satisfies the second of our three conditions on 
the IMO and regulates it. 

3.2 The Restriction to Produce Cutoff-Independent Physical 
Quantities 

Unfortunately, any regulator in Hamiltonian LFFT breaks Lorentz covariance, and 
in gauge theories it also breaks gauge covarianceQ. This means that the operators in 
the IMO are not required to respect these physical principles. Instead, the IMO must 
be allowed to contain all operators that are consistent with the unviolated physical 
principles of the theory. In addition, since there is no locality in the longitudinal 
direction in LFFT, these operators can contain arbitrary functions of longitudinal 
momenta. 

To uniquely determine the operators that the IMO can contain, as well as the 
coefficients of these operators, we place a number of restrictions on the IMO. The 
first restriction is that it has to produce cutoff-independent physical quantities. We 

impose this restriction by requiring the IMO to satisfy 

^°Our regulator breaks these symmetries because the mass of a free state is neither gauge-invariant 
nor rotationally invariant (except for rotations about the 3-axis). 
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M^{A) = U{A, cA) M^{cA) U\A, cA), 



(3.4) 



where c > 1 and is otherwise arbitrary, and ?7 is a unitary transformation that 
changes the IMO's cutoff. Note that we are considering A4'^{A) to be a function of 
its argument; i.e. A4^{cA) has the same functional dependence on cA that 7V1^(A) 
has on A. 

To see that Eq. ( p.4| ) imphes that the IMO will produce cutoff-independent phys- 
ical quantities, note that the IMO on the right-hand side (RHS) of the equation can 
be replaced by iterating the equation: 

M\A) = U{A, cA) M\cA) U\A, cA) 
f/(A,cA) U{cA,c^A)] M\c^A) [u\cA,c^A) f/^(A,cA) 
f/(A, cA)U{cA, ^A)U{^A, ^A)\M^{e'A) [u\^A, ^A)U\cA, ^A)U\A, cA) 

(3.5) 



Therefore, since c > 1, Eq. ( ^.41 ) implies that A^^(A) is unitarily equivalent to 
limA^oo -^^(A). This means that A^^(A) will produce cutoff-independent physical 
quantitiesj^. 

For simplicity, we define 



A' = cA, 



(3.6) 



and then Eq. 



takes the form 



^^If the couplings in Al^(A) diverge as A — > oo, then this argument cannot be used in perturbation 
theory. This is why we assume that the theory is asymptotically free. On the other hand, if the 
couplings are small until the cutoff is astronomical, as in QED, then we can use perturbation theory 
to prove that the IMO at a cutoff of physical relevance is unitarily equivalent to itself at a huge 
cutoff, and will thus produce physical quantities that are almost cutoff- independent. 
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M\A) = U{A, A') M\A') U\A, A'). (3.7) 

To calculate the IMO perturbatively, we need to implement this equation per- 
turbatively. In our approach, the cutoff that regulates an IMO also serves to define 
the scale for the couplings that appear in the IMO; so the couplings in 7V1^(A) are 
defined at the scale A. To use Eq. ( |3.7D perturbatively, we have to perturbatively 
relate the couplings at the scale A to the couplings at the scale A'. This perturbative 



relationship is well-defined only if A' is not very large compared to A p4|. To fulfill 



this requirement and Eq. ( p.6|) (recall that c > 1), we choose A' to satisfy 

A < A' < 2A. (3.8) 

To make Eq. ( |3.7| ) a complete statement, we have to define the unitary transfor- 
mation. The transformation is designed to alter the cutoff implemented in Eq. ( |3.2| ), 
and is a simplified version of a transformation introduced by Wegner [l^, modified 



for implementation with the IMO. The transformation is uniquely defined by a linear 
first-order differential equation: 

dU{A,A') 



d(A-4) 

with one boundary condition: 



r(A)f/(A,A'), (3.9) 



[/(A,A) = 1. (3.10) 

We show in Appendix ^ that U (A, A') is unitary as long as T(A) is anti-Hermitian 
and linear. We define 

T(A) = [A^Le,-^'(A)], (3.11) 
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which is anti-Hermitian and hnear. The transformation introduced by Wegner uses 



T(A)= MUA),M\A) 



(3.12) 



where A4^{A) is the part of A4'^{A) that is diagonal in the free basis. (It includes the 
diagonal parts of all interactions.) Our transformation regulates changes in eigenval- 
ues of A^free) ^^^1 Wcgucr's transformation regulates changes in eigenvalues of /^"^{A). 

To solve for 7V1^(A) perturbatively, we need to turn Eq. ( ^.7] ) into a perturbative 
restriction on the reduced interaction, V{A). We begin by taking the derivative of 
Eq. (P): 

dM\A) 



d{A- 



ML^^M^A) ,M\A) 



(3.13) 



where Eq. (|3.9| ) and its adjoint and Eq. (|3.11|) are used[5 Since the free part of the 
IMO commutes with itself and is independent of the cutoff, 
dMUA) 



ML,,MUA) ,M\A) 



(3.14) 



d{A-^) 

We assume that A4'^{A) can be written as a convergent (or at least asymptotic) 
expansion in powers of Ai?^^{A'): 



(3.15) 



a=0 



where V'-"-' is O {[J^int{A')]'^j , V^^^ = A1|ee; must determine the higher-order 

V(°)'s. Substituting Eq. ( p7[5|) into Eq. (|3A^ ), matching powers of A4f^^{A') on both 
sides, and doing some algebra yields 



d{A 



-4\ 



(3.16) 



6=1 



-•^^The restrictions that we place on the IMO in the next section prohibit AI^(A') from depending 
on A; so dM^{h')/dK = 0. 
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where the sum is zero if a = 0. This equation tells us how the O ([A^^nt(^')]'^) 
contribution to 7W^(A) depends on the cutoff in terms of lower-order contributions, 
and can be used to iteratively calculate A^^(A) in powers of M.f^^{A'). 

Taking a matrix element of Eq. between free states |/) and \F) for the case 

a = 1, we obtain 

'% = -V?i^h- (3.17) 



Solving this equation with the boundary condition 



M\A)\^^^,=M\A') (3.18) 



leads to 



VS = exp 



A'4 A4 



) A|,] MU^')fj. (3.19) 



A^fnt(^') is regulated such that Aif^^{A')Fi is proportional to exp(— A^.j-A'""'), which 
means that Eq. (|3.19| ) shows that Vpj is proportional to exp(— A|,jA~^). 

The higher-order V^'^'''s can be found the same way that V(^) was found: by taking 
matrix elements of Eq. ( |3.16|) and using the lower-order V^°'''s and the boundary 



condition, Eq. (|3.18|) . When products of A4f^^{A') are encountered, complete sets of 
free states have to be inserted between adjacent factors. The V*-"'''s that result from 
this procedure are all proportional to exp(— A|.jA~^), which shows that U (A, A') does 
indeed take an IMO with a cutoff A' and produce one with a cutoff A. 

If one computes V*^^^ and V*^^-*, then Eq. ( |3.15| ) gives A^^(A) in terms of J^f^^.{A') 
to O {[Ml,{A')f). Using this and 

{F\ -ML(A) |/) = {F\ V{A) \I) , (3.20) 
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which follows from Eq. (|3.2| ), one can show that the perturbative version of Eq. (|3] 
in terms of the reduced interaction is 

V{A)-V{A') = 6V, (3.21) 

where 5V is the change to the reduced interaction and is a function of both A and A': 

{F\6V\I) = WiF\ViA')\K){K\ViA')\I)Ti^'^'\F,KJ) 

^ K 

+ ^ E (^1 ^(A') 1^) (^1 V^(A') \L) (L| y (A') |/) (F, L, /) 

^ K,L 

+ O ([K(A')]') . (3.22) 

In this equation, the sums are over complete sets of free states and the cutoff functions 
are defined by 



and 



Tt''''\F,K,L,I) = 

.2A'-4A,.,,A„ / 2A-4A^^A,,, „2A'-4A^,,A 



1 1 \ / 1 1 \ 9A'-4a a / OA-^A A 2A'-4A a \ 

( e ^FK^KI — g^^^ ^FK^KI I 



e 



KL LI 



+ — — ^Fi^ ^ "^IK / 2A'-4(A^^A^,+A,,^A„) 

\^KL ^Lll ^KL^LI + ^FK^KI ^ 
_ ^2A-4{A^^A^,+A^^A„)\ , / _J: \_\ I _^ ^_ \ 2A'-^A^^A^^ 











1 




FK 





^ " ^ ' \^FK ^KlJ ^FK^KL + ^FL^l 

X (^e^^'~'(^FK^KL+^FL^Ll) - e^^-^^^FK^KL^^FL^Ll)\^ . (3.24) 

The above definitions for the cutoff functions assume that none of the A's that appear 
in the denominators is zero. In the event one of them is zero, the appropriate cutoff 
function is defined by 

t/^'^')(A = 0) = hm^r/'^'^')(A). (3.25) 
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3.3 Restrictions from Physical Principles 

Eq. ( p.21| ) is the first restriction on the IMO. To uniquely determine the IMO, we 
need to place additional restrictions on it, and we do this using the physical principles 
of the theory that are not violated by the cutoff. (See Footnote ^.) 

3.3.1 Symmetry Principles 

Any LFFT should exhibit manifest momentum conservation, boost covariance, 
and covariance under rotations about the 3-axis. Our cutoff does not violate any of 
these principles; so we restrict the IMO to conserve momentum and to be invariant 
under boosts and rotations about the 3-axis. 

3.3.2 Transverse Locality 

Ideally, the IMO should be local in the transverse directions, and thus each of its 
matrix elements should be expressible as a finite series of powers of transverse mo- 
menta with expansion coefficients that are functions of longitudinal momenta. In our 
case, the cutoff suppresses interactions that have large transverse-momentum transfers 
and replaces them with interactions that have smaller transverse- momentum trans- 
fers. This is equivalent to suppressing interactions that occur over small transverse 
separations and replacing them with interactions that occur over larger transverse 
separations; so we do not expect our interactions to be perfectly transverse-local. 
Nonetheless, we expect that interactions in Ai'^{A) should appear local relative to 
transverse separations larger than or, equivalently, to transverse momenta less 
than A. This means that for transverse momenta less than A we should be able to 
approximate each matrix element of Ai'^{A) as a finite power series in p±/A. We 
enforce this by assuming that transverse locality is violated in the weakest manner 
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possible, i.e. that any matrix element of the IMO can be expressed as an infinite 
series of powers of transverse momenta with an infinite radius of convergence. In 
other words, we assume that the matrix elements of the IMO are analytic functions 
of transverse momenta. 

3.3.3 Cluster Decomposition 

Cluster decomposition is a physical principle of LFFT that is partially violated 
by our cutoff; however, we can identify the source of the violation and still use the 
principle to restrict the IMO. To see how to do this, we consider a generic LFFT in 
d spacetime dimensions. Let us call the particles of this theory particles. Each of 
these particles has a momentum degree of freedom that we label with the index p 
and a set of discrete quantum numbers that we label with the index a. We denote 
the complete set of quantum numbers of particle i hy qf 

qi = {pi,ai}. (3.26) 

Momentum conservation implies that any matrix element of the reduced interac- 
tion can be written as a sum of terms, with each term containing a unique prod- 
uct of momentum-conserving delta functions. For example, in our generic LFFT, 
(^3^4! V{A) |^i^2) can be written in the form 

(661 V{A) 166) = [nM+Pt)3^''~'KPi +P2-P3-P4)] F^'\q,,q2,qs,q4,A) 
naptS^'-^^ipi - Pz)] [ndPt8^''~^^ip2 - Pa)] F^'\qi, q2, qs, q^, A) 
ndPtS^'~'\pi-P4)] [ndPtS'^''~'\P2-P3)] F'^'\qi,q2,q3,q,,A). (3.27) 



+ 
+ 



We have included a factor of Ud and a longitudinal momentum factor with each 
delta function because our normalization of states produces these factors naturally. 
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Cluster decomposition implies that the F^'^'s in Eq. ( p. 27] ) cannot contain delta 
functions of momenta It also implies that each term in Eq. ( p.27| ) has exactly one 
delta function associated with the conservation of momenta of interacting particles, 
with any additional delta functions associated with the conservation of momenta 
of spectators. This means that the first term on the RHS of Eq. ( pj.27D is the 
contribution from four interacting particles, and the second term on the RHS can have 
two contributions: one for which and ^3 are spectators, and one for which ^2 and 
^4 are spectators. Similarly, the third term on the RHS can have two contributions: 
one for which ^1 and ^4 are spectators, and one for which ,^2 and ,^3 are spectators. 

Normally, contributions to matrix elements depend on spectators only through 
quantum-number-conserving delta functions (Dirac delta functions for continuous 
quantum numbers and Kronecker deltas for discrete quantum numbers) and the cor- 
responding factors of longitudinal momenta, in which case we can write 

(^3^41 V^(A) 166) = \nM + pt)S^'~'\pi +P2-P3- P4)l F^'\q,, q,, q^, q^, A) 



+ 



nMS^'~'\P2-p^)] {5,,,3F(2'i)(g2,g4, A) 
ndptS^'-'^Pi-p^)] \ndptS^''-'\p2-P3) 



+ 5a„a,F^'''\quq3,A)} + 
X {6,^^,^F^^''\q,,q,,A) + 6a,,a,F'^'''\qi,q4,A)} 



(3.28) 



where we explicitly show the contributions from different sets of spectators and the 
Kronecker deltas that conserve their discrete quantum numbers. However, we use a 
cutoff on differences of free masses of states, and this cutoff violates cluster decom- 
position. To see this, note that the change in free mass for some process |/) — \F) 
is 



^FI 



(3.29) 
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where the p^'s are the momenta of the particles in |/), the p^s are the momenta of the 
particles in \F), and is the total longitudinal momentum of each state. The minus 
momenta of any spectators cancel in this difference, but their longitudinal momenta 
still contribute to the overall factor of P"^. Thus our cutoff partially violates the 
cluster decomposition principle, such that the F^^'-^^'s in Eq. ( |3.28D can depend on 
P^, the total longitudinal momentum of each state: 

{^3^41 V{A) 1^16) = \nM+Pt)S^'''\Pi+P2-P3-PA)] F^'\q,,q2,qs,q,,A) 



+ 



ndPi 



S^'-'\pi 



P3] 



ridPi 



5^'~'\P2-P,)\ {Sa,,a,F^'''\q2,q,,P^,A) 



+ 6a,,a,F^'''\quq3,P^,A)} + [naptS'^'''\Pi-P4)] [ndPtS'^'''\P2 - Ps) 
X {6a„a,F^''''\q2, qs, A) + Sa„a,F^''''\qi, q^, P+, A)} • 



(3.30) 



The result of this analysis can be used to formulate a general restriction on the 
IMO. When any of the IMO's matrix elements is written as an expansion in the pos- 
sible products of momentum- conserving delta functions, the coefficient of any term 
in the expansion is restricted as follows. (If there is more than one possible set of 
spectators for a given product of momentum-conserving delta functions, then the co- 
efficient, e.g. F^^-* in Eq. ( 3.27 ), has to be broken into a distinct part for each possible 
set, and these restrictions hold for each part separately.) It can depend on the cutoff, 
the quantum numbers of the interacting particles, and the total longitudinal momen- 
tum. It must be proportional to a quantum-number-conserving Kronecker delta for 
each discrete quantum number for each spectator. It can have no other dependence 
on the quantum numbers of spectators, and it cannot contain delta functions that fix 
momenta. 
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It is worth mentioning that if one uses a cutoff on free masses of states rather 
than differences of free masses, then the coefficients of the delta-function expansion 
can be functions of all the momenta in the matrix element, which is a signal of a 
much more severe violation of cluster decomposition. If one uses a cutoff on free 
energy differences, then the factor of in Eq. ( ^.29| ) is not present and cluster 
decomposition is not violated at all. However, longitudinal boost covariance is lost in 
this case because energy differences are not invariant under longitudinal boosts. It is 
possible to use cutoffs that maintain both cluster decomposition and boost covariance 



Tl| , but these complicate the process for computing the IMO and do not seem to 
provide any significant advantages. 

3.3.4 Representation of the Theory of Interest 

The preceding restrictions on A^^(A) are valid for any LFFT in more than two 
dimensions. In order to represent a particular theory, we must place additional re- 
strictions on A^^(A). 

We assume that we can compute the IMO perturbatively, which means that we 
can expand ^(A) in powers of the couplings at the scale A. Our cutoff has no effect 
in the noninteracting limit; so our IMO must reproduce free field theory in this limit. 
According to Eq. (|3.2|) , this means that V{A) vanishes in the noninteracting limit. 

In both theories that we are considering, the only fundamental parameter is a 
single coupling; so we require the IMO to depend only on it and the scale. (For an 
example of the application of our method to a theory with more than one parameter, 
see Ref. 0.) In this case, the expansion of ^(A) takes the form 



oo 

nA) = E^?:^^^nA), (3.31) 
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where g^^ is the couphng at the scale A. We refer to V^'^\A) as the O(y^) reduced 
interaction, although for convenience the coupling is factored out. 

Qj^ is the correct fundamental parameter for a theory if and only if its definition is 
consistent with the canonical definition of the coupling. The canonical definition in 
massless 0^ theory is 



9 



ndptS^'~'^iPl -P2- P3)l (02031 101 



(3.32) 



In this theory we choose 



9a 



napt5^''-^\pi-P2-P3)\~ (02031 M\K) |0 



IP2±=J'3±;P2 =P3 



naptS^''-'\pi-P2-P3)\~ (02031 V{k) I0i)|,^,,=,^3,;4=,+ - (3.33) 



In pure-glue QCD, the canonical definition of the coupling is 



9 = 



ridPt^^'^ ^\pi-P2-P3) ^92gM9iV\ {929z\M'^\gi) . (3.34) 



The denominator removes all dependence on momentum, spin, and color in the canon- 
ical matrix element for gluon emission, and thus isolates the coupling. For this theory 



we choose 



9a = 



ndPt^^"^ ^\PI-P2-Pz) ^929M9i~k\ {929z\M'^{K) \gi) 



Sn=l; Cn=n; e=0 



ndPt5^''-''\pi-P2-Pz)^929M9i-:^] (^2^3! ■^^(A) [ ^ ,(3.35) 

P2X=P3X;Pj=P3 



where n— 1, 2, 3. 

Note that in each theory, momentum conservation and boost invariance imply 
that the matrix elements that define the coupling can depend only on the transverse 
momentum of particle 2 in the center-of-mass frame and the ratio pt /Pi- The re- 
strictions p2i_ — P3± and P2 — pt fix these quantities to be and 1/2, respectively. 
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In each case, our definition of tlie coupling is consistent witli tlie canonical definition 
because the conditions on the matrix elements in Eqs. ( p.33| ) and ( p.35| ) have no 



effect on the right-hand sides of Eqs. ( p.32| ) and ( p.34| ), and these conditions do not 



force (02031 |0i) or {g2g3\ M"^ \gi) to vanish. According to Eq. ( p.31| ), M'^{A) is 
coupling coherent WR, K7\, EH] in each theory because the couplings of its noncanonical 



operators are functions only of the fundamental parameters of the theory and they 
vanish in the noninteracting limit. 

We have assumed that the IMO obeys approximate transverse locality, which 
means that we can expand any matrix element {F\ V"*^^'-'(A) |/) in powers of trans- 
verse momenta. Each term in this expansion is either cutoff-dependent or cutoff- 
independent. We define Vq^I^) be the cutoff-dependent part of V^^'"-*(A), i.e. the 
part that produces the cutoff-dependent terms in transverse-momentum expansions 
of matrix elements of V^^'\A). We define Vqi to be the cutoff-independent part 
of V^'^\A), i.e. the part that produces the cutoff-independent terms in transverse- 
momentum expansions of matrix elements of V^^\A). Then 

V^^\A) = V^i{A) + vi\\ (3.36) 

This separation is necessary because the procedures for computing Vq^{A) and V^j'' 
differ. 

If A4'^{A) is to reproduce the second-order scattering amplitudes of massless 0^ 
theory, then it is necessary and sufficient for the reduced interaction to satisfy the fol- 
lowing conditions[^. First, we prohibit V^'^\A) from having a three-point interaction 

unless r is odd. Second, we require 

'^'^A proof of this statement exists, but the proof is tedious and is awkward to fit into our presen- 
tation; so we omit it. 
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+ 44M^'"'Hp2+P3-Pi)]- (3.37) 

On the other hand, we beheve that if Al^(A) is to reproduce the perturbative 
scattering ampUtudes of pure-glue QCD, then it is necessary and sufficient for the 
reduced interaction to satisfy the following conditions'^. We prohibit V'^'^\K) from 
having a three-point interaction unless r is odd, and we prohibit V^'^\h) from having 
a four-point interaction unless r is even. We require 

+ / D^D^D^a\a^a25^''-^\p^+p2-pz) ^9M9i92-^ (3-38) 

and 



y(2) _ 



4 

jDiD2DsD^alalaia2S^''-'\pi + P2 - Ps - Pa) J2^[%s,4 ^9s9M9i92:^i 
■ ■ 1=1 

4 

+ V^"^ (d^D^D^D^ 444ai 5^"-^^ {p, -P2-P3- Pd E ^Jl3,4 ^52^354kbi:f i 

4 

+ 7^+^ [D,D2DsD^alaia2as6^''-'\p, + p2 + P3 - P4) E ^ilw 4:54^^152^3:^^, 

i=l 

(3.39) 

where 

'^l,2;3,4 — -L, 

^''Any longitudinal regulator that is consistent with the physical principles that we use to restrict 
the IMO is sufficient. It is not necessary to use our e cutoff in order to reproduce the perturbative 
scattering amplitudes of pure-glue QCD. 
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and 



and 



fl(2) 
^l,2;3,4 

^(3) 

^^3,4 = 1, (3.40) 



^(1) 

1,2,3;4 
1,2, 3;4 



,2,3;4 



0{\pt 


-pt\ 




e{\pt 


-pi\ 




1, 






J-, 






d{\pt 


-pt\ 




e{\pt 


-pt\ 




e{\pt 


-pt\ 




1 

J-, 








-pt\ 






~pt\ 






-pt\ 





/](!) 
^1;2,3,4 

n(2) 
Cl;2,3,4 

/l(3) 
1;2,3,4 

(4) 



^i:^,3,4 = Oi\pt-pt\-eV'-), (3.41) 



Ci3;4 = ^dP^-pJl-eP"')- (3.42) 

The presence of e in these step-function cutoffs ensures that we will avoid divergences 
from exchanged gluons (either instantaneous or real) with infinitesimal longitudinal 
momentum. In Chapter § we show how we can take e to zero before we diagonalize 
M\A). 

We have not yet proved that the above conditions on the reduced interaction 
are necessary and sufficient to reproduce the perturbative scattering amplitudes of 
pure- glue QCD. 
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3.4 The Recursion Relations for the Invariant- Mass Opera- 
tor 

In each of the theories that we are considering, the restrictions that we have placed 
on the IMO are sufficient to allow us to derive recursion relations that determine the 
IMO to all orders in perturbation theory. In this section, we derive these recursion 
relations. We want to again point out that the restrictions that we have placed 
on the IMO are based on a subset of the physical principles of the theory. The 
remaining principles, such as Lorentz covariance and gauge covariance (in a gauge 
theory), must be automatically respected by physical quantities derived from our 
IMO, at least perturbatively. If they are not, then they contradict the principles we 
use and no consistent theory can be built upon the complete set of principles. The 
reason that this process is possible is that there are redundancies among the various 
physical principles. In pure-glue QCD for example, gauge covariance is related to 
our requirements that the IMO produce the perturbative scattering amplitudes of 
pure-glue QCD and be a function only of and A. 

3.4.1 Removal of the Coupling 

To begin, we consider the restriction that forces the IMO to produce cutoff- 
independent physical quantities: 

V(A)-y(A') = 5V. (3.43) 

This restriction is in terms of the reduced interaction and the change to the reduced 
interaction. To derive the recursion relations, we need to remove the coupling from 
this equation. To do this, we note that since V^(A') can be expanded in powers of gf^, , 
so can 5V: 
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(3.44) 



t=2 



We refer to 5V^^^ as the 0{gj^^) change to the reduced interaction, although for con- 
venience the couphng is factored out. Note that 5V^^^ is a function of A and A'. 
Now Eq. (|3.43| ) can be expanded in powers of and g^, : 



(3.45) 



t=l 



t=l 



t=2 



This equation is a bit tricky to use because it involves the couplings at two different 
scales. To see how they are related, consider massless (p^ theory. The matrix element 
of Eq. (lO^ ) for 0i ^ 0203 is 



2031 ViA) |0i) - (02031 ViA') |0i) = (02031 6V |0i) 



According to the definition of the coupling, this equation implies that 



(3.46) 



9a - 9a, 



P2 - P3) 



(02031 5^ |0l) 



P2±=P3±; Pj=P3 



(3.47) 



Since V^^'^ changes particle number by 1, inspection of Eq. (|3.22|) reveals that 
(02031 5V-|0i) is 0{giy,so 



9A=9Ai+0{gl). 



(3.48) 



The same analysis can be done in pure-glue QCD, and it also leads to Eq. ( 3.48| ). 
Thus in each theory. 



9a=9a.+119:,Cs{A,A'1 



(3.49) 



s=3 



where the C^'s are functions of A and A'. For an integer t > 1, Eq. ( |3.49|) implies 
that 
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9i=i,+Y.9TBU^^'), (3.50) 

s=2 



where the Btjs are also functions of A and A', and can be calculated in terms of the 
Cs's by raising Eq. ( p.49| ) to the t*'^ power. 



We substitute Eq. ( ^.5(J| ) into Eq. ( ^.45| ) and demand that it hold order-by-order 
in g^,. At C ((?[,) {r > 1), this implies that 



r-1 



V-W(A) - \/M(A') = - J2 Br-s,sV^""\A), (3.51) 

s=2 

where ^V^'-^^ = 0, and we define any sum to be zero if its upper limit is less than 
its lower limit. The cutoff-independent parts of V^^\A) and V^*^'")(A') cancel on the 
left-hand side (LHS), leaving 

F(52(A) - V^iiA') = <5yM - J2Br-s,sV^'-'\A). (3.52) 

s=2 

This equation can be used to derive the desired recursion relations. 

3.4.2 The Recursion Relation for the Cutoff-Dependent Part of the Re- 
duced Interaction 

We first compute the recursion relation for the cutoff-dependent part of the re- 
duced interaction. In each theory that we consider, we already know the first-order 
reduced interaction [see Eqs. ( ^.371 ) and ( ^.381 )]; so we need to compute V^^ fo^' ^ > 2, 
in terms of the lower-order reduced interactions. 

Recall that momentum conservation implies that any matrix element of V{A) 
can be written as an expansion in unique products of momentum-conserving delta 
functions. This means that an arbitrary matrix element of Eq. (|3.52| ) can be expanded 
in products of delta functions: 
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E 



r-1 



{F\ 5V^M |/) - Y: Br-s,s {F\ V'^^-'\A) I J) 



s=2 



(i) 



(3.53) 



where the (i) superscripts denote that we are considering the i*^ product of delta func- 
tions that can occur in a delta-function expansion of (F| V^'^\A) This equation 
is equivalent to a set of equations, one for each possible product of delta functions: 



{F\V^,i{A)\lf^-{F\viti^')\lf^ 

lii) 



{F\ SV^'^ \I) 



r-1 



Ei?r-.,.(F|y(-)(A)|/) 



s=2 



(3.54) 



Cluster decomposition implies that we can write 



{F\viiW\I)^'> = 



where Sj'^ is the j*'^ delta function in the i^^ product of delta functions (it includes 
the longitudinal- momentum factor), Ng^ is the number of delta functions in the i*^ 
product, and -^ci)({Pn}; {c^nji A) is a function of the cutoff and the quantum numbers 
of the particles in the matrix element, but does not contain delta functions that fix 
momenta. We define A^part to be the number of particles in state |/) plus the number of 
particles in state |F), and n = 1,2,3,..., Apart- Here p„ is the momentum of particle 
n and a„ is its set of discrete quantum numbers: a„ = c„} in pure-glue QCD, and 
a„ is the null set in massless (f)^ theory. We define N-^^^ to be the number of particles 
in the matrix element that participate in an interaction for the i^^ product of delta 
functions. In order for the IMO to have the dimensions (mass)^, must have the 
dimensions (mass)**^*', where b^^^ — 6 — 2N^^^^ in massless 0^ theory and b^^^ — A — N-^^^ 
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th 



^^?)({Pn},K},A), 



(3.55) 



in pure-glue QCD. Note that we are suppressing the dependence of the RHS of Eq. 
( p.55| ) on r. 

We have assumed that any matrix element of the IMO can be expanded in powers 
of transverse momenta, not including the momentum-conserving delta functions; so 



{■mnt} 



■({p^},K}) n n (¥) ' (3-56) 



n=l t=l 



A 



where t denotes a component of transverse momentum and mnt is a non-negative 
integer index associated with transverse-momentum component t of particle n. The 
sum is over all values of each of the rant's, subject to the constraint that 



(3.57) 



n,t 



which is necessary to avoid terms in the expansion that are cutoff-independent. The 
zt"'^'s are the coefficients for the momentum expansion. They depend on the m„t's 
and are functions of the particles' longitudinal momenta and discrete quantum num- 
bers. 

Since the RHS of Eq. ( |3.54| ) has the same product of delta functions as the LHS, 
we can write 



r-l 



s=2 



N 



(i) 



(0 



X G«({p„},{a„},A,A'), 



(3.58) 



where G^*^ has dimensions (mass)^*'', and by inspection of the LHS of Eq. ( ^.541 ) 
and Eq. ( 3.55|) is a function of the quantum numbers of the particles, A, and A'. 
Substitution of Eqs. (|3.58|) and ( p.55|) into Eq. (p.54|) yields 



Fg{{pn},{an},A)-F^I,i{pn},{an},A') = G^'\{pn},{an},A,A'), (3.59) 
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where the momenta in this equation are constrained by the delta-function conditions. 

Since the LHS of Eq. ( |3.59| ) is the difference of a function of A and the same 
function with A A', G*^*^ must be as welL Since the LHS of Eq. ( |3.59| ) can be 
expanded in powers of transverse momenta, G*^*-* must have the form 



G«({p„},{a4,A,A') = Y: zt"'\{pt},Wn} 



{m„t} 



(3.60) 



where the sum is restricted by Eq. ( |3.57|) . 

Substituting Eqs. ( |3.6CI| ) and ( p.56| ) into Eq. ( |3.59| ), we find that 



{mnt} 



{■m„t} 



{{pth{an})Y[r-f 

n,t \ 

Ct \ ™nt 



{rant} 



A' 



.{3.61) 



Matching powers of transverse momenta on both sides of this equation gives 



(3.62) 



The factor in brackets cannot be zero because A 7^ A' and Eq. (|3.57|) holds. Thus 
Eq. (|3.62|) implies that 



{m„t} 



z 



{rrint} 



(3.63) 



Then Eqs. (|336|) , (ICTp , and (|333|) imply that 



Fl^UiPn}, {a J, A) = G« ({p4, {a4, A, A') I 



A terms 



(3.64) 
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where "A terms" means that G'-*-' is to be expanded in powers of transverse momenta 
and only the terms in the expansion that are proportional to powers or inverse powers 
of A contribute. From Eqs. ( ^35]) , ( p38D , and (^), 



(F|r^^:)(A)|/)« 



r-l 



(3.65) 



•5—2 J A terms 

where it is understood that the momentum-conserving delta functions are ignored for 
the purposes of transverse-momentum expansions. Since a matrix element is the sum 
of the contributions to it from different products of delta functions, both sides of this 
equation can be summed over i to obtain 



(3.66) 

s=2 J A terms 



This equation tells us how to calculate the cutoff-dependent part of the 0{g^^ reduced 
interaction in terms of lower-order contributions. 

3.4.3 The Recursion Relation for the Cutoff-Independent Part of the Re- 
duced Interaction 

To complete the solution for the IMO, we need to specify how to compute the 
cutoff-independent part of the reduced interaction. In massless 0^ theory, we have 
specified y*^^-*; so we need to determine Vqi for r > 2. In pure-glue QCD, we have 
specified V"*^^-* and V^i''; so we need to determine V^i^ for r > 3. Before we derive 
the recursion relation, we first consider which contributions to V^"^"^ (A) can be cutoff- 
independent. 

A matrix element of the cutoff- independent part of V"*^''^(A) can be expanded in 
products of delta functions and in powers of transverse momenta just as was done for 
the cutoff-dependent part. Thus we can write 

{F\V^'hl)=T.^F\yS\lf\ (3.67) 

i 
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where 

I T\ (*) 



{F\V^,'\I)' 

and 



j 



F'Si{Pn},{an}) (3.68) 



FmiPn}, {an}) = A''' E ^t"'\{pt}^ M) n n ¥ , (3-69) 

{m„t} ri=l t=l \ / 

where the sum is over all values of each rrint, subject to the constraint 

6»-E"^ni = 0. (3.70) 

n,t 

Eq. (|3.70| ) ensures that all the terms in the expansion of F^l are cutoff-independent. 

Eq. (13.701 ) places constraints on the possible cutoff-independent contributions to 
the reduced interaction. Any contribution to a matrix element of ^^'^''"^(A) has an 
> 2, but Eq. (ITfOl) can only hold if N^^^ < 3 in massless (j)^ theory and if 



Ai2 < 4 in pure-glue QCD. 

Suppose that N^^l = 2. In this case, Eq. ( p.70|) implies that Fq} is quadratic 
in transverse momenta. Due to approximate cluster decomposition, only interacting 
particles' transverse momenta can appear in Fq^. So any contribution to F^l can 
depend on the transverse momenta of two interacting particles. Thus F^j can be 
written as a sum of terms, where each term corresponds to a distinct pair of interacting 
particles. The momentum dependence of each term in FqI is limited to dependence 
on the momenta of the interacting particles and the total longitudinal momentum: 

m 

= E4rHA;„,P+,{a„}), (3.71) 

m 

where and k'^ are the momenta for the initial and final particles in the m^^ 
interacting pair, and where we have used the fact that for N^^^ = 2, momentum 
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conservation implies that km = k'^- Fq(^^ must be quadratic in kmi. or it must be 
zero. The matrix elements of the IMO are boost-invariant, as is the delta-function 
product in Eq. ( |3.68| ). This means that FqI must be boost-invariant, but it cannot 
be if Fqj"*^ is quadratic in so Fq^^'^ must be zero. Thus the reduced interaction 
does not contain any cutoff-independent two-point interactions. 

Note that two-point interactions are self-energies, and they change the particle 
dispersion relation. If they change the dispersion relation such that the coefficients 
of the free relation become modified by interactions, then this can be viewed as 
renormalization of the field operators, i.e. wave-function renormalization. This effect 
is absent unless either Fq\ or FqI^ can be quadratic in transverse momenta for A^/*'' — 



int 



2. We have just shown that boost invariance prevents this for Fq\^ and according to 
Eq. ( |3.57| ), Fq\^ cannot be quadratic in transverse momenta for N^^l = 2; so there is 
no wave-function renormalization at any order in g^^ in our approach. 

At this point the analyses of the cutoff-independent parts of the reduced inter- 
actions in the two theories diverge. We first consider massless (f)^ theory and then 
pure-glue QCD. 

According to Eq. ( 3.701) , in massless (fy^ theory, the only other possible contri- 
bution to Vqy is from transverse-momentum-independent three-point interactions. 
According to the discussion in Subsection p.3.4| , V^'^\K) cannot have a three-point 
interaction unless r is odd. This means that if r is odd, then V^j^ can contain only 
three-point interactions that are independent of all transverse momenta, and if r is 
even, then Vqi = 0. 
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Since V^^ = if r is even, for the remainder of this massless </)^ discussion, we 
assume that r is odd. Then we need to consider only N^^^^ = 3. To calculate Vqi , we 



consider Eq. ( p.54| ) with 



r + 2: 



{F\ |/) - X: 5.+2-.,s {F\ ^(^•+2-^)(A) |/)' 



s=2 



(3.72) 



We expand Eq. ( |3.72| ) in powers of transverse momenta and keep only the constant 
term: 







r+1 



s=2 



(0 



p'j' term 



(3.73) 



where we have used the fact that {F\ Vq^^^(A) |/)'~*'' has no constant part when N-^ 



r(0 

int 



3 [see Eq. ( p.57| )]. We move the first term in the sum on the RHS to the LHS: 



Br,2 {F\ V^\> \I) 



(r) I T\{i) 



r+1 



r+2-s) 



(A) |/) 



s=3 



(i) 



(3.74) 



Now we can sum over all values of i corresponding to three-point interactions. Then 
for any even r, Vqi = 0, and for any odd r > 3, 



{F\ VS \I) 



r,2 



r+1 



s=3 



3— point 



. (3.75) 



term 



This equation tells us how to compute the cutoff-independent part of the 0{g2) re- 
duced interaction for odd r > 3 in terms of the cutoff-dependent part and lower-order 
reduced interactions. To use the equation, the RHS has to be expanded in powers 
of transverse momenta and only the constant (transverse-momentum-independent) 
three-point interactions contribute. This is an integral equation^ because V^j"* ap- 
pears on the RHS of the equation inside integrals in SV^'^~^'^\ A cursory examination 

^^It is very difficult to prove that integral equations of this type have a unique solution; so we 
simply assume it is true in this case. 
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of the definition of 5V^^~^'^^ may lead one to believe that Eq. (|3.75|) is useless because 
it requires one to know V^^^^'^{A). However, since r is odd, V^'''^^\A) = Vqy)^\A), 
and the dependence of 6V^^~^'^^ on V^j^^^ (A) can be replaced with further dependence 
on V^j^ and Vq^{A) (and lower-order reduced interactions) using Eq. ( ^.66| ) with 
r — > r + 1. 

To conclude our discussion of massless 0^ theory, we would like to deduce a bit 
more about the relationship of to g^, . According to Eq. ( p.47|) and the surrounding 
discussion, this relationship is determined by the matrix element (0203| which 
can be expanded in powers of g^, : 

oo 

(02031 5V |0i) = J29l, (02031 5V^'^ |0i) . (3.76) 

t=3 

Recall that (5V^*^*^ is built from products of V^*^''^(A')'s. This implies that SV^^"^ can 
change particle number by 1 only if t is odd, and thus Eq. (|3.76|) implies that the 
coupling runs at odd orders; i.e. Cg is zero if s is even [see Eq. ( p.49|) ]. 

We now turn to the case of pure-glue QCD. In this theory, it is possible for there 
to be both three-point and four-point interactions in Vqi . According to Eq. ( 3.7CI| ), 



if N^^^ = 3, then FqI has to be linear in transverse momenta, and if N^^^ = 4, then 
FqI has to be independent of all transverse momenta. According to assumptions that 
we made in Subsection 3.3.4 , if r is odd then V^j'' has no N^^^ = 4 part, and if r 
is even then Vqi has no N^^^ = 3 part. This means that if r is odd, then V^i'' can 
contain only three-point interactions that are linear in transverse momenta, and if r 

(r) 

is even, then V^j can contain only four-point interactions that are independent of all 
transverse momenta. 

To calculate V^j'', we consider Eq. ( |3.54|) with r ^ r + 2: 

(^1 v![^'\a) |J)« - {F\ V^''^'\a') |/)« = 
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r+1 



{F\ |/) - Y: Br+2-s.s {F\ V^'-+''^\A) \I) 



s=2 



(3.77) 



For the remainder of this discussion, we assume that r is odd. Then we need to 
consider only A^j^^ = 3 initially. We expand Eq. ( |3.77| ) in powers of transverse 
momenta and keep only the linear term: 







r+l 



{F\ |/) - ^ Br+2-s,s {F\ V^^+'-^\A) \I) 



s=2 



p I term 



(3.78) 



where we have used the fact that {F\ Vqy)'^\A) has no linear part when A^- 



[see Eq. ( |3.57| )]. We move the first term in the sum on the RHS to the LHS: 



Br,2 (FIV^C^^ |/)^^^ 



-r+l 



s=3 



. (3.79) 



-I p I term 



Now we can sum over all values of i corresponding to three-point interactions: 



{F\ V^f,^ \I) 



Br' 



r+l 



s=3 



3— point 



(3.80) 



term 

.(r+l) 



To use this equation, we also need to use Eq. ( p.77| ) to solve for V^j • Since r is 
odd, r + 1 is even. Thus V^i^^'' will contain only transverse-momentum-independent 
four-point interactions. Using steps analogous to those that led to Eq. ( p.80|) , we find 
that 



{F\ Vt'^ \I) 



B. 



r+1,2 



r+2 



s=3 



4— point 



(3.81) 



term 



To use these equations, the right-hand sides have to be expanded in powers of trans- 
verse momenta. Only three-point interactions that are linear in transverse momenta 
contribute to Vqi , and only four-point interactions that are independent of all trans- 
verse momenta contribute to V^i''^^^ 

These equations are coupled integral equations because both V^j'' and V^j^^^ ap- 
pear on the RHS of Eq. ( CT) inside integrals in 5V"(^+2), and V^j^^'^ appears on the 
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RHS of Eq. (|;8T|) inside integrals in SV^'^+^l It would seem that V^'j^^^ also appears 
on the RHS of Eq. ( p.Slj ) inside integrals in 6V^^~^^\ but Vqi^^'' cannot couple to V'^^^ 
to produce a transverse- momentum-independent four-point contribution to 5V^'^~^^\ 
This is because the cutoff function Tg^'^ ^ vanishes when the intermediate state has 
zero free mass and all the external transverse momenta are zero. This means that 
since we specified V^^'' and in Subsection |3.3.4| , we can use Eqs. ( |3.8U| ) and 
( ^.81| ) to solve for V^j and V^i^ simultaneously, and V^j'* and Vqj"* simultaneously, 
and so on. Note that before we can use these equations to solve for V^j'' and V^i'''^'' 
simultaneously, we must first use Eq. (|3.66| ) both to compute Vqj^{A) in terms of 
lower-order interactions and to express Vq^)^'' (A) in terms of lower-order interactions 
and 1/W(A). 

To conclude our discussion of pure-glue QCD, we note that just as in the case of 
massless 0^ theory, ^V^*^*) can change particle number by 1 only if t is odd, and thus 
the coupling runs at odd orders; i.e. Cg is zero if s is even. 
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CHAPTER 4 



CALCULATION AND ANALYSIS OF MATRIX 
ELEMENTS IN MASSLESS 0^ THEORY 

In this chapter we demonstrate how to use the recursion relations for the IMO by 
calculating and analyzing some second- and third-order matrix elements of the IMO 
in massless 0^ theory. We work in six dimensions explicitly. 

We first calculate some second-order matrix elements and show that they obey 
the restriction of approximate transverse locality. We then show how our IMO yields 
cutoff-independent perturbative scattering amplitudes, even though it explicitly de- 
pends on the cutoff. We proceed to calculate the emission matrix element of the IMO 
to third order (neglecting the cutoff- independent part of the reduced interaction), in 
terms of five- dimensional integrals that we compute with Monte Carlo methods. In 
the process of deriving this matrix element, we derive the cutoff dependence of the 
coupling and show how it relates to the standard result that one finds with dimen- 
sional regularization. 

4.1 Second- Order Matrix Elements 

The matrix elements of the IMO are defined in terms of the matrix elements of the 
reduced interaction by Eq. ( p.2|) . To calculate matrix elements of the IMO to second 
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order, we need to calculate the corresponding matrix elements of V^^^^(A), the second- 
order reduced interaction. As mentioned in the previous chapter, for V"*^'"''(A) to have 
a cutoff-independent part, r must be odd. According to Eq. ( p.66| ), this means that 



V"*^^-'(A) is given in terms of the second-order change to the reduced interaction by 

{F\V^'\A)\I) = (i^|5V^^'M/)L,_ ' (4-1) 

where the second-order change to the reduced interaction is defined in terms of the 
first-order reduced interaction and the cutoff function: 

{F\6V^^'>\I) = -Y,{F\V^^^\K){K\V^^^\I)Ti^'^'\F,K,I). (4.2) 
2 K 

4.1.1 Example: (02|A^^(A) |0i) 

As a first example, we compute the matrix element (02| A^^(A) \(f)i) to 0{g^). To 
compute this matrix element, we use |J) = and \F) = \4>2)- Then the matrix 
element of the second-order change to the reduced interaction is 

{^2\5V^'^\<Pi) = W{h\V^'^\K){K\V^'^\<P^)Tt'^'\F,K,I). (4.3) 

^ K 

Since V^^^ changes particle number by 1, \K) has to be a two-particle state and Eq. 
( [4.3| ) becomes 

(02| SV^'^ = y D,D, V(i) 10304) (03041 V^^'^ |0i) T^^'^i^, I), (4.4) 

where the completeness relation in Eq. ( |2.57| ) is used and li^) = |0304). We represent 
{(f)2\5V^'^^ |0i) diagrammatically in Figure 4.1. Using the expression for V^^^ in Eq. 



X Tf'^')(F,i^,J). (4.5) 
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Figure 4.1: A diagrammatic representation of (02 1 SV^"^"^ |</)i), a matrix element of the 
second-order change to the reduced interaction. The numbers label the particles. 



The integral over can be done with one of the delta functions. At this point it is 
convenient to change variables from p^ and p4± to x and r±: 

Pa = {P4,Pa±) = {xpt,xpii_ + fi_), (4.6) 

where x is the fraction of the total longitudinal momentum that is carried by particle 
4 and r± is the transverse momentum of particle 4 in the center-of-mass frame. We 
only display the longitudinal and transverse components of the momentum. (Since 
the states that we use are free states, and since the momentum of a free massless scalar 
particle satisfies p^ = 0, the minus component of the momentum of 0j is constrained 
to be given by p^^ = pf^/pt ■) The momentum- conserving delta functions imply that 



P3 



([1 - x\pt, [1 - x]pi^ - f^) . (4.7) 



Note that all longitudinal momenta are positive, although in this chapter we do not 
display the step functions that enforce this. This implies, for example, that < x < 1. 



With the change of variables, Eq. (|4.5| ) becomes 



= \pt5^''\p2 - Pi) I d'r^dx . ^ ^ Tf -^'^F, K, I). (4.8) 
4 J x{\ — X) 

We are using massless particles; so the free masses of the states |/) = and 
1-^) = 102) are zero: 
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Mf = M^ = 0. (4.9) 

The particles in the intermediate state li^) = |0304) are also massless, but the state 
still has a free mass from the particles' relative motion: 

Mi = (4.10) 

X{1 — X) 

Then from the definition of the cutoff function Tg^'^ \ 

Tt'''\F^ I) = -^-^^^^ [e-^^'-^^ - e-^^-^^) . (4.11) 



After the integral is done, Eq. (|4.8| ) becomes 

W^I^V'^'W = ''MYJ^yfCA'^-A"). (4.12) 
To apply Eq. (4J.), we recall that we are to expand the RHS of Eq. (|1|) m powers 



of transverse momenta and keep only those terms that are proportional to powers 
or inverse powers of A. This yields the matrix element of the second-order reduced 
interaction: 

(0,|V(^)(A)|0,) = 5i,2A^Y5^/f. (4.13) 

Using the definition of the IMO in terms of the reduced interaction, we find that to 
second order, the matrix element of the IMO is given by 

{<^,\M\A) = 5,,A^ ^fl^/f- (4-14) 

This matrix element of the IMO has been completely determined to second or- 
der, in terms of the fundamental parameters of the theory, by the restrictions that 
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the IMO has to produce cutoff-independent physical quantities and has to respect 
the unviolated physical principles of the theory. If we wanted to compute physical 
quantities, this is one of the matrix elements of the IMO that we might need. We 
would have to choose A and determine Qj^ by fitting data since the theory contains 
one adjustable parameter. Note that Eq. ( ^.141 ) is consistent with all the restrictions 
that we have placed on the IMO. 

4.1.2 Example: (03041 A<^(A) 10102) 

Next we calculate the matrix element (0304| A^^(A) |0i02) to 0{g'^. We use 
|/) = 10102) and \F) = |0304). Then the matrix element of the second-order change 
to the reduced interaction is 

(03041 10102) = ('^3041 V^''^ \K) {K\ |0i02)Tf J). (4.15) 

Since V^^^-* changes particle number by 1, \K) has to be either a one- or three-particle 
state. When V^^^ is substituted into the RHS of this equation and all the creation 
and annihilation operators from the possible intermediate states and the interac- 
tions are contracted, we see that (0304 1 5 V^^-* |0i02) has contributions from a num- 
ber of different processes. We define (0304 1 ^V^*-^-* |0i02)s to be the contribution to 
(03041 5V^^'^ 10102) from the process shown in Figure 4.2s, where s — cL, b, c, . . . , i. 
Then 

(03041 SV^^^ 10102) = E (<^304| 5^(2) 10102), . (4.16) 

Eq. ( [4. ID suggests that we define a contribution to the matrix element of the 
second-order reduced interaction for each contribution to (0304| SV^"^^ |0i02): 



(03041 V^(') (A) 10102), = {M5V^'^\M2)s\a,,,^,, (4.17) 
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(a) 







Figure 4.2: A diagrammatic representation of the contributions to (0304| dV^"^^ \<Pi<p2), 
a matrix element of the second-order change to the reduced interaction. The numbers 
label the particles. 
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and then 

(03041 V^^'^A) 10102) = E('/'304|V^^'nA) 10102),. (4.18) 

For each contribution to (0304| 5V^'^^ |0i02), we use the change of variables 
P2 = ([l-x]P+,[l-x]P^-r-l), 

P3 = iyP^,yP± + q±), 

P4 = ([l-y]P+,[l-2/]Px-g±), (4.19) 
where P = pi + p2. Then the initial and final free masses are 



x{l — x) 



Ml = (4.20) 



The Annihilation Contribution 



The contribution to (0304| 5V^'^'> |0i02) from the annihilation process, shown in 
Figure 4.2a, is 

(03041 5V^('M0102)a = \ j (03041 V^^'M05)(05|V^^'M0102)rf''''^(F,K,/) 

= 32.W)(p,+p.-,3-P4)(^ + ^) 

^ (^g-2A'-A4M| _ g-2A-A4Mf) ^ (4 21) 



where the intermediate state is \K) = |05). According to Eq. (|4.17|) , the correspond- 



ing part of the matrix element of the second-order reduced interaction is 

(03041 V^'\A) 10102), = 327t'P+6^'\pi+P2-P3~P4) (^-^ + ^ 



^ I g-2A'-4M|, Af2 _ g-2A-4M|,Af2 
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(4.22) 

A terms 



If we expand everything multiplying the delta function in powers of transverse mo- 
menta, the lowest-order terms from the two exponentials cancel and leave two types of 
terms: those that are proportional to inverse powers of A and those that are propor- 
tional to inverse powers of A'. We isolate the terms that are proportional to inverse 
powers of A without altering the cancellation of the lowest term by replacing the first 
exponential with a 1: 

ihM V^'\A) \M2), = 327r^P+5(^)(pi +P2 -P3 -P4) + ^) 

X (i _ e-2A-'^l^f ) . (4.23) 

The Two-Particle Self-Energy Contributions 

The rest of the contributions to (0304| dV^"^^ |0i02) have a three-body intermediate 
state: \K) — 1 050607)- The contribution from the two-particle self-energy, shown in 
Figure 4.2b, is given by 

(03041 51^^') |0i02)b = \— [ D,D,DrTi''''''\F,K,I)6A7r'ptd^'>\p3-p5-p,) 

X 647rV<^^'HPi - P6 - Ps) <^7,4 <^7,2 ■ (4.24) 

We use the change of variables 

Pb = (^Pi",^Pi± + ^±), (4.25) 

and do the integrals over pe and pj using the delta functions. Then we find that 

P6 = ([l-z]p+,[l-2:]m-^±)- (4.26) 

The free mass of the intermediate state is 

Ml = + , ^ , , (4.27) 

^ x{l-x) z{l-z)x' ^ ' 
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and 



— _2A'-4a2 



According to Eq. (|4.17|) , the corresponding part of the matrix element of the second- 
order reduced interaction is 

(03041 ^^'nA)|0102)b = ^1-3 52,4 A' (4.29) 



Note that Eqs. (|4.13| ) and (|4.29|) indicate that to 0{g^) in this theory, the only 



effect a spectator has on the self-energy is to produce an extra delta function. If our 
cutoff respected cluster decomposition, this would be all that could happen. However, 
our cutoff partially violates cluster decomposition. This means that in principle the 
self-energy could depend on x, the fraction of the total longitudinal momentum that 
is carried by particle 1. This requires a violation of cluster decomposition because 
X = pt / {pt ~^pt)y P2 is the momentum of a spectator (see Figure 4.2b). Although 
this manifestation of the violation of cluster decomposition is not seen at second order 
in this theory, it is in pure-glue QCD [see Eq. ( p. 741 )]. 
According to Figure 4.2, 

(03041 V^^\K) 10102), = (03041 V^^\K) |0102)b|3^4 = '^1,4 

(03041 ^(^^A) |0102)d = (03041 V(')(A) 10102)^13^4 = '^1.3 '^^,4 A^J^/f, 



(03041 V^^\^) 10102), = (03041 V^^\K) |0102)d|3„, = ^1.4 ^^,3 A^ ^^'^^^ 
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The Exchange Contributions 

The exchange contribution to the matrix element of the second-order change to 
the reduced interaction, shown in Figure 4.2f, is 

(03041 |0i02)f = ^^(Y^ / D,D,DjTt'^'\F,K,I) 
X 647r^p+5(^) (p4 - P6 - P7)647r^p^5(^) (pi - - pe) ^7,2 ^5,3 



-327r5p+5(^)(pi+P2-P3-P4)- 



1 



x-y 



1 1 

+ 



X 



— 2A'~4a a — 2A^4a a 



(4.31) 



where the free mass of the intermediate state is 



Ml 



1 — X y 



x-y 



(4.32) 



Eq. ( [4.17] ) implies that the corresponding part of the matrix element of the second- 
order reduced interaction is given by 



(03041 V^^'^A) |0102)f = -327r^P+5(5)(j9i+p2-p3-p4) 



1 



X-y 



1 1 

+ 



X 



1 _ e-^^'^^KF^KI 



(4.33) 



According to Figure 4.2, the remaining parts of the matrix element of the second-order 
reduced interaction are given by 

(03041 V^(')(A) |0102)g = (03</'4| V'^'\A) |0102)f |3^4 

= (03041 V(^)(A)|0102)4^,_^^,^_,^, 



304| V^^'^A) |0i02)h = (03</'4|V^^'HA) 10102) 



I 1^2 
f l3^4 



(03041 V^(')(A) 10102) l^^l-^^^^^ 



64 



(03041 V^^'nA)|0102)i = (03041 V^^'HA) 10102) 



h|3^4 



(03041 V^'\A) |0i02)f . (4.34) 



The Complete Matrix Element and Transverse Locality 



Using Eqs. (|3.2|) and ( |3.31| ), we can write the full second-order matrix element of 
the IMO: 

(03041 M\A) 10102) = ( 5i,3 (^2,4 + (5i,4 §2,3 )Mj, 

+ e-^-'^^fE (03041 V^(')(A) 10102),. (4.35) 

By inspection, it is clear that this result respects all the restrictions we have placed 
on the IMO, except perhaps the requirement that it can be represented as a power 
series in transverse momenta with an infinite radius of convergence. To see how to 
check this, consider the first term in the sum on the RHS: 



A 



X e 



-a-HmI-m;)' + ^) (l - . (4.36) 



The presence of the mass denominators might seem likely to cause a problem with the 
convergence of a power series expansion. Everything multiplying the delta function 
can be expressed as a power series in transverse momenta if it can be expressed as a 
power series in Mp and Mf, since it is a function only of these and since Mp and Mf 
are themselves power series in transverse momenta. We can write 

e-^"^-^7f (03041 V^'\A) 10102), = 327r'P+5('\p,+p2-ps~p,)glj^ 



X 



^(M| + Mf) - A(M^ + M^) + -l^jl (M^M]^ + M'pMf: 



- (m|m2 + m|m|)| + 



(4.37) 
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(a) (b) (c) 



Figure 4.3: A diagrammatic representation of the contributions to (020304| 5V^'^^ 

a matrix element of the second-order change to the reduced interaction. The numbers 

label the particles. 



If we approximate the RHS of this equation by keeping only terms, the differ- 
ence between this finite sum and the RHS of Eq. (|4.36|) can be made arbitrarily 



small, for any values of Mp and Mf, by making large enough. This means that 
g-A *Al,jg2 (1030^1 V"*^^^(A) |0i02)a Can indeed be represented as a power series in trans- 
verse momenta with an infinite radius of convergence. The other contributions to 
(0304| A^^(A) 10102) can be analyzed similarly, and the conclusion is the same for 
each. 

4.1.3 Example: (0203041 A<^(A) |0i) 

The final second-order matrix element that we compute is (020304| A/i^(A) |0i). 
To compute this matrix element, we use |/) = |0i) and \F) = 1020304)- The matrix 
element of the second-order change to the reduced interaction is 

(0203041 SV^'^ |0i) = W ('^203041 V^'^ \K) {K\ V^'^ |0i) "^'^ (F, K, I). (4.38) 

^ K 

Since V^'-^^ changes particle number by 1, \K) has to be a two-particle state. 

We define (0203041 '^V"^^-' |0i)^ to be the contribution to (020304 1 <5V"^^-' |0i) from 
the process shown in Figure 4.3s, where s = a, b, c. Then 
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(0203041 5V^^^ |0l) = E (0203041 8V^^'^ |0l), • (4.39) 

Eq. ( [4.1| ) suggests that we define a contribution to the matrix element of the second- 
order reduced interaction for each contribution to (020304 1 5V^'^'^ |0i): 



(0203041 l^(')(A)|0l), = (0203041 5^^'M0l).L™ ' (4-40) 

and then 

(0203041 V(^HA) |0l) = E(<^20304|V^^')(A)|0l),. (4.41) 

For each contribution to (020304 1 SV^"^^ |0i), we use \K) = |0506) and the change 
of variables 

P3 = iypt,ypi± + Qi.)^ 

P4 = izpf,zpi± + w±), (4.42) 

where 



X + y + z = 1, 
f± + q± + w^ = 0. (4.43) 



Then the relevant differences of free masses are 

A, 



^FK 





y 




X 


X + y 


= 


( ' 




\x + y 


X 


+ — + 
y 


z 



A„ = ^ + ^ + ^, (4.44) 
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and the contribution to the matrix element of the second-order change to the reduced 
interaction shown in Figure 4.3a is 

(0203041 5V^^^ |0i), = 327r V^^'Hpi - P2 -Pz -P4)^ - ^) 

X-\-y \L^FK ^KI ) 

and the corresponding contribution to the matrix element of the second-order reduced 
interaction is 

1/1 1 \ 



(0203041 V^^'HA) I0l)a = 327r ^^^^^Pl " P2 " PS " P4) " 



^^V\^FK ^Kl) 

X [e^^'-^^FK^Ki - e^^~^^FK^Ki\ . (4.46) 

^ /A terms 

As before, if we expand everything multiplying the delta function in powers of 
transverse momenta, the lowest-order terms from the two exponentials cancel and 
leave two types of terms: those that are proportional to inverse powers of A and 
those that are proportional to inverse powers of A'. We can isolate the terms that 
are proportional to inverse powers of A without altering the cancellation of the lowest 
term by replacing the first exponential with a 1: 

(0203041 l^('nA)|0i), = 327r^p+5(5)(pi-p2-p3-P4)^f-r^--r^') 

X -\-y \L^FK ^KI ) 

X (x-e^^-'^FK^Kiy (4.47) 

According to Figure 4.3, the other parts of the matrix element of the second-order 
reduced interaction are given by 

(0203041 ^('HA)|0l)b = (0203041 ^('HA)|0i)J3^, 

= (0203041 W^HA)|0i)J^_^.^_.^, (4.48) 

and 
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= {hhMV^'\A)\<Pi)S-r'l'^'-'''^f . ..(4.49) 

Using the definition of tlie IMO in terms of tfie reduced interaction, we can write the 
full second-order matrix element of the IMO: 

(0203041 ^'(A)|0i) = e-^'*%ff: (0203041 V^('HA)|0i),. (4.50) 

s=a 

This result is consistent with the restrictions that we have placed on the IMO. 

4.2 The Removal of Cutoff Dependence from Physical Quan- 
tities 

One of the restrictions on the IMO is that it has to produce cutoff-independent 
physical quantities. To see how this requirement is fulfilled, consider as an example 
the scattering cross section for 0i02 — > 0304. At second order, the T matrix has 
contributions from the s, t, and u channels. The different channels are linearly inde- 
pendent functions of the external momenta; so each contribution should individually 
have all the properties required of the full T matrix. For simplicity, we limit ourselves 
to consideration of the s (annihilation) channel. 

In our approach, the scattering matrix S is defined by 

(F| S \I) = {F\I) - 2m6{P^ - Pf) {F\ T{Pr) \I) , (4.51) 
where the T matrix is defined by ^1 



T{p-) = i7i,t(A) + i/int(A) ^_ . r(p-), (4.52) 

and ifint(A) is the interacting part of the Hamiltonian cut off by A: 
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^int(A) = (4.53) 



To second order, 



1 



^3' 



h\ H,^tiA) + ^int(A) _ , _ ^- , ■ -gint(A) 10102) . (4.54) 



If we neglect all noncanonical operators, the first term does not contribute and the 
second term can be calculated by inserting a complete set of free states. This leads 
to 

SiPf - Pf) {<f),<p,\T{Pf) 10102) = ^e-2^"'^'647r5(5(6)(pi +P2 -P3 -P4), (4.55) 

where |/) = |0i02) and \F) = |0304). The cross section is cutoff-independent only if 
the T matrix is cutoff- independent. Thus Eq. ( [4.55|) shows that the cutoff appears 



in physical quantities if noncanonical operators are neglected. 

To test whether our procedure for computing A^^(A) leads to cutoff-independent 
physical quantities, we now include noncanonical operators. This means that we have 
to include the first term in Eq. ( [4.541 ): 



(03041 i/int(A) 10102) = ^^e-^-'^-(7^2 ^03041 K(2)(A) 10102), 

Pi +P2 

= ^ (l-e-2^-'^^')647r^5(^)(pi +P2 -P3 -P4),(4.56) 

where the result for the annihilation contribution to the second-order reduced inter- 
action is used, and we use the fact that we only need the on-energy-shell part. This 
contribution shifts the T matrix so that the total T matrix is 

6{Pp-Pr) (03041 T(P7) 10102) = ^6471^5(6) (pi +P2 -P3 -P4). (4.57) 
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Since g"^ is cutoff-independent at second order [dg'^/dA = 0{g^)], tliis is tlie correct 
cutoff-independent result to Oi^g"^). 

Our restrictions on tlie IMO liave led to its matrix elements being such that they 
compensate for the presence of the cutoff in physical quantities. No proof exists, but 
we expect that our procedure leads to a renormalized IMO that produces exactly 
correct perturbative scattering amplitudes order- by-order. When this IMO is used 
nonperturbatively, however, there will be some cutoff dependence in the resulting 
physical quantities, as we have mentioned in previous chapters. 

4.3 (02031 A42(A) to Third Order 

We wish to calculate (02031 A^^(A) |0i) to third order to demonstrate how a 
higher-order calculation proceeds and to derive the cutoff dependence of the cou- 
pling in our approach. In Chapter ^ we postulated that matrix elements of even- 
order reduced interactions cannot contain three-point interactions. This means that 
(0203| A^^(A) \(j)i) has only odd-order contributions. The first-order contribution is 
defined in terms of the first-order reduced interaction, which we have defined in Eq. 
( p. 371 ). The third-order contribution is defined in terms of the third-order reduced 



interaction. In this section, we compute the third-order contribution, neglecting the 



cutoff-independent part of the third-order reduced interaction. Eq. ( p.75|) shows that 
a calculation of the cutoff-independent part would be a fifth-order calculation, which 
we avoid at this time. 

According to Eq. ( p.66| ), the cutoff-dependent part of the third-order reduced 
interaction is given by 



,03| Vi?i{K) |0i) = [(02031 5V^^^ |0l) - 5l,2 (02031 ^^'^ l'^l)]Aterms " ^^"^^^ 
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(02031 SV^^^ |0i) is defined by 

(02031 <5r(=^M0i) = Wi<p2h\V^'^\K){K\V^'\A')\<P,)Ti'''^'^ 



K 



{A,A') 



^ K 

+ ^ E (</'203| V^^^ \K) {K\ \L) (L| |0i) rf^') (4.59) 

where we are suppressing the dependence of the cutoff functions on the states. Note 
that i?i 2 = C3 and is given by 



<oA'K^pt5^^\pi -P2 



P3) 



(02031 |0l) 



P2X=P3±;P2 =P3 



(02031 |0l)]" (02031 5V^'^ m,,,^,^,.,,.^,. 



(4.60) 



Eq. ( |4.58|) then becomes 



(02031 Vg^iA) |0i) = (02031 5V('^ |0i) - (02031 5V^'^ \<l)l)L_^.^^.._. 



.(4.61) 

A terms 



This equation indicates that (0203 1 ^cd(^) I0i) be computed solely in terms of 
{Msl^V^^^ |0i), which we now calculate. 

All the matrix elements of the second-order reduced interaction that can appear 
on the RHS of Eq. ( |4.59| ) were calculated in Section 4.1. When (0203| 5^*-'^^ |0i) 
is expanded in terms of the possible intermediate states and possible contributions 
to matrix elements of V^^^ and V^'^\A'), we see that it has contributions from a 
number of different processes. We define (0203 1 5V^^^ |0i)^ to be the contribution to 
(02031 5 V^(^) |0i) from the process shown in Figure 4.4s, where s = a, b, c, . . . , o. In 
Figure 4.4, the plain vertices represent matrix elements of V^^^ and the boxes represent 
contributions to matrix elements of V^'^^ (A') corresponding to the processes depicted 
in the boxes. For example, the box in Figure 4.41 represents (0203 1 V^'^\A') 10405)^ 
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(a) (b) (c) 




(m) (n) (o) 



Figure 4.4: A diagrammatic representation of the contributions to (<^2<^3 1 <5^*^^^ l^i)- 
The numbers label the particles, the bare vertices represent matrix elements of V^^\ 
and the boxes represent contributions to matrix elements of V^'^\A') corresponding 
to the processes depicted in the boxes. 
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[see Figure 4.2f and Eq. ( 4.33 )]. Note that not every contribution to every ma- 



trix element of V"*-^''(A') that appears in Eq. ( ^.59| ) appears in Fig 4.4. For ex- 



ample, (02031 V^'-^n^') I0405)i [see Figure 4.2i and Eq. ( |4.34| )] does not appear in 



Figure 4.4, because when it is paired with (0405 1 V^^^ |0i) and we integrate over 
and J95, its contribution is identical to that of (0203| V^'^\A') |0405)f. For this rea- 
son, (0203| SV^^^ |0i)i5 shown in Figure 4.41, includes contributions from both. Then 

(02031 5^(3) |0i) is 

(02031 SV^^^ |0l) = E (<^203| SV^^^ |0l), . (4.62) 

For each contribution to (0203 1 5V^^^ |0i), we use the change of variables 

P2 = iypt,ypi± + Q±), 

P3 = ([i-2/]p^[i-2/]pi±-gl), (4.63) 

which leads to the initial and final free masses 

Mj = 0, 

4.3.1 The One-Particle Self-Energy Contributions 

The contributions to the matrix element of the third-order change to the reduced 
interaction from one-particle self-energies, shown in Figures 4.4a, 4.4b, and 4.4c, can 
be combined because they have similar structures. For each of them we use 

1^) = 10405), 

\K) = 106), (4.65) 

and the change of variables 
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P4 = {xPi,Xpi± + f±), 

P5 = ([1 - x]pt, [1 - x]pi± - fj_), (4.66) 



which leads to the intermediate free masses 

Mi = 0, 

and Pe — Pi- The contribution from Figure 4.4a is 

(02031 SV^'^ = ^ 1 D^D,De (02031 V^'^ |06) (06| l^(') 10405) (0405| ^^'^ |0l) 

X Tf'^')(F,i^,L,/) 

= ^fd'r^dx—^Tt'^'\F,K,L,I), (4.68) 
ZTT^ J a;(l — a;) 

where 

N^^ptS^'\pi-p2-P3)- (4.69) 
The contribution from Figure 4.4b is 

(02031 sv^'^ |0i), ^Ifo, (02031 1^(^) 106) (06| v^'\A') |0i) ri^'^')(F, /r, /) 

.^A'2t(^>^\f,K,I), (4.70) 



2 

iV 



12 



and the contribution from Figure 4.4c is 



(02031 SV^'^ |0i), = ^ 1 Z)4i?5 (02031 X^^'HA') IM^), (0405| -^"^'^ |0l) 

X Ti^^^'^F, L, I) 
N 



J x{l - x) \Ml Ml) ^ ^ 



27r2 

X ri^'^')(F,L,7). (4.71) 
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It is easiest to combine these three contributions before doing the integrals. This 
leads to 



j:{<P2h\SV^'^\<P,)^ = -l d{Ml)le- 



■2A'-4m4 



1 



+ 



+ 



Ml Ml Ml 



g2A'-4M£M|, 



1 



1 



Ml 



1 



Ml MlAp^ 



Ml 



- [A' - A] 



(4.72) 



This one- dimensional integral cannot be done analytically for arbitrary Ml; so we 
consider it to be a function of Ml that must be computed numerically. 

4.3.2 The Two-Particle Self-Energy Contributions 

The contributions to the matrix element of the third-order change to the reduced 
interaction from two-particle self-energies, shown in Figures 4.4d, 4.4e, and 4.4f, can 
be combined. For each of these we use 



1^) = 10405), 
\K) = 1060708), 



(4.73) 



and the change of variables 



P6 = {xpt,xp2± + r±), 

P7 = {[i- x]pt,[^- x]P2±-r±), 



(4.74) 



which leads to the intermediate free masses 



Ml = 








yy 










y) 



+ 



fl 



(4.75) 



and 
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P4 = P2, 



P5 = P8 = P3- (4.76) 
Then these contributions are given by 



(1 - x)y'^ 

KL Ml) 



{cl>,h\SV('^ = ^y|A'2rf'^')(F,L,7). (4.77) 

Writing their sum as an integral leads to 



f AT foo ( 

J:{M3\SV^'^\<I>,)^ = -1 d(A^J e 



Ml 2 



+ g-2A'-4A^^M|, f \ ^_ ^ ^KL 



^ (02031 (4.78) 



Prom Figure 4.4, 

E (02031 51^^'^ = i: (02031 5y(^M0l).|2^3' (^'^S) 

s=g s=d 

but since Z^s=d (0203| SV^^^ \4>i)s ^ function only of M|, and N, which are invariant 
under 2 3, 

E (02031 ^V^^'M^l). = E (02031 ^V^^'M^l).- (4.80) 

s=g s=d 
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4.3.3 The Exchange Contributions 

For the contributions to the matrix element of the third-order change to the re- 
duced interaction from exchange processes, shown in Figures 4.4j, 4.4k, and 4.41, we 
use 



\K) = 1060708), 



(4.81) 



and the change of variables 



P4 = {xpt,xpi± + r±), 

P5 = {[l-x]p'l,[l-x]pi±-f±), 



(4.82) 



which leads to the intermediate free masses 



Ml 



and 



Mi 



x{l — x) ' 

Q± ^ ^± ^ {r± - q±) 
y (l-x) x-y 



(4.83) 



Pe = P2, 

P7 = ([x - y]pt, [x - y]pi± + [f± - q±]), 
P8 = P5- 



(4.84) 



Then 



1 N r 1 



X 



{l-x){x-y) 

2A'-*A^^Ml 



Ml 
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1 


1 \ 


1\ Li 


Ml ) 


1 


1 




^FK 


1 


1 


^FK 





- [A'^A]|. (4.85) 
This five- dimensional integral also cannot be done analytically. From Figure 4.4, 

s=m s=j 

= E(</>203|5^(^M0i).U_i„,;,,_,,. (4.86) 

4.3.4 The Cutoff Dependence of the Coupling 

We now have the complete matrix element of the third-order change to the re- 
duced interaction in terms of numerical functions of q±, the transverse momentum of 
particle 2 in the center-of-mass frame, and y, the fraction of the total longitudinal 
momentum that is carried by particle 2. We first use the matrix element to compute 
the cutoff dependence of g^. According Eq. (f4.60|), we need (0203 1 5V^^^ |0i) in the 



limit p2 = pt and p2± = P3±. In this limit, q± = Mp = 0, and all the integrals in 
(02031 5V^^^ |0i) can be done analytically. However, care must be used when taking 
this limit because there are delicate cancellations of divergences among the different 
parts of each integrand. The result of applying Eq. (|4.6CI| ) is 



which from Eq. ( |3.49| ) implies that 



9. = 9a, + ^,gl log ^ + 0{gl). (4. 
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This equation tells us how the couplings at two different cutoffs are related. Since 



A' > A, Eq. ( ^4.88| ) shows that > (?^, . This means the coupling grows as we 
reduce the amount by which free masses can change. This shows that the theory is 
asymptotically free, as expected. 

In conventional covariant perturbation theory with the minimal subtraction (MS) 
renormalization scheme, one obtains, for the scale dependence of the coupling ||28| , 

9, = V + ^^9l' log ^ + 0{gl,), (4.89) 

where fj, and /x' are two different renormalization points. We see that our coupling 
changes with our cutoff in the same manner that the MS coupling changes with the 
renormalization point. This is what one expects due to the scheme independence of 
the one-loop beta function. However, since there is no simple connection between our 
scheme and the MS scheme, we expect there to be no simple relation between the 
respective couplings at higher orders (although we could derive the relationship by 
equating physical quantities calculated in the two schemes). 

4.3.5 Numerical Results for the Three-Point Matrix Element 

We now proceed to calculate the matrix element of the cutoff-dependent part of 
the third-order reduced interaction. From Eq. (4.61), 



(02031 V^S(A) |0l) = [(02031 SV^'^ |0l) - (02031 SV^'^ l^l) l.-^J ' (^-SO) 



A terms 



where we use the fact that the gl = limit of (0203 1 5V^^^ |0i) is independent ofp2 and 
p^. The IMO is invariant under rotations about the 3-axis and q± is the only trans- 
verse vector on which (0203| 5V^^^ |0i) depends. This means that (0203| 5V^^^ |0i) is 
a function of ql- By inspection, we can see that (02031 SV^^^ |0i) is the product of 
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and some dimensionless quantity, and is the difference of a function of A and the same 
function with A — > A'. Also, by our assumption of approximate transverse locahty, 
the RHS of Eq. ( |4.90| ) can be expanded in powers of gl. Thus 

oo 

(02031 sv'^'^ |0i) - (02031 sv^'^ |0i)| .^^^ = NY.My)q± {^~'' - a'"^^) , (4.9i) 

i=l 

where the /j's are unknown dimensionless functions of y. This implies that 

(02031 V^(5'l](A)|0i) = [(0203|5V('M0l)- (02031 5V^^'M0l)L-.,-ol,, 

L ' y±— 'JJ A terms 

oo 
i=l 

= U<l>2<l>3\SV^'^\<l^i)-{hh\SV^'^\<jy,)\ ^ .(4.92) 

Using this equation and our result for (0203 1 5V^^^ |0i), it is straightforward to write 
(0203| Vcd(A-) |0i) as a numerical function of q±, y, and A. 

It is clear by inspection that (0203| SV^^^ |0i) is symmetric under y 1 — y and 
each (02031 |0i), for s < i is a function only of Mp, the free mass of the final 
state. It is not clear from inspection if the rest of (0203 1 SV^^^ |0i) is a function only 
of Mp and there is no reason that it should be. This means that the matrix element 
of the IMO is not necessarily a function only of Mp. 

Our result for the three-point matrix element of the IMO through third order, 
neglecting any cutoff-independent contribution to the third-order reduced interaction. 



IS 

(02031 M'iA) |0i) = e-^"^^ [647r5p+5(5)(pi -^2-^3)^?^ 

+ (02031 ^(5S(A)|0i)]. (4.93) 

To get an idea of the size of the noncanonical contribution to the matrix element, 
we can compare it to the canonical contribution. To do this, we need to choose a 
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value for the coupling. We would like to choose a large coupling so that we can get 
a pessimistic estimate of whether or not the expansion for the IMO is converging. 
When the coupling is large, the second term in Eq. ( [4.88| ) is nearly as large as the 
first. The value of log(A'^/A^) is our choice, but the natural value is 1, because then 
the range of scales over which off-diagonal matrix elements of the IMO are being 
removed is comparable to the range of scales that remain. Thus we estimate that a 
large coupling is g"^ = 5127r^/3. 

Using this coupling. Figure 4.5 compares the regulated canonical and noncanonical 
contributions to the matrix element as functions of the free mass of the final state. The 
plot is the result of a numerical computation using the VEGAS Monte Carlo algorithm 



for multidimensional integration |2^. The contributions to the matrix element are 
plotted in units of the matrix element of the unregulated canonical interaction, and 
we consider only the situation in which the two final-state particles share the total 
longitudinal momentum equally. The solid curve represents the canonical contribution 
and the diamonds represent the noncanonical contribution. The statistical error bars 
from the Monte Carlo integration are too small to be visible. 

Figure 4.5 shows that for the case we consider, the noncanonical corrections to 
the matrix element of the IMO are small compared to the canonical part, even when 
the perturbative expansion of the coupling is breaking down. However, this does 
not necessarily imply that corrections to physical quantities from the noncanonical 
part of this matrix element would be small. An investigation of this would require a 
fifth-order calculation in this theory; so we do not pursue it. 

For the same coupling. Figure 4.6 shows the noncanonical part of the matrix 
element as a function of the fraction of the total longitudinal momentum that is 
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carried by particle 2 for fixed free mass of the final state. It is plotted in units of the 
matrix element of the unregulated canonical interaction for three different values of 
the free mass of the final state. This plot shows that the noncanonical part of the 
matrix element is a function of y for fixed Mp, and thus is not a function only of Mp. 
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Figure 4.5: A comparison of the regulated canonical and noncanonical contributions 
to (0203 1 A^^(A) |0i) as functions of the free mass of the final state. The contributions 
to the matrix element are plotted in units of the matrix element of the unregulated 
canonical interaction, and we consider only the situation in which the two final-state 
particles share the total longitudinal momentum equally. The coupling used in the 
plot is = 5127r^/3. The solid curve represents the canonical contribution and the 
diamonds represent the noncanonical contribution. The statistical error bars from 
the Monte Carlo integration are too small to be visible. 
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Figure 4.6: The noncanonical part of (0203| A^^(A) as a function of the fraction 
of the total longitudinal momentum that is carried by particle 2 for fixed free mass 
of the final state. It is plotted in units of the matrix element of the unregulated 
canonical interaction for three different values of the free mass of the final state. The 
couphng used in the plot is = 5127r^/3. The statistical error bars from the Monte 
Carlo integration are too small to be visible. 
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CHAPTER 5 



CALCULATION OF GLUEBALL STATES AND MASSES 

IN PURE-GLUE QCD 



In this chapter we calculate part of the pure- glue QCD IMO and diagonalize it, 
obtaining glueball states and masses. We work in four dimensions explicitly and begin 
by defining a basis in terms of which we expand all physical states. We construct 
the physical states to be simultaneous eigenstates of the IMO, the three- momentum 
operator, and the part of the generator of rotations about the 3-axis that generates 
rotations of the internal degrees of freedom of states (as opposed to their centers of 
mass). We approximate all physical states as two-gluon states, which means that 
they are relatively simple single-glueball states. We analytically calculate the color 
parts of the states by requiring the states to be color singlets, and then we expand 
the states' spin and momentum degrees of freedom in a complete set of orthonormal 
basis functions. After deriving the eigenvalue equation for the IMO in this basis, we 
use the recursion relations for the IMO to compute its two-gluon to two-gluon matrix 
element to second order in the running coupling. This is the only free-state matrix 
element that we need to solve the eigenvalue equation. 

We compute the matrix elements of the IMO in our basis, in terms of integrals 
that must be evaluated numerically. To avoid roundoff error, we evaluate the integrals 
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for the kinetic energy and the two-point interaction (the self-energy) by writing them 
as sums of gamma functions. The remaining integrals are five-dimensional, and we 
compute them with Monte Carlo methods. However, before we can do this, it is 
necessary for us to make manifest the cancellation of the infrared divergences from 
exchanged gluons with infinitesimal longitudinal momentum. This is nontrivial, and 
even after the divergences are gone, we must manipulate the integrals quite a bit to 
get them into a form that is amenable to Monte Carlo integration. We present some 
of these details and other technical issues in Appendix C. 

We proceed by discussing how we assign quantum numbers to our numerical results 
for glueball states and the procedure that we use to compute these results. We 
derive the nonperturbative cutoff dependence of the coupling and discuss the cutoff 
dependence of our glueball masses. We use this analysis to choose the value of the 
cutoff that minimizes our errors. We then present the spectrum that we find with 
this optimal cutoff and compare it to recent lattice results. The last results that we 
present are the probabihty densities for our glueballs. We conclude this chapter by 
discussing the sources of error in the calculation and estimating the uncertainty in 
our spectrum from these errors. 

5.1 The Basis for the Expansion of Physical States 
5.1.1 Preliminaries 

The physical states of pure-glue QCD are eigenstates of A^^(A). Since our cutoff 
preserves translational covariance and covariance under rotations about the 3-axis, 
we would like the states to also be simultaneous eigenstates of the generators of 
these symmetries, but this is impossible because translations do not commute with 
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rotations. However, a rotation about the 3-axis separates into a part that rotates the 
centers of mass of states and a part that rotates states' internal degrees of freedom, 
and translations do commute with these internal rotations. 

To be precise, we recall from Chapter 2 that Jz is the generator of rotations about 
the 3-axis, and we define to be the part of J-^, governing gluons' momenta in the 
center-of-mass frame and spin polarizations. Al^(A), 7^+, Pj,, and Jf^ are a set of 
commuting observables; so an eigenstate of A4^(A) can be labeled by their eigenvalues. 
We choose to write a physical state as |^-'"(P)), where P is the three-momentum of 
the state, j is the eigenvalue of Jf^ for the state, and n labels the mass eigenvalue 
of the state {n — \ has the smallest mass, n — 1 has the second-smallest mass, etc.). 
Note that because \W'^{P)) will be determined by 7W^(A), it will implicitly depend 
on A and g^^. 

An examination of the matrix elements of A4^(A) leads us to believe that the 
light-front infrared divergences (see the discussion below) will not cancel unless the 
physical states are color singlets, although we do not have a rigorous proof of this. 
Therefore, we assume that the physical states are color singlets. Using the Fock-space 
expansion for the identity operator, we can expand a physical state in terms of the 
number of gluons: 



where there is no one-gluon component because there is no color-singlet gluon. We 
neglect contributions to the states with more than two gluons, which is a severe 
approximation. This means that we are approximating all physical states as relatively 
simple single-glueball states, and j is then the projection of the glueball's spin onto the 
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3-axis. From now on, we refer to the approximate physical states simply as glueball 
states. 

Any eigenstate of ^A'^{A) can be written as a superposition of cross-product states: 

= ElTrn.)®|rr(n), (5.2) 

i 

where |T„ j) is the color part of the i^^ contribution to the state and r^"(P)^ is the 
momentum/spin part. 

The free state \g1g2) is defined as a two-gluon state in which one particle has 
quantum numbers pi, Si, and Ci, and the other has quantum numbers p2, S2, and C2. 
\g1g2) is symmetrized with respect to particle labels. Thus we can write 

1 



\g1g2) 



V2 



1^1; ^2) + 1^2; ^1) 



(5.3) 



where the semi-colon in \gi; g2) indicates that particle #1 has quantum numbers pi, 
si, and ci, and particle #2 has quantum numbers p2, S2, and C2. Then we can write 
the wave function as the sum of two permutations: 



gig2 



Vl>^'^(P)) = —[{g,;g,\ + {g,;g,\\^^-{P)y 



(5.4) 



This implies that the wave function {gig2 |\1/''"(P)) is symmetric under exchange of 
particles 1 and 2. Using Eq. ( p.2|) , we can write the first permutation as 



5-1; 5-2 



^'"(^)) = E(ci;c2|T„,)(pi,Si;p2,S2|rr(P)). (5.5) 

i 

(ci; C2 |T„ j) and (pi, Si;p2, S2 T{"'{P)'^ are the color and momentum/spin wave func- 
tions for the first permutation of the i^^ contribution to the total wave function. 
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5.1.2 Calculation of the Color Wave Function 

In our two-gluon approximation, the fact that the glueball states are color singlets 
uniquely determines their color wave functions. Using the completeness relation for 
color states, the color part of the i*^ contribution to the glueball state can be written 

l^n.i) = XI l^i^^a) (ci;c2 |T„,j) . (5.6) 

C\C2 

Let Y be the unitary operator that rotates two-gluon states in color space. Then 
under a color rotation, we have 



|T„,i) T^^i) = |ci;c2) (ci;c2|T„,i) 

ClC2 

= |c'i; 4) (ci; C2 |T„,i) , (5.7) 



ClC2 



where the primes on the basis vectors denote rotated basis states. Since |T„^j) is a 
color singlet, it does not change under a color rotation: 



\^n,i) = ^n,i} = 4) (Ci; C2 |T„,,) . (5.^ 

ClC2 

A gluon color-rotated state is given by B 



= U |Ci) = Y |C2) (C2| U |Ci) = Y Uc,c, |C2) , (5.9) 



C2 C2 

where 

U = e-''^^^^' (5.10) 

and 

F^u = -^r"''"'- (5.11) 

Here 6c is a vector of A*"^ — 1 real numbers that parameterize the color rotation. Since 
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14; 4) = \c[) ® 14) 



(5.12) 



Eqs. (|]8|) and (U) imply that 



C1C2C3C4 



As a guess, let us try 



(Ci;C2 |T„^j) — — 5ci,C2; 



where is some constant. Then Eq. (|5.13|) becomes 



Note that 



^ ] Uc^dUc^d 1^3, C4) . 
C1C3C4 



(5.13) 



(5.14) 



(5.15) 



(5.16) 



Then 



C3C4 



E (ci; C2 |T„,j) |ci; C2) 



C1C2 



IT 



(5.17) 



Since this is true, our guess for a color-singlet wave function was correct. This solution 



is unique [£T|; so we can drop the subscript i in Eq. ( p.5| ). is determined by 
normalization to be \J N"^ — 1. 



Now Eq. (|5.4| ) becomes 
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9i92 



2(iV2 - 1) 



+ (P2,S2;Pl,Si 



-'Cl,C2 



Pi, Si;p2, S2 



(5.18) 



Since the IMO conserves momentum, the factor in brackets must be proportional to 
a momentum-conserving deha function. Thus we can write 

1 



2(iV2 - 1) 



-.6r 



Pi Pi 



(5.19) 



where $^"^^(^1,^2) is the momentum/spin wave function with the momentum-conser- 
ving delta function removed, and all the extra factors in this equation are present to 
simplify the normalization of $^"^^(^1,^2)- We must solve for $^"^^(^1,^2), and this 
equation indicates that it is symmetric under exchange of particles 1 and 2. 

5.1.3 Jacobi Variables 



Using Eq. (|5.19|) , the Fock-state expansion for a glueball state in Eq. ( 5.1 
becomes 

1 1 ^ . r d?pi±dpi 



^J'^(P) 



VlQ^ ^2{N^ - I) 



- I) S1S2C1C2 JpiP2 



(5.20) 



where momentum conservation implies that P2 = P ~ Pi- It is useful to separate the 
motion of the center of mass of the state from the internal motions of the gluons. To 
do this, we change variables from pi to the Jacobi variables x and k±: 



Pi = {xP^,xP± + k±), 

P2 = i[l-x]P+,[l-x]P^-kj_). 



(5.21) 
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Here x is the fraction of the total longitudinal momentum that is carried by particle 
1, and k± is the transverse momentum of particle 1 in the center-of-mass frame. We 
only display the longitudinal and transverse components of the momenta. (Since the 
glueball state is a superposition of free-particle states, and since the momentum of a 
free gluon satisfies p"^ = 0, the minus components of the momenta of gi and g2 are 
constrained to be given by pj^ = pf^/pt ■) 

In terms of the Jacobi variables, the glueball state is 

^"(^)> = ft 3 / ^ S 6c,,J^ML=e{x-e)e{l-x-e) 

X 



^l^,s^{x. P) 9{x, k^, P; si, ci)g{l - X, -k^, P; Ss, Cs)) , (5.22) 

where we explicitly show the dependence of the RHS ket on the Jacobi variables 
and the total momentum. Recall from Chapter 2 that each particle's momentum 
transforms the same way under a boost, regardless of whether the particle is point- 
like or composite. We can use this to show that ^i^^^i^x, k±, P) is independent of P. 
To do this, we apply a boost operator to both sides of Eq. ( |5.22|) : 



^^P')) = ,1 , ,\ , E 5c^,cJ ^Ml=e{x - e)e{l - X - e) 
I V16^y^2(Ar2 - 1) sis1^ic2 i ^x{l - x) 

X ^l^^^^{x,k^,P) g{x,k^,P'-SuC,)g{l-x,-k^,P'-s,,C2)), (5.23) 



where the boost takes the glueball's momentum from P to P'. Note that the boost 
does not affect the wave function, only the kets. Since Eq. ( p.22| ) holds for all P, it 
holds in particular for P': 

^^"(n) = ,t , ,\ , E 5c,,cJ ^Ml=e{x - e)e{l - X - e) 
' p(N^ - I) s^t;^,c2 J ^x{l - x) 

X 



$i;^,^(x,fc^,P') ^(x,fc^,P';si,Ci)^(l-x,-fc^,P';s2,C2)). (5.24) 
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Eqs. ( |5.23|) and (|5.24|) contradict each other unless the wave function ^^"^^(x, k±, P) 



is independent of P: 

^i-^,^{x,k^,P) = ^Z^{x,k^). (5.25) 
Thus we can write |\1/-'"(P)) as 

X , (5.26) 

where 6*^ = 6'(x — e)^^(l — x — t). 

5.1.4 The Momentum and Spin Wave-Function Bases 

— * 

To solve for ^^"^^(x, fc^), we expand it in a complete orthonormal basis for each 
degree of freedom. Since ^^"^^(x, k^ is symmetric under exchange of particles 1 and 
2, 

$ir.,(x,M = $i:.,(l-x,-^^). (5.27) 

To use this, we define 

9=1 

where the spin wave functions are 

Xl — fsi,l^S2,l5 

X2^''^ = (^51,1^52,15 

= ^ [^^i.K^Sa.l + '^Si,l'5.2,l] , 

xr^^ = ^ [^.i,i'5s-2,i - '5s-,i5s2,i] , (5.29) 
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where s = —s, and the momentum wave functions satisfy 



Qr{x,k^) = -(]f (5.30) 



Note that 



Exr"xr = V- (5.31) 

SlS2 

We define k = and we define the angle by 

fc_L = cos (/) X + sin ^. (5.32) 
We expand the momentum wave function in complete orthonormal bases: 

oo oo oo 

^"(^,^±) = EE E B';taL\''\^)T!'\k)AM), (5.33) 

1=0 t=0 a=-oo 

where Li'\x), Tl'^\k), and A^(0) are basis functions for the longitudinal, transverse- 
magnitude, and transverse-angular degrees of freedom, e and d are parameters that 
govern the widths of the longitudinal and transverse-magnitude basis functions, re- 
spectively. We can adjust these widths to optimize the bases. Note that if we do 



not truncate the sums in Eq. (|5.33| ), then the RgitaS depend on e and d such that 



Q^g"'{x,k±) is independent of e and d, although we do not indicate this dependence 
explicitly. 

We define the transverse-magnitude basis functions Tj;'^\k) by 

Ti'\k) = de-^'''Y.ot,sk'd\ (5.34) 
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where the cxj^s's are constants that we have computed numerically using the Gram- 
Schmidt orthogonalization procedure and are such that 

°° dkkTt^'^\k)Tif\k) = 6t,t'. (5.35) 

Adjusting d allows us to adjust the width of the Gaussian weight function in Eq. 
( p.34| ). d cannot be allowed to pass through zero, because when d is zero, Tj;'^\k) is 



zero. We choose d to be positive without loss of generality. When doing numerical 
computations, we work with a dimensionless form of these basis functions: 

T,{kd) = -J^'\k). (5.36) 

Under exchange of the two particles, k is unaffected; so Tl;'^\k) is unaffected. 
We define the longitudinal basis functions Lf\x) by 

Lt\x) = [x{l-x)rj:\^t!m^-, (5.37) 

m=0 

where 



Ae)_, 1 r(l + 4e + / + m) 



= (-1)' 



+ (5,38) 



i\{l-m)\ r(l + 2e + m) ^ r(l + 4e + /) ' 

These definitions imply that 

dxLi'\x)Llf\x) = 5^,, (5.39) 



1 



as long as e > —1/2. (If e < —1/2, then the state is not normalizable.) We can 
adjust the width of the weighting function in Eq. (|5.37| ) by adjusting e. It is often 
more convenient to work with 

Lf\x)= I / lP(x). (5.40) 
'x(l — x) 
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Under exchange of the two particles, x — > 1 — which means that Lf\x) 

We define the transverse-angular basis functions Aa{(t)) by 

M^) = ^e-^ (5.41) 

and then 

/■27r 

/ d<t>Al,{<t>)A,{<t>)^5,,a'. (5.42) 

These basis functions are useful because they are eigenfunctions of C^, the part of the 
generator of rotations about the 3-axis that governs gluons' momenta in the center- 
of-mass frame. Under exchange of the two particles, — > + tt, which means that 

Aa{4>) ^ {-lYAa{4>). 

Using the above definitions of the bases, \^^'^{P)) can be written 

' Vl67r-y2(iV2 - 1) ,,,,eiC29ita ^ 

X ^„(0)|^iy2). (5.43) 

Since the glueball state is symmetric under exchange of the two gluons, the behaviors 
of the spin and momentum wave functions under exchange of the two gluons imply 
that if g = 4, then / + a must be odd, and if g 7^ 4, then / + a must be even (so that 
the spin and momentum wave functions have the same symmetry under exchange). 

5.1.5 Rotations about the 3-Axis 

We want to ensure that |^'-^"(P)) is an eigenstate of with eigenvalue j. The 
action of Jz on the state is given by 
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X 



J3 "^'""{P) 
'■ 2 d 



X ^a(0)M5l^2). 



d 



d 



d 



'2± 



+ Si + S2 



Lf{x)Ti''\k) 



(5.44) 



Using the definitions of the Jacobi variables, we can separate the center-of-mass and 
internal degrees of freedom: 

X 



^9 d „i 9 9 



± 



— «7T7 + Si + S2 



Lf\x)Ti'\k)Aa{ci>) 



X |5'i5'2) 



(5.45) 



This imphes that the action of Jf - on the glueball state is given by 

1 1 



X 



+ Si + S2 

0(f) 



Vl6^ ^2(Ar2 _ 1) 

^ 1 + Lj^)(a;)Ti'^)(A;)A.(0)| to) 

i— E <^c.c.^iLx^^^^ / [a + si + S2] i;^^ (x) 

(5.46) 



Vi6^y2(Ar 



Since |^'^"'(P)) is an eigenstate of with eigenvalue j, it must be the case that 
a + Si + S2 = j for all values of a, Si, and S2 that contribute to the sums in this 
equation. This implies that for some set of coefficients i?^^, 



Pqlta - P-M [^q,l^a,j-2 + (^g,2(^a,j+2 + 5q,2,K,j + ^qA^J 



(5.47) 



This means that [^^"(P)) can be written 
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*^"(^)) = Y.R'^t\Q,l,t,3), (5-48) 

qlt 



where 



(0) biy2> . (5.49) 



We are going to calculate the matrix elements of the IMO in the \q,l,t,j) basis 
and diagonalize it. Since A1^(A) commutes with J^, we can do this for each value of 
j separately. The diagonalization procedure will yield mass eigenvalues and the R■'^^ 
coefficients. As long as the coefficients satisfy 

^3,3' Kt'^^qlt = ^3,3'Kn' , (5-50) 
qlt 

the glueball state will have a plane-wave normalization: 

(^/-'(P')I^^-(P)) = 16^P^6^'\P - P%,fdn,n'. (5.51) 

Due to the symmetry of the glueball state under the exchange of particles 1 and 
2, the basis state |g, I, t,j) is zero if Z + j is odd and g 7^ 4, or if Z + j' is even and 
q = 4. To take advantage of this, we consider only the subspace in which l+j is even 
if g 7^ 4, and Z + j is odd if g = 4. 

5.1.6 The Eigenvalue Equation 

The IMO's eigenvalue equation is 

M^{A) ^^"(P)) = ^^"(P)) , (5.52) 
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where the lower-case subscript on indicates that it is an eigenvalue of Al^(A) 
rather than an eigenvalue of A^l-ee- (An example of an eigenvalue of A^lee ^k-) 
In the basis that we have defined, the eigenvalue equation takes the form 

M\A) E \q, I, t, j) ^MlY. \q, I, t, j) . (5.53) 

qlt qlt 

If we project (g', I', t',j\ onto the left of this equation and use the identity 

{q',l',t',j\q,l,t,j) = 167r3p+5(3)(P-P')VW> (5-54) 

then we find that the eigenvalue equation in this basis is 

^ {q',r,t',j\M^{A) \q,l,t,j) _ 2 pin 

^ 167r3P+5(3)(P-P0 - ^^-^uRq'l'f 

5.2 Calculation of the Free- State Matrix Element of the In- 
variant-Mass Operator 

5.2.1 Preliminaries 

Before we can solve the eigenvalue equation, we must calculate the matrix elements 
{q',l',t',j\A4'^{A) \q,l,t,j). These can be written in terms of the free-state matrix 
element {g'ig2\ A1^(A) \g1g2), which is specified by our renormalization procedure. We 
are going to calculate this matrix element to second order in perturbation theory. 

Before continuing, we would like to point out three simplifications that we repeat- 
edly use in this section. To make these simplifications clear, we note that the matrix 
elements {q' ,1' ,t' , j\ M.'^{h)\q,l,t, j) can be written as integrals of wave functions 
times the free-state matrix element {gig'2\ A1^(A) \gig2)- To make the first simplifica- 
tion, we observe that {g'ig2\ Al^(A) \g1g2) has step functions on the particles' longitu- 
dinal momenta, and these step functions are redundant because they also appear in 
the integrals that define (g', V , t' ,j\ Al^(A) |g, I, t,j). We drop these step functions in 
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the formulas that we present for {gig'2\ A1^(A) \gig2)- To make the second simphfica- 
tion, we point out that the matrix elements (g', t', j| M.^{h) |g, /, t, j) are symmetric 
under exchange of the two initial-state gluons and also under exchange of the two 
final-state gluons. Because of this, when we are computing {g[g2\M.'^{Iy) \g1g2), we 
combine terms that are related by exchange of the two initial-state or two final-state 
gluons. Finally, due to the color-singlet nature of the glueball states, certain parts 
of {g[g'2\ J^'^{^) \g1g2) do not contribute to {q' ,1' , i\ M.'^{K) \q,l,t,j), and we drop 
these terms. 

We begin by defining Jacobi variables for the final state of the free matrix element: 



Then using the definition of the IMO in terms of the reduced interaction in Eqs. (|3.1|) , 



(PD, and O), 

{g[g'2\ M\A) \g1g2) = {g'Al -^L \9192) + {g'Al -^int(A) \g1g2) 

= M] {g[g'2\gig2) + {g'M ^(5? \g1g2) 



([l-a:']P+,[l-x']Px-A;l). 



(5.56) 



■FI 



{9W2\V^'M)\9i92) 



(5.57) 



where the free masses of the final and initial states are given by 



x'{l — x') 



(5.58) 




5.2.2 The Cutoff-Independent Part of the Reduced Interaction 



Based on the definition of V^j in Eq. (|3.39|) , we find that 
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{g'Mvi:^ \9i92) = WP+5(^)(P-P')E^!:2;1',2' ^gW2\v\9l92-^^. (5.59) 

We can divide this matrix element into contact (momentum-independent) and instan- 
taneous-exchange interactions: 

{9'i9'2\ ^ci \9192) = {9'i92\ Vci \9i92)coN + (^^^2! ^c? \9i92)ij^ , (5.60) 

where 

{9W2\ ^c? \9i92)coN = 16vr3p+5(3)(p _ P')<2;i',2' Wi92\v\9i92:^i 
= 327r3p+5(3)(P - p')fcic[cfc',c,c ^5^^^^,Ss,,s', - 5^;,^- 5^,,.,) , (5.61) 

and 

{9W2\ \9192U = lG7i'P^6^'\P - P') E<^^2' Wl92\v\9l92:^^ 

i=2 

X ^—-{x + x'){l-x + l-x'). (5.62) 

[x — x'Y 

(The remaining term in Eq. ( p.59| ) vanishes because the \q,l,t,j) states are color 
singlets.) The contact interaction is displayed in Figure 5.1 and the instantaneous- 
exchange interaction is displayed in Figure 5.2. 

5.2.3 The Cutoff-Dependent Part of the Reduced Interaction 
Preliminaries 

According to our recursion relation for the cutoff-dependent part of the reduced 
interaction, Eq. (|3.66| ), {g[g'2\ V^^^{A) \g1g2) is given by 

(^l^^lV^S(A)te) = {9W2\SV^'^\9i92)\^,,,^,. (5.63) 
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Figure 5.1: A diagrammatic representation of (5'i5'2l lfl'i5'2)coN' two-gluon 
contact interaction. The numbers label the particles. 




Figure 5.2: A diagrammatic representation of {g[g2\ V^j \gig2)ij^, the two-gluon 
instantaneous-exchange interaction. The numbers label the particles. 



Using the definition of the 0{gj^^) change to the reduced interaction in Eqs. ( p.22|) 



and (|3.44| ), we find that 



{9W2\vi"iW\gm) = lj:(9W2\V^'^\K){K\V^'^\g,g2) 

X Tf'''^(F,i^,/)L,_. (5.64) 

The intermediate state can be either a one-particle or three-particle state. The con- 
tribution to the eigenvalue equation from the one-particle-intermediate-state part is 
zero because the \q,l,t,j) states are color singlets. This means that 

{g[g'2\ vgliA) \g,g2) = ^ / ^3^41^5 {g'Al V^'^ Igsg.g,) {gsg.g.l V^'^ \gig2) 

X rf-^')(F,ir,/)|,^_, (5.65) 

where Mj^ is the mass of the state Igsg^g^)- Substituting the definition of V^*^^^ in Eq. 
( p.38| ) into this equation, and simplifying, we find that 
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3 



1 




5 

2 TTTTTTTTTTm 2' 



Figure 5.3: A diagrammatic representation of {g'ig2\ Vcd{^) |fl'ifl'2)sE5 the self-energy 
interaction. The numbers label the particles. 



{g'Al vSliA) \9i92) = {g'Al VgiiA) \9ig2)sE + {9W2\ vg^iA) Igmhx^ (5-66) 
where the self-energy interaction is given by 

{gW2\ vg^iA) \g,g2)s^ = P+\l6^f J D;D,D,t!^''^'''\f, K, I) 62,, 52',5 Pi 
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X /^i'34L™' (5-67) 
and the exchange interaction is given by 

{gW2\vSi{A)\g^g,)^^ = 2P^\l6n'yJ D3D,D,Ti''''''\F, K, I) 62,, Sr^s 

X /^2'45L™' (5-68) 

and we have defined 

= ^{pt)0{ptnpt)S^''\p^-PJ-Pk) ■\:gjgk\vh:^. (5.69) 

The self-energy interaction is displayed in Figure 5.3 and the exchange interaction is 
displayed in Figure 5.4. 

The Self-Energy Interaction 

After some simplification, the self-energy interaction takes the form 
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Figure 5.4: A diagrammatic representation of {g'ig2 \ ^cd(A) |5'ifi'2)EX' exchange 
interaction. The numbers label the particles. 



S3S4C3C4 Ps 



1 



X -^9{pt-eP^)Ti'''^'\F,K,I) j:gs9M9i:^ ^^?>l^?3^74^ L^erms ' (5-70) 

where = pi — Ps- At this point, it is useful to change variables from p^ to Jacobi 
variables y and r^: 



Ps = iypt,ypi± + r±). 



Then the free mass of the intermediate state is 



(5.71) 



Ml 



{Ps +P^ +P5)P^ -Pi 



fc2 



+ 



a;(l — x) xy{l — y) 
and after a lot of additional simplification, 



(5.72) 



'BE 



27r2 



X 



A2_A'2 
11 



^2,2' Oix - 2e) 



log X — log e — 



12 



(5.73) 



A terms 



Recalling that "A terms" means that we are to expand the RHS of this equation 
in powers of transverse momenta and keep only the terms that are proportional to 
powers or inverse powers of A, we find that 
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{gW2\ v^im |^iff2>sE - ^/l^' '^lA' e{x - 26) 



logx-loge- ^ 



. (5.74) 



This result shows that our cutoff violates cluster decomposition. The evidence is that 
the self-energy depends on x = pt /ipt Pt)i even though p2 is the momentum of a 
spectator (see Figure 5.3). 

The Exchange Interaction 

For the exchange interaction, the free mass of the intermediate state is 



Ml = (P3 +P4"+P5)^+-^± 
k'"" (k^ - k[Y 



X' 



X — X' 



1-x 



Then the changes in free mass are 



A;'2[l - xf + - xf - 2kk'{l - x'){l - x) cos 7 
{I - x'){x - x'){l - x) ' 
k^x'"^ + /c'^x^ — 2kk'xx' cos^ 
xx'{x — x') ' 



where 7 = — 0'. 
Using the identity 

- ^ -1 .^^ 

the exchange interaction becomes 

\I\fK ^IkJ ^ i^rn=\ 

where 



(5.75) 



(5.76) 



(5.77) 



,/ Q{x — x' — e) 



x — x' 

m) 



,(5.78) 



A terms 
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g(i'^) = -^sis[6s,,s', i-ke''^^ - k'e'^'A [k' e-'''^^' [I - x] - ke-'''^^[l - x'] 

X OC I ^ I 

g(i.2) = sis'^5s,,s',^^ke''^'^ -k'e'^^'^'^^ke-'''^^ 

Q{i.3) ^ -siS2<5.-' l-ke''''*' - k'd^'A [k'e'^^'t'' - l—^ke'^A , 

^' ^ \ X 1 — X 



— /cA;' — x) + x(l — x')] COS7 >, 



2) 



1 


X 




x' 




1 




X 




x' 




1 




X 




x' 



g(2.3) = S2s'^6s,,s' Ixk'e-'"'^'^' - x'ke-''H [k'e'''^' - ^-^ke''^'*' 



— X 



In Eq. ( ^.781) , if we expand everything multiplying the delta function in powers of 



transverse momenta, the lowest-order terms from the two exponentials cancel and 
leave two types of terms: those proportional to inverse powers of A and those propor- 
tional to inverse powers of A'. We isolate the terms that are proportional to inverse 
powers of A without altering the cancellation of the lowest term by replacing the first 
exponential with a 1: 

{9i92\vt?i{^)\gig2)^^ = QA7i^P^5^'\P-P')r<n^'-j-^^9{x-x' -e) 



'^FK ^IkJ i^ra=\ 
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5.2.4 Combining the Interactions 

In order to get the infrared (e — * 0) divergences to cancel, it is useful to combine 
the interactions in a particular manner. From Eqs. ( |5.57| ), ( |5.6(J| ), and (5.66), 



{g'ig'2\ M'^{K) \g1g2) = {g[g2\ M^{A) \ + {g[g'2\M\A) {g^g^) 

CON 

+ {9W2\ M\A) \gig2)m + {9W2\ M\A) \gig2)sE 

+ {g[g'2\M\A)\g,g2)^^, (5.81) 



where 



{g[g'2\M\A)\g^g2)^^ = {g[g'^\g,g2) , 
{g[g'2\M\A)\g,g2)co^ = ^>-^"^^M^?i^72l ^c? te)coN > 
{g[g'2\M\A)\g,g2),^ = ^?,V^"^" te)iN , 



{9W2\M\A)\g,g2),^ = ^K^^^" (^?;^?2l K5S(A) te)sE , 



W2\M-\A)\g,g2)^^ = 9^'"''-' {9W2\vi1lW\9ig2h^. (5.82) 



We break the instantaneous interaction into two parts: a part that is "above" the 
cutoff, i.e. a part that would vanish if we took A 00, and a part that is "below" 
the cutoff, i.e. a part that would survive if we took A — > 00: 

{g[g'2\ M\A) \g,g2U = {9W2\ M\A) \g,g2)t^ + {g[g'2\ M\A) \g,g2)f^ , (5.83) 

where 

{9W2\ M'{A) \g,g2)t^ = (l - e-^^-^^-^-) {g[g',\ M'{A) \g,g2\ 



'IN ' 



{g[g'2\ M\A) \g,g2)'^^ = e-'^'*^--^'- {g[g'2\ M\A) \g,g2),^ . (5.84) 

Next, we break the self-energy and exchange interactions into finite and divergent 
parts: 
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{g[g'^\ M\K) \g,g2)^^ = {gW,\ -M'(A) \gig2)l^ + {g[g',\M'{A) {g.g^)^^ , 

{gW2\ M\A) |^?i(72)ex = {9W2\M'{A) \g,g2)l^ + {gW2\ M\A) \gig2)^^ , (5.85) 

where the divergent part of the self-energy interaction consists solely of the term 
containing the loge, and the divergent part of the exchange interaction consists solely 
of the term containing Q^^'^\ Finally, we define an interaction that is a combination 
of the instantaneous interaction "above" the cutoff and the divergent part of the 
exchange interaction: 

{g[g',\ M\A) {gmU+EX = {9W2\M\A) \g,g2)t^ + {gW,\M\A) \gig2)Ex ■ (5-86) 



Then 



+ 



{g[g'^\M^{A) \g1g2) = {gW2\M^{A) \gig2)KE + {9W2\ M^A) |(?i(?2)se 
{gW,\M\A) \gig2)coN + {9W2\ M\A) \g,g,)l^ + {gW,\M\A) 
+ {g[g',\M'{A) \g^g2)f^ + {g[g',\M'{A) {g^g^)^^ . (5.87) 

Perry showed that with a suitable definition of long-range interactions, a renormal- 
ization method that is similar to ours yields a logarithmically confining potential for 
quark- ant iquark bound states at O^g"^) [0. His calculation uses sharp step-function 
cutoffs and is based on an analysis of the part of the two-body interaction that is 
most singular in the limit in which the exchanged gluon has infinitesimal longitu- 
dinal momentum. The corresponding part of our interaction, which is contained in 
{gW2\M\A)\ 

5'ifi'2)iN+EX {9i92\-^^i^) l5'i5'2)iN) is similar to what Perry found 
in the quark-antiquark case. However, to determine whether or not our interaction 
is truly confining, we would have to do a careful analysis of the complete two-body 
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potential, not just the most singular part. This analysis would be complicated by the 
smooth cutoff that we employ, and we leave it for future consideration. 

5.3 Calculation of the Matrix Elements of the Invariant- 
Mass Operator in the Basis for Physical States 

The matrix elements {q',l',t',j 1^4"^ (A) \q,l,t,j), which appear in the eigenvalue 
equation, can be divided into contributions corresponding to the different terms in 
Eq. (^M) : 



(g', /', t',j\M\A) |g, /, t,j) = (g', t\j\M'{A) \q, I, t,j)^^ 
+ {q', I', t',j\M\A) \q, I, t,j)l^ + {q', /', t',j\M\A) \q, I, t, j)con 
+ (g', /', t', j\M\A)\q,l,t,j)l^+ {q', /' , t' , j\ M\A) \q,l,t, 
+ (g', r, t',j\M'{A) |g, /, t,j)f^ + (g', t',j\M'{A) |g, /, t,j)^^ . (5.88) 

In this section we express these different contributions in terms of integrals that can 
be computed numerically. These integrals fall into two classes: two-dimensional and 
five-dimensional. We treat each class of integral separately. 

Each of the terms in Eq. ( ^.88 ) is proportional to the plane- wave normalization 



factor 167r3p+5(3)(p _ p'). To make the remaining equations that we present sim- 
pler, we suppress this factor. We also take the limit e — > in any contribution to 
{q' ,1' ,t' , j\ A4'^{A) \q,l,t,j) that is finite in this limit. 

5.3.1 The Two-Dimensional Integrals 
The Kinetic Energy 



Using the definition of our basis in Eq. ( p. 491 ) and the expression for the free-state 



matrix element of the kinetic energy in Eq. ( ^.821 ), we find that 
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{q',l',t',j\M\A) \q,l,t,j) 



KE 



= ^...'^ dxL^\x)Lf\x) 

POO 

X / drr^Tt'{r)Tt{r). 
Jo 



(5.89) 



Note that the kinetic energy is infinite unless e > 0, whereas normahzabihty requires 
only e > -1/2. 

These integrals can be computed with standard numerical integration routines, 
but the results can have large errors for large values of the function indices /, /', t, 
and t', due to the oscillating nature of the basis functions. It is better to rewrite the 



integrals as sums that can be computed numerically with Mathematica ||32[ to any 
desired precision. Using the definitions of the basis functions, 



J',t',j\M\A)\q,l,t,j)^^ = 6, 



r(2e) 



X 



s + s'\ 



V A^'^ V A^'^) T{2e + m + m') 
2^ \m 2^ '^'''™'r(4e + m + m') 



m=0 m'=0 



(5.90) 



.s=0 s'=0 V ^ / 

The Finite Part of the Self-Energy Interaction 



The self-energy interaction conserves each particle's momentum, thus M| = Mp 
for this contribution. This means that the Gaussian cutoff factor in Eq. (|5.82|) has 
no effect on the self-energy. For the finite part of the self-energy, we use the same 
method for evaluating integrals that we use for the kinetic energy. This yields 



{q'J',t\j\M\A)\qJ,t,j) 

ir 

12. 

Nc9' nr 



SE- V5t,*'^^/oA' L dxLif\x)L\^\x)x{l - x) 







logx 



11 



(5.91) 

r(l + 2e + m + m') 



rra=0 m'=0 

X [?/'(l + 2e + m + m') - ^(2 + 4e + m + m')] 



(5.92) 
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where the digamma function ip{z) is given by 

dr(z) 



^(z) = (5.93) 



5.3.2 The Five-Dimensional Integrals 
The Contact Interaction 



Using the definition of our basis in Eq. ( ^.49|) and the expression for the free-state 



matrix element of the contact interaction in Eqs. (|5.61| ) and ( |5.82|) , we find that 



N a- 

{q'J',t',j\M\A)\q,l,t,j)^^^ = -^ 



dkdk'kk'e{k)e{k')Tf{k')Ti''\k) 



X j dxdx'e{x)e{l - x)e{x')e{l - x')Lf{x')Lf\x)e-^~'^'^'. (5.94) 

The reader may have noticed that this integral is not five- dimensional as we have 
implied, but rather four-dimensional. However, when we numerically compute the 
integrals that have more than two dimensions, it is most efficient if we combine them 
into one integral; so we want their integration variables and their ranges of integration 
to be identical. Thus we increase the number of dimensions of this integral by one by 
introducing an extra integral over 7 = — 0' using the identity 

l = ^/c?7^(7)^(27r-7). (5.95) 

This will help us to combine this integral with others that contain integrals over 7 
that cannot be done analytically. 

Since the integration domain of the exchange interaction is restricted so that 
X > x' when e — * [note the step function in Eq. (|5.8CI|) ], we would like to enforce 
this restriction in the other contributions. Using the identity 
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1 = e{x-x') + e{x' -x), 



(5.96) 



we break the longitudinal integral in Eq. ( |5.94| ) into two parts: 



J dxdx'e{x)e{i-x)e{x')e{i - x')Llf\x')Ll'\x)e-^''^'^i 
X [e{x - x') + e{x' - x)]. 



(5.97) 



In the second term, we let a; ^ 1 — a; and x' 1 — x', and then the longitudinal 
integral becomes 

dxdx'e{x)e{l - x)e{x')e{x - x')Lf{x')L'f\x)e-^''^'^' [l + (-1)'+''] . (5.98) 

Recall that we are restricting ourselves to the subspace of the states |g, Z, t, j) in which 
I + j is even if g 7^ 4, and I + j is odd if g = 4. Then since q = q' for the contact 
interaction, I and V must both be even or both be odd for this interaction. This 
means that the two terms in Eq. ( p. 98] ) are equal. Thus we can write the contact 
interaction contribution as follows: 



(g',/',t',j|A^2(A)|g,/,t,j) 



CON 



X Die 



(5.99) 



where 



D = dxdx'dkdk'djkk'e{x)e{l-x)e{x')e{x-x')e{k)e{k')e{j)e{2'K --f), (5.100) 



and 



i = LPix')L^f\x)T^f\k')T,^'\k) 



(5.101) 
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The Finite Part of the Exchange Interaction 

In order to simplify the contribution to {q' , I', t',j\ A^^(A) \q, I, t,j) from the finite 
part of the exchange interaction, we wish to change variables from to 7 = — 0': 

(fk^cPk'^ = dkdk'd(j)d(p'kk'e{k) e{k') e{(j)) 6(271 - 0) 0(0') e{27c - 00 
= dkdk'd'yd(f)'kk'e{k)e{k')e{-f + (f)')e{27c--f-(f)')e{(f)')e{27c-(f)'). (5.102) 

All the contributions to {q', I', t',j \ A1^(A) \q, I, t,j) depend on 7 only through depen- 
dence on cos 7 and sin 7. This means that we can use the identity 

/2'K — (j}' r27v 
(i7/(cos7, sin7) = / (i7/(cos7, sin7) (5.103) 
-<l>' Jo 

to write 

d^k^d^k'^ = dk dk' d^ d(f)' k k' e {k) e {k') e e {27C --f)e{<f)') 6(271 - 0'). (5.104) 

Inspection of Eq. (|5.78|) implies that the integrals in {q', I', t',j \ A^^(A) \q, I, t,j)^-^ 
depend on complex exponentials of 0' and 7. However, using the identity 

/ (i7/(cos7) sina7 = 0, (5.105) 
Jo 

where a is an integer, we can do the 0' integral in (g', i:', j| A^^(A) \q,l,t, j)^^ and 
write the remainder as a real quantity with integrals that depend on 7 only through 
cos 7 and sin 7. The result is 

{q',l',t',j\M\A)\qJ,t,j)l^ = -^ [ D^e-^-'^h^ (l^ + l^) 



X 



(1 



^/ + ^^^S + ^.n - x)x'(\ - x'^ 



(5.106) 



where some of the 51^1 's and S'^^l's are given by 
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S['l = cos(L7--2]7), 

S[^l = -cos{\j -l]-f)[x{l-x') + x'{l-x)], 

s[]l = ^cos{n)[x" + {1 - x'n 

Si'l = -^cos(L7--l]7)(^[l-^1+^'[l-^])(^'' + [l-^f), 

si]} = -^cos(j7)(l-2a;0, 

Sl'l = ^cos{[i-l]^){l-2x'){x[l-x']+x'[l-x]), 

S^l = ^cos{[j-2Mx' + {l-xn 

4S = cos(j7)[x'+(l-^r-2x'(l-x')], 

43 = -cos(j7)cos7(a;[l-x']+x'[l-x])(l-2x[l-x] -2x'[l-x']), 

'-'3,4 — 

45 = sm-fsm{n){2x^ -2x'^[l + x']+x'[l-2x']+x[l-2x' + 4:x'^]), 

Si]l = -Lcos(L7--2]7)(l-2x), 

'-'4,3 — 

Si]l = cos(j7)(l - 2a; - 2a;' + 4a;a;'), 

= -cos7COs(i7)(l-2a;)(l-2a;')(x + a;'-2a;a;'), (5.107) 



and 



^S='^?iL,,.r (5-108) 



The rest of the Si^l's and -^fi's are given by 

q,q q,q o j 



q{i) _ oil) _ 
*-'l,2 — ^2,1 — ^: 
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5« 



4,2 



"^4,1 



(5.109) 



where i = 1, 3. The 'S'^^^vs are given by 

:(2) _ c(l) 



(5.110) 



The Instantaneous and Exchange Interactions Combination 



Using similar methods for the contribution to (g', /', t' ,j\ A1^(A) |g, /, t,j) from the 
combination of the instantaneous interaction above the cutoff and the divergent part 
of the exchange interaction, we find that 

M 

87r3 

X (l-e~'^^ "i.K".K , (^x^x'){l-x^l-x')^ 



X (x{l - x)k''^ + x'{l — x')k^ — kk' [x'{l — x) + x{l — x')] cos 7 J 





1 


( 


X — x'y 


1 


- ' ] 


— 4 

FK 


1^1 k) 



(5.111) 



where 



Wq^q> = V,i cos(7[7' - 2]) + 5g,2V,2 cos(7[j + 2]) + 5^,3^,3 cos(7j) 

+ 5g,459',4Cos(7j). (5.112) 

The reader should note that the divergences from the two interactions that comprise 
(g', V, t', j\ Al^(A) \q, I, t, j)iN+Ex cancel, allowing us to take e — > in this contribution. 
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The Instantaneous Interaction Below the Cutoff 



The contribution to (g', t',j\ A^^(A) \q, I, t,j) from the instantaneous interaction 
below the cutoff is divergent. In this subsection, we extract the divergence, show 
that it cancels the divergent part of the self-energy, and compute the remainder of 
{q'J',t',j\M\A)\qJ,t,j)f^. 

After simplification, the complete contribution to {q', I', t',j \ 7W^(A) \q, I, t,j) from 
the instantaneous interaction below the cutoff is 



{q',l',t',j\M'{A)\q,l,t,j) 



B 
IN 



167r4 72. 



dx / dx' -F{x,x'), (5.113) 



X — X 



where 



F{x,x') 



X — X' 
X (1 — X +1 — x' 



-W,,g>e~^^ ^^^e-^^^^^'^^^'^ix + x') 



(5.114) 



To extract the divergence, we integrate Eq. ( |5.113|) by parts with respect to x': 



{q',l',t',j\M'{A) \q,l,t,j) 



IN 



I dx 



167r4 J2. 



\og{x - x')F{x,x')\ 



+ \og{x — x')F{x, x')\ ,_ + / dx'\og{x — x 



Bi + B2 + -B3. 



,.dF{x,x') 



dx' 



(5.115) 



The first contribution to {q' ,1' ,t' , j\ Ai^{A) \q,l,t,j)ff^ is 



5i 



levr^ e J2. 
X (1 — X + 1 — x') I , 



dx I d2A;xd2fc^eW,,,.e-^"'^^ie-2^"'^^^^^^(a; + a;') 



(5.116) 



To simplify this, we change variables from x to y: 



x = y{l- 3e) + 2e, 
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(5.117) 



and from k± and to Q± and N±: 



2 



k', = (5.118) 



Then 

X (2y[l-3e]+3e)(2-2|/[l-3e]-3e)|^,^^_^. (5.119) 



As e — s> 0, the only contribution that survives is 

Bi = ^,5,,,>\ogef^dyy{l-y)Lf{y)Lf\y)jSQ^Tif^{Q/2)T^^^ 

5,,,,8t,t'^^j\^' f^dxLf{x)Lf\x)x{l-x)\oge. (5.120) 



Since {g[g'2\ M'^ {K) |5'i5'2)se is the part of {g[g'2\ M"^ {K) |5'i5'2)se with the loge [see Eq. 
( ^.74| )], from inspection of Eq. (|5.91| ), we see that (g', t\j \ M.'^{A) |g, /, j)gE is just 



(g', j| A<^(A) |g, l,t,j)g^ with the [logx-11/12] factor replaced with - loge. This 
means that 

B, = - {q\r,t',j\M'iA) |g,/,t,j)sE- (5.121) 

Thus the divergence in the instantaneous interaction below the cutoff cancels the 
divergence in the self-energy. 

The second contribution to {q' ,1' ,t' , j\ Ai'^{A) \q,l,t,j)fj^ is 

B2 = f^" dxiog(x-x') f dkdk'd^kk'e(k)e(k')eh)e(27c --f)^- 

StT-^ J2e J X — X' 

X Wg^q,e-^'''^'^ie-^^''^^^^'^{x + x'){l-x + l-x')\^,^^. (5.122) 
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To evaluate this, we change variables from k' to s — k'/ ^Je. Then the leading term 
as e — > is 



^2 = -e"+5i^(i(jt.,oA!J /^'da;(2-x)loga; ldkdsdjks9{k)9{s)9{j)9{2Ti--f) 
X Ll^\x)Ti''\k)W,,,,e 





A-4 



A A. Ill 

x^(l-x)^ 



0, 



(5.123) 



since e > —1/2. 

To simplify the third contribution to {q',l',t',j\Ai'^{A) \q,l,t,j)f^, we take the 
derivative of F{x, x') and take e — > 0. Then 



87r3 



/ i^log(x - x')W,,,.Tif^{k')Lt\x)Ti''\k)i 



X n 

i=l m=l; m^i 



where 



(5.124) 



and 



E2 
E, 



E[ = 



X — X 

1, 



L]f>(x'), 
X + x' , 

1-x + l-x', 



(X — X 



>\2- 



Li^\x'), 



(5.125) 
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E'^ = -1, (5.126) 



{E[ = dEi/dx', except for i = 2) and 
dUf^ix') 



LfHx') 



X 



dx' 
2 

^/\(e) _/m,'-l 



e - [x'{l - x')Y-^-'[l - 2x'] i: \%x'-^' + {x\\ - x')X-'^ 



m,'=0 

i: m'AP^,x'-'-! (5.127) 



m'=l 



and 



A', 



dx' (\ — x'Y' [x — x'Y' ' 
dx' (x — x' 



^ ^ = ^ - ■ (5.128) 



This means that we can write the contribution to (g', t', j| A^2(A) |g, Z, t, j) from the 
instantaneous interaction below the cutoff as a divergent part and a finite part: 



(g',/',t',j|A<2^A)|g,/,t,j)fN = - (g', t', j| ^^^(A) |g, /, t, j) 



SE 



+ (g',/',t',j|A^2(A)|g,/,t,j)f/, (5.129) 

where 

(g',/',t',j|A^2(A)|g,/,t,jC = -^\D\og{x~x')W,^,,Tf{}^)Lf{x)Tt\k) 

5 5 

X e-^"*(^^^+^?^)^E: n E^. (5.130) 

2=1 m=l; m^j 



Using the results of this section, the expression in Eq. (|5.88|) for the matrix 
elements of the IMO becomes 
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I', t', j I M\A) \q, I, t, j) = {q', I', t', j I M\A) \q,l,t, j)^^ 
+ {q', I', t',j\ M\A) \q, I, t,j)l^ + {q', I', t',j\M\A) \q, I, j)coN 
+ {q',l',t',j\M\A)\q,l, t, {q', I', t' , j \ M'{A) \q,l,t, j)in+ex 

+ {q',l',t',j\M\A)\q,l,t,j)f^. (5.131) 

Each of these terms is finite and we have taken e ^ everywhere. We have written 
the first two terms as sums that can be computed numerically, and the four remaining 
terms as five-dimensional integrals that can be grouped into one integral suitable for 
numerical calculation. (See Appendix C for a discussion of some of the technical 
issues involved in the numerical calculation of these matrix elements.) Once we have 
computed these matrix elements, we can diagonalize the matrix to obtain glueball 
states and masses. 

5.4 Results and Error Analysis 

In this section we diagonalize the IMO matrix, obtaining glueball states and 
masses. We then discuss the sources of error in the calculation. We begin by dis- 
cussing the procedure that we use to calculate the results. 

5.4.1 The Procedure for Calculating Results 

We represent a state using the notation J^^*^, where J is our best guess for the 
spin of the state, P is our best guess for the parity of the state, C represents the 
charge-conjugation eigenvalue of the state (it is always -|- because we have an even 
number of gluons), and j is the projection of the state's spin onto the 3-axis. We 
need to distinguish states with identical J's and P's and different j's because we do 
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not have manifest rotational symmetry. If J = 0, we omit the subscript j in the 
state notation. We use an asterisk in the state notation next to the value of C to 
denote an excited state with the given quantum numbers. We will base our guesses 
for J and P on the numerical degeneracies of the states that have identical J's and 
P's and different j's, and on the ordering of the states according to lattice data (see 
the discussion below). We consider only the five lightest glueballs (not counting as 
distinct those states that differ only in their value of j). The five lightest glueballs 
have spins J < 2. This means that we need to consider only \j\ < 2. For a given |j|, 
the states with j = \j\ and j = are degenerate and simply related (see Appendix 
C); so we explicitly consider only j = 0, 1, 2. We consider nine values of the coupling: 
•^A = S'a/I^^) ~ 0.1, 0.2, 0.3, 0.9. To calculate our results, we implement the 
following four-step procedure. 

We execute the first step for all pairs {j,a^). In this step, we define A = 1 and 
diagonalize {q',l',t',j\A4'^{A) \q,l,t,j) with all 4 spin basis functions {q — 1,2,3,4), 
but with only the lowest transverse-magnitude basis function {t — 0) and the two 
lowest longitudinal basis functions (/ = 0, 1). (It is necessary to use an even number of 
longitudinal basis functions so that symmetry and antisymmetry of the wave function 
under x ^ 1 — x are equally represented.) We perform this diagonalization as a 
function of the basis-function parameters d and e that determine the widths of the 
transverse-magnitude and longitudinal wave functions, respectively, and we find the 
values of d and e that minimize the ground-state mass. This yields what we consider 
to be the optimal wave-function widths for each pair {j, ) . 

We also execute the second step of the procedure for all pairs {j,a^). In this 
step, we fix d and e to be their computed optimal values and again define A = 1. 
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We diagonalize the matrix with all 4 spin basis functions, transverse-magnitude 
basis functions, and Ni = 2Nt longitudinal basis functions, for a total of SA''^^ basis 
functions, with A^t = 1,2,3,..., 10. We use twice as many longitudinal functions as 
transverse-magnitude functions because \q, I, t,j) is zero if / + j is even and g = 4, or 
a 1+ j is odd and q ^ 4:. We want to use as many basis functions as possible, but we 
find for all pairs (j, a^) that when A^t > 7, the statistical errors from the Monte Carlo 
integrations of the matrix elements become overwhelming and the spectrum and wave 
functions become unreliable. The evidence of the breakdown is sudden contamination 
of the low-lying wave functions with high-order components. (See Appendix C for a 
more complete discussion of this topic.) Thus in the remainder of our procedure, we 
analyze the results that we find in this step when A^t = 7, which corresponds to 392 
basis functions. 

In the third step of the procedure, we use the mass of the state (our most 
numerically reliable state) to determine the value of A for each a^. To do this, we 
note that A is the only mass scale in the problem. This means that the mass of the 
0~+ state, Mq-+, can be written 



where 6 is a dimensionless function of . Since we defined A = 1 in the second step 
of our procedure, the diagonalization of the IMO as a function of yielded &(q;^). 
In this third step, we consider A to be a parameter and define Mq-+ to be a constant. 
Then for a given coupling, we can use the results of the second step of our procedure 
to write the cutoff in units of Mo-+: 



Mo-+ = b{aJA, 



(5.132) 



A/Mo-+ 



1 



(5.133) 
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Figure 5.5 shows the result for the third step of our procedure: a plot of the 
coupling as a function of the cutoff. When > 0.7 it is not a single- valued function 
of the cutoff. This is an indication that the coupling is too large. For this reason, 
we consider only < 0.7 in the remainder of our procedure. When < 0.7, 
the coupling decreases as the cutoff increases, as expected. However, depends 
on A more strongly than one may expect. We expect that perturbative pure-glue 
QCD would indicate that ~ 1/lnA, but the result that we get is much closer to 

~ exp(— aA), where a is a constant. The reason for this is that the truncation of 
the perturbative series for A^^(A) and the truncation of the free-sector expansion of 
the states introduce spurious cutoff dependence in our results for physical quantities. 
Mq-+ is one such physical quantity. The spurious cutoff dependence of Mq-+ is 
manifested through incorrect dependence of 6(0^) on a^. This means that A has to 
compensate by depending on incorrectly in order to keep Mq-+ a constant function 
of a^. This results in the strong dependence of on A that is shown in Figure 5.5. 

Using the recent anisotropic Euclidean lattice results of Morningstar and Peardon 
||33| , we can make a rough estimate of the range over which our cutoff is varying in 
Figure 5.5. They found that the mass of the 0"+ state is Mo-+ = 2.590 ±0.040 ±0.130 
GeV. This means that our cutoff is varying from about 3.1 - 6.0 GeV in Figure 5.5. 

The fourth step of our procedure is to determine the optimal value of the cutoff, 
or equivalently, the optimal value of the couphng. We use two criteria to determine 
this. First we determine the value of the cutoff for which the computed masses are 
most independent of the cutoff. Figures 5.6, 5.7, and 5.8 show the masses of our states 
with j = 0, 1,2, respectively, as functions of the cutoff. The masses and the cutoff 
are displayed in units of Mq-+. (Recall that Mq~+ was defined to be independent of 
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the cutoff in the process of defining the cutoff as a function of the couphng.) The 
seven values of the cutoff at the points that we display correspond, from right to left, 
to = 0.1, 0.2, 0.3, . . . , 0.7. It is difficult to tell from these plots where the cutoff 
dependence is weakest (more points are needed), but we see that the dependence is 
relatively weak from = 0.5 to = 0.7, which corresponds to A/Mq-+ = 1.33 to 
A/Mo-+ = 1.20. 

The second criterion that we use to determine the optimal cutoff is the degree 
to which the states with a given J and P and different j's are degenerate. This 
determines the cutoff that minimizes the violation of rotational symmetry. We find 
that these degeneracies are best when = 0.5. Given this and the fact that Figures 
5.6 - 5.8 indicate that the cutoff dependence of the masses is weak when = 0.5 to 

= 0.7, we determine that = 0.5 is the optimal coupling, and thus the optimal 
cutoff is A/Mo-+ = 1.33. Using the result of Morningstar and Peardon for the mass 
of the O""*" state, we estimate that this cutoff is about 3.4 GeV. 

5.4.2 Results 

Now we present our main results. Our glueball masses for = 0.5 are summa- 
rized in Table 5.1, in units of the mass of the ground state (the 0"''+ state). Table 
5.1 also shows an average of lattice results from a number of different calculations for 



the sake of comparison [^. The uncertainties in our results that we report in Table 
5.1 are only the statistical uncertainties associated with the Monte Carlo evaluation 
of the matrix elements of A^^(A). The full uncertainties are much larger (see the 
discussion of sources of error below). We list three values of the masses for the 2"*"+ 
and 2++* states for our calculation, corresponding to j = 0, 1, 2. In each case the 
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three masses would be degenerate if our calculation were exact. Our results agree 
with the lattice results quite well, perhaps better than we should have expected. 

We display our spectrum graphically in Figure 5.9. The masses are plotted in 
units of the mass of the 0"''+ state and the vertical widths of the levels represent the 
statistical uncertainties in the masses. The black lines connect the states that we 
believe should be degenerate. We see that the and 2q glueballs are relatively 
degenerate with their 2^'"^ and 2^^* counterparts, and the 22"^* is not too bad, but the 
glueball is much too hght. Our labehng of the states and subsequent assignment 
of the expected degeneracies are based on the ordering of the lattice states and the 
apparent degeneracies of the 20""*" and 2^^* states with the 2^"*" and 2^"^* states, 
respectively. 

We want to show some of the features of the glueball wave functions. Rather 
than presenting the spin-dependent wave functions themselves, we present more illu- 
minating spin-independent probability densities. A glueball state has the plane-wave 
normalization shown in Eq. ( |5.51| ) as long as the wave function ^^"^^(x, kj_) satisfies 



(fk^dxe{x)e{l -x)Y, M'ls^ix, kA.)\^ = 1. (5.134) 



This implies that 



poo I 

Jo ^[a)J ^^n(x,A;/A) = 1, (5.135) 
where we define the dimensionless probability density Il{x,k/A) by 

U{x,k/A) = 27TAkJ2\<!>Z,{x,k^)f. (5.136) 



S1S2 



We show the probability densities for our glueballs in Figures 5.10 - 5.18. The 
masses of the states tend to increase as the probability densities move away from the 

126 



region x ~ 1/2 and more towards the edges. This is what we expect based on the 
form of the kinetic energy of a free state. Notice that the probabihty density for the 
22"*" glueball is peaked around the region x ~ 1/2 and looks similar to the probability 
density for the 0++ glueball. This is consistent with its small mass. On the other 
hand, the probability density for the 22~^* glueball is away from the region x ~ 1/2, 
like the probability densities for its j = and j — 1 counterparts, and its mass is 
close to theirs. 

5.4.3 Error Analysis 

We now turn to a discussion of the sources of error in our calculation. The sources 
of error are: 

• truncation of the renormalized IMO at 0{g^); 

• truncation of the free-sector expansion of physical states at two gluons; 

• truncation of the basis-function expansion of wave functions; 

• numerical approximation of the matrix elements of A4^{A). 

We do not know how to estimate the size of any physical effects that require 
nonperturbative renormalization. However, we can naively estimate the size of the 
effects of higher-order perturbative renormalization. We have calculated the matrix 
elements of A1^(A) through C(a^); so the corrections to these matrix elements should 
be 0(q;^). This translates to corrections to the mass spectrum of 0{a'^/2). Since we 
have used = 0.5, we estimate that the uncertainty in the mass spectrum from the 
effects of higher-order perturbative renormalization is about 13%. 

We do not know how to estimate the size of any physical effects that require an 
infinite number of particles. In fact, until we include at least two free sectors, it is 
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impossible to directly estimate the size of corrections from higher free sectors. How- 
ever, we can use the lack of degeneracy of the state with the 20"*" and states 
to estimate these corrections. According to Table 5.1, the discrepancy in the various 
2++ states is about 33%, if we believe the quoted lattice result. Since the uncertainty 
in the mass spectrum from effects of higher-order perturbative renormalization is 
around 13%, an uncertainty of 30% due to the truncation of the free-sector expansion 
is necessary to explain the lack of rotational symmetry in the spectrum (neglecting 
the other sources of error, which we expect to be small). 

As we mentioned, when we increase the number of basis functions that we use 
to represent the wave functions, we find that we reach a point where the statistical 
errors from the Monte Carlo evaluations of the matrix elements are overwhelming and 
cause our results to become completely unreliable (see Appendix C). For this reason, 
we have to truncate our basis-function expansion for the wave functions at A^t = 7 
transverse-magnitude functions. (A^i = 2A"t = 14, and there are 4 spin basis functions, 
for a total of 392 basis functions.) This truncation results in additional errors in our 
results. In Figures 5.19, 5.20, and 5.21, we show the convergence of the masses of the 
states with j = 0, 1, 2 respectively, as functions of At- The masses do not decrease as 
rapidly as functions of the number of states as one might expect. This is primarily 
because we have already optimized the states quite a bit by determining the widths 
of the transverse-magnitude and longitudinal basis functions that minimize the mass 
of the 0"*"+ state (using At = 1)- Our best guess for the uncertainty that we introduce 
into the spectrum when we truncate the basis-function expansion, based on Figures 

5.19 - 5.21, is a few percent^. 

-•^^ Technically, this is not an uncertainty because improving the states can only reduce their masses, 
according to the variational principle. However, our discussion of errors is not meant to be rigorous. 

128 



We can estimate the uncertainty in our results associated with the Monte Carlo 
evaluation of our matrix elements. To do this, we compute our results with = 0.5 
four times, obtaining statistically independent results, and we compute the standard 
deviations of the masses that we obtain. This leads us to estimate that the uncertainty 
in the spectrum from the Monte Carlo routine is 1 - 2%. This is the uncertainty that 
we report in Table 5.1 and Figure 5.9. Combining this uncertainty with the others 
leads us to estimate that the total uncertainty in our results is about 33%. 
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Figure 5.5: The coupling as a function of the cutoff. We show the cutoff in units of 
the mass of the state, and we use 14 longitudinal basis functions, 7 transverse- 
magnitude basis functions, and 4 spin basis functions, for a total of 392 basis functions. 
Using the recent anisotropic Euclidean lattice result of Morningstar and Peardon for 
the mass of the ^ state [Q, we estimate that the cutoff is roughly varying from 
about 3.1 - 6.0 GeV in this figure. 



130 



M/Mo-+ 1.3 



0++* 


1 

▲ ▲ 


▲ 


▲ 


1 


1 


1 


+ + * 


o o 


o 


o 


▲ 


A 




2++ 


• • 


• 


• 


o 
• 


o 

• 


▲ 
• 


- 0-+ 


o o 


o 


o 


o 


o 


o - 












■ 


■ - 






■ 


■ 


■ 






~ 0++ 


■ ■ 

1 






1 


1 


1 



1.00 1.27 1.54 1.81 2.08 2.35 



A/Mo-+ 



Figure 5.6: The masses of the five Ughtest glueballs with j = 0, as functions of the 
cutoff. The masses and the cutoff are displayed in units of the mass of the 0^"*" 
state. The seven values of the cutoff at the points that we display correspond, from 
right to left, to = 0.1, 0.2, 0.3, . . . , 0.7. We use 14 longitudinal basis functions, 7 
transverse-magnitude basis functions, and 4 spin basis functions, for a total of 392 
basis functions. 
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Figure 5.7: The masses of the two Ughtest glueballs with j = 1, as functions of the 
cutoff. The masses and the cutoff are displayed in units of the mass of the 0^"*" 
state. The seven values of the cutoff at the points that we display correspond, from 
right to left, to = 0.1, 0.2, 0.3, . . . , 0.7. We use 14 longitudinal basis functions, 7 
transverse-magnitude basis functions, and 4 spin basis functions, for a total of 392 
basis functions. 
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Figure 5.8: The masses of the two Ughtest glueballs with j = 2, as functions of the 
cutoff. The masses and the cutoff are displayed in units of the mass of the 0^"*" 
state. The seven values of the cutoff at the points that we display correspond, from 
right to left, to = 0.1, 0.2, 0.3, . . . , 0.7. We use 14 longitudinal basis functions, 7 
transverse-magnitude basis functions, and 4 spin basis functions, for a total of 392 
basis functions. 
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State 


M/Mq++ 


Lattice 34| 


0~+ 


1.38 ±0.02 


1.34 ±0.18 


2++ 


1.58 ±0.01 
1.58 ±0.02 
1.11 ±0.01 


1.42 ±0.06 


2++* 


1.70 ±0.01 
1.68 ±0.02 
1.62 ±0.02 


1.85 ±0.20 


0++* 


1.77 ±0.02 


1.78 ±0.12 



Table 5.1: The glueball masses from our calculation compared to an average of lattice 
results from a number of different calculations |Q. We display the masses in units 
of the mass of the 0"*"^ state. The uncertainties for our results are only the statistical 
uncertainties associated with the Monte Carlo evaluation of the matrix elements of 
A1^(A). The three values of the masses for the 2"'""'" and 2"'""^* states for our calculation 
correspond to j = 0, 1, 2. We use the optimal coupling = 0.5, with 14 longitudinal 
basis functions, 7 transverse-magnitude basis functions, and 4 spin basis functions, 
for a total of 392 basis functions. 
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Figure 5.9: Our glueball spectrum. The masses are plotted in units of the mass of the 
0^^ state and the vertical widths of the levels represent the statistical uncertainties 
in the masses. The black lines connect the states that we believe should be degen- 
erate. We use the optimal coupling = 0.5, with 14 longitudinal basis functions, 
7 transverse-magnitude basis functions, and 4 spin basis functions, for a total of 392 
basis functions. 
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Figure 5.10: The probability density of the O"*""*" glueball. We use the optimal couphng 
CK^ = 0.5, with 14 longitudinal basis functions, 7 transverse-magnitude basis functions, 
and 4 spin basis functions, for a total of 392 basis functions. 
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Figure 5.11: The probability density of the 0"+ glueball. We use the optimal couphng 
CK^ = 0.5, with 14 longitudinal basis functions, 7 transverse-magnitude basis functions, 
and 4 spin basis functions, for a total of 392 basis functions. 
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Figure 5.12: The probability density of the 2q glueball. We use the optimal couphng 
CK^ = 0.5, with 14 longitudinal basis functions, 7 transverse-magnitude basis functions, 
and 4 spin basis functions, for a total of 392 basis functions. 
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Figure 5.13: The probability density of the 2^ glueball. We use the optimal couphng 
CK^ = 0.5, with 14 longitudinal basis functions, 7 transverse-magnitude basis functions, 
and 4 spin basis functions, for a total of 392 basis functions. 



139 




Figure 5.14: The probability density of the 22 glueball. We use the optimal couphng 
CK^ = 0.5, with 14 longitudinal basis functions, 7 transverse-magnitude basis functions, 
and 4 spin basis functions, for a total of 392 basis functions. 
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Figure 5.15: The probability density of the 2q^* glueball. We use the optimal couphng 
— 0.5, with 14 longitudinal basis functions, 7 transverse-magnitude basis functions, 
and 4 spin basis functions, for a total of 392 basis functions. 
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Figure 5.16: The probability density of the 2i^* glueball. We use the optimal couphng 
— 0.5, with 14 longitudinal basis functions, 7 transverse-magnitude basis functions, 
and 4 spin basis functions, for a total of 392 basis functions. 
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Figure 5.17: The probability density of the 22 glueball. We use the optimal couphng 
= 0.5, with 14 longitudinal basis functions, 7 transverse-magnitude basis functions, 
and 4 spin basis functions, for a total of 392 basis functions. 
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Figure 5.18: The probability density of the 0++* glueball. We use the optimal couphng 
= 0.5, with 14 longitudinal basis functions, 7 transverse-magnitude basis functions, 
and 4 spin basis functions, for a total of 392 basis functions. 
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Figure 5.19: The masses of the five Ughtest glueballs with j = 0, in units of the cutoff, 
as functions of the number of transverse-magnitude basis functions, Nt. We use the 
optimal couphng = 0.5, with 4 spin basis functions and A^i = 2Nt longitudinal 
basis functions, for a total of SN^ basis functions. 
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Figure 5.20: The masses of the two Ughtest glueballs with j = 1, in units of the cutoff, 
as functions of the number of transverse-magnitude basis functions, Nt. We use the 
optimal couphng = 0.5, with 4 spin basis functions and A^i = 2Nt longitudinal 
basis functions, for a total of SN^ basis functions. 
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Figure 5.21: The masses of the two Ughtest glueballs with j = 2, in units of the cutoff, 
as functions of the number of transverse-magnitude basis functions, Nt. We use the 
optimal couphng = 0.5, with 4 spin basis functions and A^i = 2Nt longitudinal 
basis functions, for a total of SN^ basis functions. 
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CHAPTER 6 



CONCLUSION 



We have presented a formalism that allows the physical states of a quantum field 
theory to rapidly converge in a free-sector expansion. In this approach, we force 
the free-state matrix elements of the IMO to satisfy three conditions to make the 
desired expansion possible. First, the diagonal matrix elements of the IMO must 
be dominated by the free part of the IMO. Second, the off-diagonal matrix elements 
of the IMO must quickly decrease as the difference of the free masses of the states 
increases. Third, the free mass of a free state must quickly increase as the number of 
particles in the state increases. 

We assume that we can use perturbation theory to derive the operators of the 
theory, and if this is valid, then the first condition is automatically satisfied. To satisfy 
the second condition, we place a smooth cutoff on the IMO. We use LFFT so that the 
effects of the vacuum are isolated in particles with zero longitudinal momentum, and 
we remove these particles from the theory with the intent of replacing their physical 
effects with interactions. This makes it reasonable to expect that the third condition 
on the IMO is satisfied automatically due to the free-particle dispersion relation of 
LFFT. 
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The cutoff tfiat we use violates a number of physical principles of LFFT. How- 
ever, by requiring the IMO to produce cutoff-independent physical quantities and 
by requiring it to respect the unviolated physical principles of the theory, we are 
able to derive recursion relations that uniquely determine the IMO to all orders in 
perturbation theory. 

We have apphed our method to two field theories. First we worked with massless 
0^ theory in six dimensions. We derived the recursion relations that determine the 
IMO and used them to calculate some second-order matrix elements of the IMO. We 
showed that these matrix elements naturally obey the restrictions that we have placed 
on the IMO and that they allow our IMO to yield cutoff-independent perturbative 
scattering amplitudes. We then derived part of the IMO's emission matrix element 
to third order in perturbation theory in terms of five-dimensional integrals, and com- 
puted it numerically using Monte Carlo methods. We used this matrix element to 
derive the scale dependence of the coupling and showed how it relates to the standard 
result that one finds with dimensional regulation. 

The second field theory to which we applied our method is pure-glue QCD. We 
derived the recursion relations for the IMO in this theory and then turned to a cal- 
culation of the physical spectrum. For this calculation, we approximated all physical 
states as two-gluon states. We calculated the color parts of the states analytically, 
and we expanded the states' momentum and spin degrees of freedom in terms of ba- 
sis functions. We designed the states to be simultaneous eigenstates of the IMO, the 
three-momentum operator, and the projection on the 3-axis of the internal rotation 
generator. 
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Using our recursion relations for the IMO, we calculated to second order in per- 
turbation theory the two-gluon to two-gluon matrix element of the IMO, which is 
required for the calculation of the spectrum. We then used it to calculate the IMO 
matrix in terms of the basis functions. We showed that the infrared divergences in 
the matrix from exchanged gluons with infinitesimal longitudinal momentum cancel 
when treated properly. 

In order to diagonalize the IMO matrix, we computed the five-dimensional inte- 
grals in the matrix elements using Monte Carlo methods. We calculated the glueball 
spectrum for a range of couplings and found that we could not use more than about 
400 basis functions without the statistical errors becoming overwhelming. We used 
the mass of our 0~^ glueball to compute the nonperturbative cutoff dependence of 
the coupling, and we analyzed the cutoff dependence of the spectrum. We found that 
the cutoff that minimizes our errors is A/Mo-+ = 1.33. The corresponding coupling 
is CK^ = 0.5. We presented the probability densities for our glueballs and found that 
our results for the spectrum compare favorably with recent lattice data. The largest 
discrepancy seems to be the 22~^ state, which is much too light. Finally, we analyzed 
the errors in our calculation from the various possible sources, and estimated the total 
uncertainty in our spectrum to be 33%. 

There are two main paths that we can take for future work with our approach. 
The first path is to further test our method with the theories that we have considered 
so far. Since the scalar theory is relatively simple, it would be interesting to use it 
compute the IMO to higher orders in perturbation theory. This would require us to 
solve the integral equation for the cutoff-independent reduced interaction and could 
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be used to further check our conjecture that our IMO leads to correct scattering 
amphtudes order-by-order in perturbation theory. 

In pure-glue QCD, we can further test our approach for calculating physical states 
by computing the IMO to higher orders in perturbation theory and by keeping more 
free sectors in the expansion of the states. However, to keep more free sectors in the 
expansion of the states, we have to be able to calculate IMO matrices that have more 
degrees of freedom, while controlling the statistical errors in the spectrum. This means 
that we need a better algorithm for determining how accurately individual matrix 
elements of the IMO have to be computed in order to get a desired uncertainty in the 
spectrum. We could also use a better basis that requires fewer momentum functions 
to represent a wave function. Overcoming these problems will be challenging, but it 
is important to test our method by studying the rate of convergence of the free-sector 
expansion of states as a function of both the cutoff and the masses of the states. 

Another test of our method that we can do with pure-glue QCD is to analyze 
the interaction that we have derived in this dissertation to test the conjecture that 
it is logarithmically confining. It would also be interesting to analyze the long-range 
parts of higher-order interactions to see if the perturbative series for the interaction 
is building towards a linearly confining potential. Analyzing the long-range parts of 
higher-order interactions may be much easier than computing the interactions in their 
entirety. 

The second main path that we can take in the future is to extend our method to 
other theories and operators. In order to compute quantities that can be compared 
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with experiment, we wish to extend our method to full QCDQ. This is complicated 
for two reasons. First, there is additional algebraic and numerical complexity from 
the vertices involving quarks. Second, quark masses complicate the method for deter- 
mining the IMO because they increase the number of reduced interactions that can 
be cutoff-independent 0| . In addition, if large and small quark masses are considered 
simultaneously, then efficient numerical representation of the states and accurate cal- 
culation of the IMO's matrix elements become more difficult. Masses also quickly 
enlarge the parameter space that must be explored to compare to experimental data. 

We can also extend our method by applying it to the computation of operators 
other than the IMO, such as the rotation generators, the parity operator, and currents. 
The rotation generators and the parity operator are of particular interest because they 
may aid in the classification of the physical states of a theory. 

In summary, there are many avenues of research that must be explored, and some 
of them are quite complex. However, all the improvements that we have discussed 
are necessary if we are to accurately represent the physical states of quantum field 
theories as rapidly convergent expansions in free sectors. 



'^^We are not thinking of including QCD effects that require nonperturbative renormahzation or 
an infinite number of particles because these would require a new method rather than an extension 
of our current approach. 
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APPENDIX A 



PARTICLE-NUMBER DEPENDENCE OF FREE MASSES 



The purpose of this appendix is to give a quahtative argument that shows that 
in LFFT, the free mass of a free state should rapidly increase with the number of 
particles in the state. For simplicity, we consider the case of massless particles and 
limit ourselves to analyzing free states that are components of physical states with a 
single physical particle (e.g. single-glueball states in pure-glue QCD). If some of the 
particles in a free state are massive, then the free mass of the state should increase even 
more rapidly with the number of particles in the state. Extension of our argument to 
free states that are components of physical states with multiple physical particles is 
straightforward. 

The free IMO satisfies the Fock-space eigenvalue equation 

MIJF)^M'^\F), (A.l) 

where |F) is any free state. We define Np to be the number of particles in |F). The 
eigenvalue Mp depends on Np: 

Ml^Y.—^ (A.2) 

i=l -^i- 
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where Xi is the fraction of the total longitudinal momentum of \F) that is carried 
by particle i, and fi± is the particle's transverse momentum in |-F)'s center-of-mass 
frameQ. 

Let us assume that rlj_ is approximately independent of i. Then 

Np ^ 

M^~(rl)typicalE-- (A-3) 

i=l ^« 

This is minimized if the particles share the longitudinal momentum equally, in which 
case Xi = l/Np. Then 

M^F ~ iV|(rl)typical. (A.4) 

This shows that the free mass of a free state should increase at least quadratically 
with the number of particles in the state. If r^^ is approximately independent of i 
and the particles do not share the longitudinal momentum equally (this is the most 
likely scenario; see the probability densities in Chapter 5, Section 4.2), then Mp will 
increase even faster with A'"^.[^ This rate of growth should be sufficient as long as 
the scale of the free mass, which is set by the proportionality constant {rj_) typical, is 
not negligible compared to the cutoff. 

In principle, we could choose the cutoff to be much larger than (rj)typicab but 
instead we choose it to be such that our errors are minimized. The dynamics of the 
system determine the value of this optimal cutoff, and we believe that it is unnatural 

for it to be much larger than (^l)typicai- In a confining theory like QCD, the reason 

^^To analyze a free state that is a component of a physical state with multiple physical particles, 
we would divide Mp into a separate contribution from each physical particle and analyze each 
contribution in terms of the associated physical particle's total longitudinal momentum and center- 
of-mass frame. 

^^Note that in equal-time field theory the free energy of a state increases only linearly with the 
number of particles, assuming they have non-negligible momenta. 
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for this is particularly clear. If the typical center-of-mass transverse momentum in a 
state is very small compared to the optimal cutoff, then the state will have a width in 
transverse position space that is much greater than the inverse of this cutoff (which 
is typically the size of the bound state) . In a confining theory, these states should be 
highly suppressed. Therefore, we do not expect (r^)typicai to be negligible compared 
to the cutoff. 

To conclude, in LFFT, the free mass of a free state should increase rapidly with 
the number of particles in the state. 
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APPENDIX B 



PROOF OF UNITARITY OF U{A,A') 



The purpose of this appendix is to prove that f/(A, A'), as defined by 

dU{A,A') 



and 



T(A)[/(A,A') (B.l^ 



f/(A,A) = l, (B.2) 

is unitary as long as T(A) is anti-Hermitian and hnear. Assume that T(A) is anti- 
Hermitian and hnear. Then a solution for U (A, A') is 

U(A, A') = 1 + / d(A"-^)T(A")U(A", A') 

= 1 + 1^^^ rf(A"-^)T(A") + J^^ ^ d{A"-^)T{A") J^^ ^ d{A"'-^)T{A"') + ■■■. (B.3) 

This solution shows that since T(A) is linear, f/(A,A') is linear^". 
Now multiply Eq. on the left by U\A, A'): 

f^^A.A')^^^^ = f/nA,A')T(A)t/(A,A') 

= - "f/^(A,A')r(A)f/(A,A')]^ 

^°We assume that the expansion of the integral equation converges. If it does not, then the 
perturbative solution for U{h.,K') in Chapter 3 is probably useless, in which case this proof is 
irrelevant. 
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'J^!^UiA AO 



(B.4) 



This is the same as 



d 



d{A- 



U\A,A')U{A,A') 



0, 



(B.5) 



which imphes that f/^(A, A')f/(A, A') is independent of A. Since Eq. (|B.2|) imphes 

that 



U\A',A')U{A',A') = 1, 



we conclude that 



U\A,A')U{A,A') = 1. 

Now multiply Eq. ( |BTTD on the right by W{A,A'): 
dU{A,A') 



d{A- 



-U\A,A') = T{A)U{A,A')U\A,A') 



and take the adjoint of this: 

,,dW{A,A') 



U{A,A' 



rf(A-4) 

Adding Eqs. (p.8| ) and (|B.9|) gives 
d 



-f/(A,A')f/^(A,A')T(A). 



d{A 



-4\ 



{U{A,A')U\A,A')) = [t{A),U{A,A')U\A,A'] 



From Eq. 



U{A',A')U\A',A') = 1. 



(B.6) 



(B.7) 



(B.8) 



(B.9) 



:b.io) 



:b.ii^ 
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f/(A, A')f/'I"(A, A') is a function of A that is uniquely determined by Eqs. (|B.10|) and 
( [B.ll| ). Therefore, since the statement 



U{A,A')U\A,A') = 1 



(B.12) 



is a solution to Eqs. (|B.10|) and (|B.11|) , it is a true statement. Since U (A, A') is linear 
and satisfies Eqs. ( [B.7|) and ( [B.12 ), it is unitary. 
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APPENDIX C 



TECHNICAL ISSUES IN THE NUMERICAL 
CALCULATION OF PURE-GLUE QCD MATRIX 

ELEMENTS 



In this appendix we discuss some of the technical issues involved in the numerical 
calculation of the pure-glue QCD matrix elements {q',l',t',j\Ai'^{A)\q,l,t,j). In 
Section C.l we discuss how we put the matrix elements into a form that is amenable 
to numerical calculation. In Section C.2 we briefly cover three topics: we show how 
the glueball state with J'^ eigenvalue —j can be written in terms of the glueball 
state with J'^ eigenvalue j, we list a few tricks that allow us to reduce the number 
of matrix elements that we must compute, and we present our method for estimating 
how numerical uncertainties in the matrix elements translate into uncertainties in the 
spectrum. 

C.l Preparation of Integrals for Monte Carlo 

There are two types of contributions to {q' ,1' ,t' , j\ J^'^{A) \q,l,t,j): finite sums 



and five-dimensional integrals. We use Mathematica to evaluate the finite sums 



to as many digits as we wish. To evaluate the integrals, we combine them into one 



integral and use the VEGAS Monte Carlo routine ||29|. It takes a bit of work to put 
the integral into a form that will converge. 
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There are two main difficulties with getting the integral to converge. The first 
difficulty is that (g', t',j\ Ai'^{A) \q, I, t,j)^-^ looks divergent: the sum in 

Eq. ( |5.13CI| ) diverges. This divergence is misleading because it actually contributes 



nothing to the integral (assuming that we calculate the integral carefully; see the 
discussion below). Left unchecked, this false divergence prevents the integral from 
converging with VEGAS. To rectify the problem, we want to subtract the false diver- 
gence from the integrand. Since it integrates to zero, this is allowed. 

The second main difficulty with getting the integral to converge is roundoff error. 
Even after we subtract the false divergence from the integrand, the integrand peaks 
around x' = x, and in this region there are large cancellations in some of the quantities 
that we have defined. To prevent these cancellations from causing roundoff error, 
we rewrite these quantities so that the cancellations are explicit, before we turn to 
numerics. 

C.1.1 Subtraction of the False Divergence 

We begin by defining a set of variables that is natural for dealing with the false 
divergence. We define 

Tj = X — x'. (CI) 

We change variables from and k'j_ to the dimensionless transverse variables f± and 
w±: 

dk\—k\ 



We define the angle between rj_ and w± to be /3: 
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fj_ • w± — rw cos p. (C.3) 

We also define 



r± — ^Jt"^ -\- rjw^ ± 2rw^ cos (C.4) 



and then we can derive a host of useful relations: 



k 
k' 



d 
d 

d ' 
1 



• /ci = ^(r -r]w ), 



cos 7 
sin 7 



2 2 

r — r^ti; 



^^^sin/5, (C.5) 



and 



dA; dk' d-f k k' 9{k) e{k') 9{-f) ^(27r - 7) = ^ dr dwd(3rw e{r) 0{w) 0{P) 0{2Tr - (3) . (C.6) 
In terms of these variables, (g', V , t' ,j\ Al^(A) \q, I, t,j)^^ takes the form 
{c[,l\t\3\M\K)\q,l,t,3)^i = -^j D'\ogr^W,,,,T,{r_)T,{r^) 

1=1 m=l; mT^i 

where 

D' = dxdx'drd«;d/?rw?7^(a;)^(l-a;)^(a;')^(a;-a;')^(r)^H^(/5)^(27r-/3) 
X (C.8) 

161 



To avoid roundoff error, we have defined simplified dimensionless versions of the 
differences of the free masses and the derivatives of these differences: 



r^?7(l — X — x') + w'^r]'^{l — x — x') — 2wr ■^{x\i — x\+ x'[l — x']) cos (5 



x{l — x)x'{l — x' 



r 1] + w (1 — ,x + 1 — x') + 2rw^{l — x + 1 — x') cos /3 
(1 -x)(l -x') ' 
r^?7 + w'^{x + x'Y — 2rWy/fi{x + x') cos P 



2 9 
A' — - A 

(1 — x'Y rj 

A',. = ^-4-. (C.9) 



Eq. ( C7 ) is now dimensionless, except for the factor in front of the integral, which 
is proportional to l/cP. 

As x' — s> X (r] — > 0), the contributions to the integrand from the first and second 
terms in the sum diverge. In the limit x' x, the contribution to the integral from 
these terms can be writtenQ 



h,,,, I Z}'logr/r,,(r)T,(r)e-=^2(Ad)-V^ 



1 V!^ 



2773^2 

X + + (C.IO) 

An examination of this integral reveals a problem: the transverse integrals are zero 
and the longitudinal integrals are infinite. To solve this problem, we consider what 
would have happened if we had not yet taken e ^ 0. In this case, the transverse 
integrals would be zero and the longitudinal integrals would be finite. Thus the 
e — s> limit of this integral would be zero. This means that this integral is actually 

zero, and we can subtract it from the full integral in Eq. (|C.7|) : 

^^We do not replace all the occurrences of x' with x because doing so hampers the convergence of 
the integral. 
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{q',l',t',j\M\A)\q,l,t,j)f^^ 



D'logrj 



X 



1=1 m=l; m^i 



X e-^2(Ac()-4«,4 ( ^ _ 64(Aci)-^^ ) Lf{x'){x + - x + 1 - 



(C.ll) 



Once we have performed this subtraction, there is no ambiguity about the value of 
the full integral, and it converges when computed numerically. 

C.1.2 Combination of the Integrals 

We now combine all the five-dimensional integrals into one integral: 



where 



27r3ci2 



D' 



B,F 



-^CON + -^EX + -^IN+EX + -^IN 



, (C.12) 



/con = [Sj,_25,',25g,2 + Sj,2S,',i5,,i-Sj,oS,'ML\f\x')Tt,{r.)^^^^ 



.(2) 



-2{Ad)-''AFKAiK 



) 



X 



r 



sri, + 



x{l-x) x'{l-x') x{l-x)x'{l-x') 



c(3) 



/iN+EX = W^,,,.4'^(x')ri,(r_)ri(r+)e 



-(Ad)-4A^ 



1 ^ g-2(A(i)-4A^KAjK^ -j^ 



X 



x(l — x)x'{l — x') r] ApK^iK 

- ri^w"^ - 2r]W{r'^ + 2w^[l - 2x]) - r/3(r^ + Ar'^w'^[l - 2x\ 



7-B,F 

'in 



+ Qw\l - 2xY) -Acos(3^rw{[-l+rfy + ■qw'^[l -x- x'Y){l - x - x') 

+ 8r^w^a;(l - a;) - 4 cos^/3r^w^(r;*^ + 2r;^[l - 2x\ + 4a;[l - a;] - 4?7^a;[l - x] 

- 2r7[l -2x\)- 4r?W(w^[l - 2xf + r^[l - 2x(l - x)\) 

+ 277(r^ - 2r'^w'^[l - 2x\ - 4w^x[-l + 3x - Ax^ + 2x^])] , 

5 5 

= logr/ l^,,,,T,,(r_)T,(r+)e-(^'^)-'(^--+^-)E^; 11 ^- 

i=l m=l; m^i 
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^-U{M)-^^]l'^\x'){x + x') 



X (1 - X + 1 - x') 



(C.13) 



We have rewritten the integrand of (g', Z', t', j\ A^^(A) |g, Z, t, j)iN+EX ehminate large 
roundoff errors, at the expense of making it more comphcated. To further avoid 



roundoff errors, we rewrite a few of the •S'^^g/'s: 



^3,3 

c(3) 
^3,4 

o(l) 
'-'4,4 



cos{j-i){2rf + 2r]{l - 2x) + (1 - 2xf), 

— sin7sin(j7)(2?7^[l — 2x\ + rj[l — Qx + 6x^] — 2x[l — 3x + 2x^]), 
cos(j7)(l-2x)(l + 2r/-2a;). (C.14) 



Note that to compute some of the trigonometric functions that appear in this integral 
[such as cos(j7)] in terms of the integration variables, it is necessary to use recursion 
relations that define these functions in terms of cos 7 and sin 7 so that we can use Eq. 



The integral in Eq. ( C.12 ) converges slowly. This is because it is strongly peaked 
when x' ~ x, even though we have subtracted the false divergence. To spread out 
this region, we change variables from x' to p where 



x' = x[l- e-n . 



(C.15) 



Now the integral is 

{q',l',t',j\M'{A)\q,l,t,j),_^ 



r L roo roo roo /'27r 

dx dp dr I dw d/Srwr]'^ 
Jo Jo Jo Jo Jo 

(C.16) 



27r3d2 



X Li^\x) 



-Z^CON + -^EX + -^IN+EX + luf 



As a final step, we note that VEGAS requires the region of integration to be finite. 
Thus we change variables from p, r, and w to i/p, i/r, and i/w: 
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and then the final expression for the contribution to the matrix elements from the 
five-dimensional integral is 

l\ t', j\ M^{k) \q, I, t, j)5_D = -^j^ dx j^^ dyp j^^ dyr dy^ dp 
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X ^^ (l + ^^)2(l + ^^)2(l+;j2 ^%^^H^) KC0N + /^X + WX + Z^N^J .(C.18) 

The integral converges nicely in this form. 
C.2 Miscellaneous Issues 

C.2.1 Rotational Symmetry: j — > —j 

To compute the spectrum, we compute the matrix (g', I', t', j\ A^^(A) \q, I, t, j) and 
diagonalize it for each value of j separately. For a given the matrices with j = \j\ 
and j — are simply related, and we can use this fact to avoid computing and 
diagonalizing both of them. By inspection, we determine that 

{q',l',t',-j\M\A)\q,l,t,-j) = E {q"',l',t',j\M\A)\q",l,t,j)f{q,q") 

q",q"' 

X (C.19) 



rF , 7- , 7-B,F 



where 

Q) = Sq,lSq',2 + Sq,2Sq',l + (5q,3(^q',3 " (^q,4(^q',4- (C.20) 



165 



This simply means that the basis \q, I, t, —j) is related to the basis \q, I, t,j) by swap- 
ping the states \l,l,t,j) and \2,l,t,j), and changing the sign of |4,/,t,j). Renaming 
basis states and changing their phases has no effect on the eigenvalues of the matrix; 
^(-j)"^(^p)\ has the same mass as |^-''*(P)). It also means that since 



so 



qlt 

Y: [Kt t,j) + |2, 1, t,j) + R^t |3, 1, t,j) + Kt |4, 1, t,j)] , (C.21) 



we have 



^(-.>(P)^ = E [Kt \2,l,t,3) + Ril \i,l,t,j) + Ril |3,M,j) - RZ |4,M,j) 

It 

q'qlt 

Thus by diagonalizing {q', I', t',j\ A4^(A) \q, I, t,j), we also obtain the eigenvalues and 
eigenstates of {q', I', t' , —j\ Al^(A) \q, I, t, —j). 

C.2.2 Reducing the Number of Matrix Elements to Compute 

There are a few facts that allow us to reduce the number of matrix elements that 
we have to compute. First, because of gluon-exchange symmetry, the basis state 
|g, /, is zero if / + j is even and g = 4, or if / + j is odd and q^ ^. Second, A1^(A) 
is Hermitian, and its matrix elements in this basis are real; so it is symmetric in this 
basis. Third, by inspection, we see that 

(l,r,t',j|Al2(A)|2,/,t,j) = 0. (C.23) 
Finally, there are some redundancies and additional zeros in the matrix when j = 0: 

(2,/',t',0| A^2(A) |2,/,t,0) = (l,/',t',0|M2(A) |l,/,t,0), 
(l,r,i',0| Al2(A) |3,Z,i,0) = (2, r, i', 0| Al^(A) |3, Z, 0) , 
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(l,/',t',0|A^2^A) |4,/,t,0) = -(2,/',t',0|A<2(A)|4,/,t,0), 
{3,l',t',0\M\A)\A,l,t,0) = 0. (C.24) 

C.2.2 Estimating Uncertainties in the Spectrum 

When we use the VEGAS Monte Carlo routine to compute the matrix elements, 
the results for the matrix elements have statistical uncertainties. In order to control 
the resulting uncertainties in the spectrum, we would like to have a method that 
allows us to estimate how accurately we must calculate any given matrix element in 
order for the spectrum to have a desired uncertainty. 

Suppose a diagonal matrix element (g, /, t, j| A^^(A) \q,l,t,j) is given by 

{q,l,t,j\M'{A)\q,l,t,j) = Z±5, (C.25) 

where Z is the Monte Carlo estimate of the matrix element and S is the associated 
absolute uncertainty. Using Mathematica with test matrices, it is straightforward to 
convince oneself that if 6 is small, then it will yield a relative uncertainty e ^ ^ 6/Z 
in the eigenvalues of the matrix^. This translates to a relative uncertainty of e^^^^ ~ 
6/{2Z) in the masses. 

Estimating the uncertainty in the spectrum due to uncertainties in off-diagonal 
matrix elements is more difficult. Using Mathematica with test matrices, the simplest 
method that we have found that is reasonably reliable is to use a type of degenerate 
perturbation theory. When we have an off-diagonal matrix element given by 

{q',l\t',j\M\A)\q,l,t,j) = Z±5, (C.26) 

^^In the development of this method, we are guided by the principles of quantum-mechanical 
perturbation theory, although we cannot legitimately use perturbation theory to analyze the 
uncertainties. 
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we diagonalize the two matrices 

/ {q,l,t,j\M^{A)\q,l,t,j) Z + 6 

[ Z + S {q',l',t',j\M\A)W,l',t',j) 

and 

/ {q,l,t,j\M\A)\q,l,t,j) Z-5 

\ Z-5 {q',l',t',j\M\A)W,l',t',j) 

and we compare their eigenvalues to the eigenvalues of the matrix 

/ {q,l,t,j\M\A)\q,l,t,j) Z 
[ Z {q',l',t',j\M'{A)\q',l',t',j) 

We then define e , to be the largest relative deviation that we have found in the 

n 

eigenvalues, and we estimate the resulting relative uncertainty in the mass spectrum to 
be Cj^^ ~ e 2/2. This estimate tends to work well unless there are too many diagonal 
matrix elements that are nearly degenerate with either {q,l,t, j\ J\4^{A) \q,l,t,j) or 
{q',l',t',j\M\A)\q',l',t',j). 

To achieve a relative uncertainty of 0{e) in the glueball masses, we require e,^^ < s 
for each matrix element. This method tends to work reasonably well. The main 
difficulty is that we have neglected to consider the highly nonlinear couplings between 
the uncertainties in different matrix elements. For this reason, as we increase the size 
of the matrix, the error in our estimate of the uncertainty eventually becomes critical. 
At this point, the spectrum that we get when we diagonalize the matrix becomes 
completely unreliable. (The evidence of the breakdown is sudden contamination of 
the low-lying wave functions with high-order components.) It should be possible 
to develop more sophisticated methods of estimating uncertainties to suppress this 
problem. 



(C.27) 



(C.28) 



(C.29) 
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